#include <bits/stdc++.h>

using namespace std;
const int modulo = 10007;

int putere(int a,int b)
{
    int rez = 1;
    while(b != 0)
    {
        if(b % 2 == 1)
            rez = (rez*a)%modulo;
        a = (a*a)%modulo;
        b /= 2;
    }
    return rez;
}

int main()
{
    freopen("matrice5.in","r",stdin);
    freopen("matrice5.out","w",stdout);
    int test;
    cin >> test;
    while(test != 0)
    {
        int n,m,p,k;
        cin >> n >> m >> p >> k;
        int ans = putere((k*p)%modulo, (n-1)*(m-1));
        ans = (ans* putere(p, n+m-1)) % modulo;
        cout << ans << "\n";
        test --;
    }

    return 0;
}