#include <iostream>
#include <fstream>
using namespace std;
const int mod=10007;
int t,n,m,p,k;
int expo(int A,int B)
{
    int ret=1,p2=A;
    for(int pu=0;(1<<pu)<=B;pu++)
    {
        if(((1<<pu)&B))
            ret=(ret*p2)%mod;
        p2=(p2*p2)%mod;
    }
    return ret;
}
int main()
{
    ifstream f("matrice5.in");
    ofstream g("matrice5.out");
    f>>t;
    for(int cnt=1;cnt<=t;cnt++)
    {
        f>>n>>m>>p>>k;
        g<<(expo(p*k,(n-1)*(m-1))*expo(p,n+m-1))%mod<<'\n';
    }
    return 0;
}