#include <cstdio>

using i32 = int;
using u32 = unsigned int;
using u64 = unsigned long long;

struct Mod64
{
    __attribute__((always_inline)) Mod64() : n_(0x0) {}

    __attribute__((always_inline)) Mod64(u32 n) : n_(init(n)) {}

    __attribute__((always_inline)) static u32 modulus()
    {
        return mod;
    }

    __attribute__((always_inline)) static u32 init(u32 w)
    {
        return reduce(u64(w) * r2);
    }

    __attribute__((always_inline)) static void set_mod(u32 m)
    {
        mod = m;
        inv = m; for (i32 i = 0x0; i < 0x5; ++i) inv *= 0x2 - inv * m;
        r2 = -u64(m) % m;
    }

    __attribute__((always_inline)) static u32 reduce(u64 x)
    {
        u32 y = u32(x >> 0x20) - u32((u64(u32(x) * inv) * mod) >> 0x20);

        return i32(y) < 0x0 ? y + mod : y;
    }

    __attribute__((always_inline)) Mod64& operator += (Mod64 rhs)
    {
        n_ += rhs.n_ - mod;

        if(i32(n_) < 0) n_ += mod;

        return *this;
    }

    __attribute__((always_inline)) Mod64 operator + (Mod64 rhs) const
    {
        return Mod64(*this) += rhs;
    }

    __attribute__((always_inline)) Mod64& operator *= (Mod64 rhs)
    {
        n_ = reduce(u64(n_) * rhs.n_);

        return *this;
    }

    __attribute__((always_inline)) Mod64 operator * (Mod64 rhs) const
    {
        return Mod64(*this) *= rhs;
    }

    __attribute__((always_inline)) u32 get() const
    {
        return reduce(n_);
    }

    static u32 mod, inv, r2;

    u32 n_;
};

u32 Mod64 :: mod, Mod64 :: inv, Mod64 :: r2;

__attribute__((always_inline)) i32 expo(i32 A, i32 B)
{
    i32 result = 1;

    Mod64 rA(A);

    while(B)
    {
        if(B & 1)

            result = (result * A) % 10007;

        rA *= rA;

        A = rA.get();

        B >>= 1;
    }

    return result;
}

#define SIZE 0x10000

__attribute__((always_inline)) void read(int &num)
{
    static char inBuffer[SIZE];

    static int p = ~-SIZE; num = 0x0;

    while(inBuffer[p] < 0x30 | inBuffer[p] > 0x39)
    {
        ++p != SIZE || (fread(inBuffer, 0x1, SIZE, stdin), p = 0x0);
    }

    while(inBuffer[p] > 0x2F & inBuffer[p] < 0x3A)
    {
        num = num * 0xA + inBuffer[p] - 0x30;

        ++p != SIZE || (fread(inBuffer, 0x1, SIZE, stdin), p = 0x0);
    }
}

char outBuffer[0x927C0];

__attribute__((always_inline)) void itoa(i32 x)
{
    static i32 p = 0x0;

    i32 digits = x > 0x1869F    ? 0x6 :
                 x > 0x270F     ? 0x5 :
                 x > 0x3E7      ? 0x4 :
                 x > 0x63       ? 0x3 :
                 x > 0x9        ? 0x2 : 0x1;

    for(i32 i = ~-digits; ~i; --i)
    {
        outBuffer[p + i] = x % 0xA + 0x30;

        x = x / 0xA;
    }

    p = p + digits; outBuffer[p++] = 0xA;
}

i32 main()
{
    freopen("matrice5.in", "r", stdin);
    freopen("matrice5.out", "w", stdout);

    i32 T, N, M, P, K;

    read(T);

    Mod64 :: set_mod(10007);

    do
    {
        read(N); read(M); read(P); read(K);

        itoa((expo(P * K, (M - 1) * (N - 1)) * expo(P, M + N - 1)) % 10007);

    }while(--T);

    puts(outBuffer);

    return 0;
}