import sys
input=lambda: sys.stdin.readline().rstrip()
mod=10**9+7
n_max=2*(10**6+1)
F,FI=[0]*(n_max+1),[0]*(n_max+1)
F[0],FI[0]=1,1
for i in range(n_max):
  F[i+1]=(F[i]*(i+1))%mod
FI[n_max-1]=pow(F[n_max-1],mod-2,mod)
for i in reversed(range(n_max-1)):
  FI[i]=(FI[i+1]*(i+1))%mod
def comb(x,y):
  return (F[x]*FI[x-y]*FI[y])%mod
t=int(input())
N=set(str(i) for i in range(10))
for _ in range(t):
  S=input()
  a,b="",""
  st=0
  for s in S:
    if s==",":
      st=1
    if s in N:
      if st==0:
        a+=s
      else:
        b+=s
  a,b=int(a),int(b)
  if S[0]=="C":
    if a<b:
      print(0)
    else:
      print(comb(a,b))
  elif S[0]=="P":
    if a<b:
      print(0)
    else:
      print((F[a]*FI[a-b])%mod)
  else:
    if a==0 and b==0:
      print(1)
    else:
      print(comb(a+b-1,a-1))