#include<bits/stdc++.h>
#define all(vec) vec.begin(),vec.end()
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
const ll MOD=1000000007;
const ll INF=1000000010;
const ll LINF=4000000000000000010LL;
const int MAX=310;
const double EPS=1e-9;
int dx[4]={0,1,0,-1};
int dy[4]={1,0,-1,0};
ll f[1000010],fi[1000010];
ll mod_pow(ll x,ll n,ll mod){
    ll ret=1;
    while(n>0){
        if(n&1){
            ret=ret*x%mod;
        }
        x=x*x%mod;
        n>>=1;
    }
    return ret;
}

void fact(ll n){
	f[0]=1;
	for(ll i=1;i<=n;i++){
		f[i]=f[i-1]*i%MOD;
	}
	fi[n]=mod_pow(f[n],MOD-2,MOD);
	for(ll i=n-1;i>=0;i--){		
		fi[i]=fi[i+1]*(i+1)%MOD;
	}
}

ll nck(ll n,ll k){
	if(k==0){
		return 1;
	}else{
		return f[n]*fi[k]%MOD*fi[n-k]%MOD;
	}
}

ll hck(ll n,ll k){
    return f[n+k-1]*fi[k]%MOD*fi[n-1]%MOD;   
}
int main(){
	int q;scanf("%d",&q);
	fact(1000000);
	while(q--){
		char c;int n,k;
		scanf("%*c%c(%d,%d)",&c,&n,&k);
		if(c=='C'){
			cout<<nck(n,k)<<endl;
		}else if(c=='P'){
			cout<<nck(n,k)*f[k]%MOD<<endl;
		}else if(c=='H'){
			cout<<hck(n,k)<<endl;
		}
	}
	return 0;
}