#include <bits/stdc++.h>
#define MOD 10007
using namespace std;
int n,m,p,k,t,p1,p2;
int power(int n,int p)
{
    int sol=1;
    if(p==0)return 1;
    for(int i=0; i<32; ++i)
    {
        if(p&(1<<i))sol=(1LL*sol*n)%MOD;
        n=(1LL*n*n)%MOD;
    }
    return sol;
}
int main()
{
    ifstream f ("matrice5.in");
    ofstream g ("matrice5.out");
    f>>t;
    ++t;
    while(--t)
    {
        f>>n>>m>>p>>k;
        p1=(n-1)*(m-1);
        p2=n*m;
        g<<(1LL*power(k,p1)*power(p,p2))%MOD<<'\n';
    }
    return 0;
}