#include <cstdio>

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

u32 T, N, M, P, K;

u32 expo(u32 A, u32 B)
{
    u32 result = 1;

    while(B)
    {
        if(B & 1)

            result = ((u64)result * A) % 10007;

        A = ((u64)A * A) % 10007;

        B >>= 1;
    }

    return result;
}

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

    scanf("%d", &T);

    while(T--)
    {
        scanf("%d %d %d %d", &N, &M, &P, &K);

        printf("%d\n", ((u64)expo(P * K, (M - 1) * (N - 1)) * expo(P, M + N - 1)) % 10007);
    }

    return 0;
}