#include<cstdio>
using namespace std;
const int MOD=1e9+7;
inline int qp(int b,int f){int res=1;for(;f;f>>=1,b=b*1ll*b%MOD) if(f&1) res=res*1ll*b%MOD;return res;}
int L,n,K,c;
struct Matrix{
	int a[2][2];
	Matrix operator*(const Matrix&b){
		Matrix res;
		for(int i=0;i<2;++i) for(int j=0;j<2;++j)
			for(int k=res.a[i][j]=0;k<2;++k) res.a[i][j]=(res.a[i][j]+a[i][k]*1ll*b.a[k][j])%MOD;
		return res;
	}
}A,B;
int main(){
	for(scanf("%d",&L);L--;printf("%lld\n",(c*(c-1ll)%MOD*B.a[0][0]+c*(c-1ll)%MOD*(c-2)%MOD*B.a[1][0])%MOD
	*qp((c-1+(c-2ll)*(c-2))%MOD,K*(n-1ll)%(MOD-1))%MOD)){
		scanf("%d%d%d",&n,&K,&c);
		A.a[0][0]=0;A.a[0][1]=c-1;A.a[1][0]=1;A.a[1][1]=c-2;
		B.a[0][0]=1;B.a[0][1]=0;B.a[1][0]=0;B.a[1][1]=1;
		for(int i=K-2;i;i>>=1,A=A*A) if(i&1) B=B*A;
	}
	return 0;
}