#include <iostream>
#include <fstream>
#define MOD 10007
using namespace std;
ifstream f ("matrice5.in");
ofstream g ("matrice5.out");
long long t,n,m,p,k;
long long putere(long long b,long long p)
{
    long long rez = 1;
    while(p)
    {
        if(p%2 == 1)
            rez = (rez * b)%MOD;
        b = (b*b)%MOD;
        p/=2;
    }
    return rez;
}
int main()
{
    f >> t ;
    for(int i= 1; i<=t; i++)
    {
        f >> n >> m>> p >> k;
        long long rez;
        /// in primul rand umplem matricea pe toate pozitiile in afara de ultima si linie si ultima coloane cu toate numerele
        rez = putere(k*p,(n-1)*(m-1));
        /// apoi umplem ce a ramas cu p, deoarece avand k*p numere, p numere vor putea transforma suma de pe linie/coloana intr o suma valida
        rez = rez * putere (p,n+m-1);
        rez%=MOD;
        g << rez<<'\n';
    }
}