#include <bits/stdc++.h>

using namespace std;

ifstream in ("matrice5.in");
ofstream out ("matrice5.out");

const int M=10007;
int E(int a, int b)
{
    int r=1;
    while (b)
    {
        if (b&1)
            r=r*a%M;
        a=a*a%M;
        b>>=1;
    }
    return r;
}
int main()
{
    int T;
    in>>T;
    while (T--)
    {
        int n,m,p,k;
        in>>n>>m>>p>>k;
        out<<E(p*k%M,(n-1)*(m-1))*E(p,n+m-1)%M<<'\n';
    }
    return 0;
}