#include #include #include #include #define REP(i,a,b) for(i=a;i0){ q=r0/r1; r2=r0%r1; a2=a0-q*a1; b2=b0-q*b1; r0=r1; r1=r2; a0=a1; a1=a2; b0=b1; b1=b2; } *c=r0; *a=a0; *b=b0; } ll get_inv(ll n, ll p){ ll a,b,c; extended_euclid(n,p,&c,&a,&b); if(a a) return 0; return (fact[a] * fact_inv[a-b])%MOD; } int main(){ int i, j; int N, M, K; int T; ll f[MAX], ff[MAX], res; ll C[MAX][MAX]; REP(i,1,MAX) inv[i] = get_inv(i, MOD); fact[0] = fact_inv[0] = 1; REP(i,1,MAX) fact[i] = (fact[i-1] * i)%MOD; REP(i,1,MAX) fact_inv[i] = (fact_inv[i-1] * inv[i])%MOD; rep(i,MAX) C[0][i] = 0; rep(i,MAX) C[i][0] = 1; REP(i,1,MAX) REP(j,1,MAX) C[i][j] = (C[i-1][j-1] + C[i-1][j])%MOD; scanf("%d",&T); while(T--){ scanf("%d%d%d",&N,&M,&K); REP(i,1,N+1){ if(K < M) ff[i] = (P(i, K) * pw(i-K, M-K, MOD))%MOD; else ff[i] = P(i, M); REP(j,1,i) ff[i] = (ff[i] - ff[j]*C[i][j]) % MOD; if(ff[i] < 0) ff[i] += MOD; f[i] = (ff[i] * fact_inv[i])%MOD; } res = 0; REP(i,1,N+1) res += f[i]; res %= MOD; printf("%d\n",(int)res); } return 0; }