#include <fstream>

using namespace std;

ifstream cin("matrice5.in");
ofstream cout("matrice5.out");

int v[1005];

long long exponent(long long x, long long n)
{
    if(n == 0)
        return 1;
    long long a = exponent(x, n/2) % 10007;
    if(n % 2 == 0)
        return (a * a) % 10007;
    else
        return (((a * a) % 10007) * x) % 10007;
}

int main()
{
    long long n,m,p,k,t,i;
    cin >> t;
    for(i = 1; i <= t; i++)
    {
        cin >> n >> m >> p >> k;
        cout << exponent(k*p % 10007, (n-1)*(m-1) % 10007) * exponent(p, n+m-1) % 10007 << "\n";
    }
    return 0;
}