#include <fstream>

const int MOD = 10007;

using namespace std;
ifstream f("matrice5.in");
ofstream g("matrice5.out");

int T,N,M,P,K,X;
long long sol,s1,s2;

long long pow(int x, int p)
{
    long long rez;
    if (p == 1)
        return x;
    rez = pow(x,p/2);
    rez = (rez*rez) % MOD;
    if (p % 2 == 1)
        rez = (rez*x)%MOD;
    return rez;
}

int main()
{
    f >> T;
    while (T--)
    {
        f >> N >> M >> P >> K;
        s1 = pow(P*K, (N-1)*(M-1) );
        s2 = pow(P,N+M-1);
        sol = s1*s2 % MOD;
        g << sol << '\n';
    }

    f.close();
    g.close();
    return 0;
}