#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=(sol*n)%MOD;
        n=(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;
        g<<power(k,(n-1)*(m-1))*power(p,n*m)%MOD<<'\n';
    }
    return 0;
}