#include <bits/stdc++.h>
using namespace std;
const int mod=10007;
ifstream f("matrice5.in");
ofstream g("matrice5.out");
int pow(int a, int b) {
    long long n,p,x=1;
    n=a;p=b;
    while (p) {
        if (p%2==0) {
            p=p/2;
            n=(n*n)%mod;
        }
        else {
            p--;
            x=(x*n)%mod;
        }
    }
    return x;
}
int main() {
    int t, i, n, m, k, p;
    f>>t;
    for(i=1; i<=t; i++) {
        f>>n>>m>>p>>k;
        g<<((pow(p*k,(n-1)*(m-1))%mod)*(pow(p,n+m-1)%mod))%mod<<'\n';
    }
    return 0;
}