#include <iostream>
#include <fstream>
#define MOD 10007
using namespace std;
ifstream f ("matrice5.in");
ofstream g ("matrice5.out");
int t,n,m,p,k;

int putere(long long x, int p) {
    long long rez = 1;
    for (; p; p >>= 1) {
        if (p & 1)
            rez = (rez * x) % MOD;
        x = (x * x) % MOD;
    }
    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';
    }
}