#include <cstdio>

using i32 = int;
using i64 = long long;

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

    while(B)
    {
        if(B & 1)

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

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

        B >>= 1;
    }

    return result;
}

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

    i32 T; scanf("%d", &T);

    while(T--)
    {
        i32, N, M, P, K;

        scanf("%d %d %d %d", &N, &M, &P, &K);

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

    return 0;
}