mo,mx=[1000000007,2000000]

I,F,R=[[0 for i in range(mx+1)]for i in range(3)]
I[1]=F[0]=R[0]=1
for i in range(2,mx+1):
    I[i]=(mo-mo/i*I[mo%i]%mo)%mo
for i in range(1,mx+1):
    F[i]=(F[i-1]*i)%mo;
    R[i]=(R[i-1]*I[i])%mo;

n=input()
s=[raw_input() for i in range(n)]

n1,n2=[0,0]
for i in range(n):
    n1,n2=[int(s[i][2:s[i].index(",")]),int(s[i][s[i].index(",")+1:s[i].index(")")])]
    if s[i][0]=="C":
        print 0 if n1<n2 else (F[n1]*R[n1-n2]%mo)*R[n2]%mo
    elif s[i][0]=="P":
        print 0 if n1<n2 else F[n1]*R[n1-n2]%mo
    elif s[i][0]=="H":
        print 1 if n1==n2 else (F[n1+n2-1]*R[n1-1]%mo)*R[n2]%mo