#include <bits/stdc++.h>
using namespace std;
using ll = long long;
template<class T=ll> using vec = vector<T>;
#define BE(x) x.begin(), x.end()

const ll mod   = 1000000007;
const ll max_n = 10000000;
vec<ll> inv(max_n+1,1), fact(max_n+1,1), finv(max_n+1,1);
bool pre_f = true;

void comb_pre () {
    for (ll i = 2; i <= max_n; i++) {
        inv [i] = mod-inv[mod%i]*(mod/i)%mod; // inverse for i
        fact[i] = fact[i-1]*i%mod;            // i!
        finv[i] = finv[i-1]*inv[i]%mod;       // inverse for i!
    }
    pre_f = false;
}

ll comb (ll n, ll k) {
    if (n<k || n<0 || k<0) return 0;
    if (pre_f) comb_pre();
    return fact[n]*(finv[k]*finv[n-k]%mod)%mod;
}
    // set max_n and mod value
    // n and k are up to about 10^7
    // pre:O(max_n), query:O(1)


int main() {
    int q;
    cin >> q;
    for (int _ = 0; _ < q; _++) {

        string s;
        cin >> s;
        char c = s[0];
        string sn, sk;
        s.erase(0,2);
        bool f = false;
        while (true) {
            if (!f) {
                if (s[0] == ',') f = true;
                else sn += s[0];
                s.erase(0,1);
            }
            else {
                if (s[0] == ')') break;
                sk += s[0];
                s.erase(0,1);
            }
        }
        ll n = stoll(sn), k = stoll(sk);

        ll ans;
        if (c == 'C') ans = comb(n,k);
        if (c == 'P') ans = comb(n,k)*fact[k]%mod;
        if (c == 'H') ans = comb(n+k-1,k);
        cout << ans << endl;
    }
}