#include #include #include template struct comb_util { int sz; std::vector mfact, minv_fact; comb_util(int N) : sz(N + 1), mfact(sz), minv_fact(sz) { mfact[0] = 1; for (int i = 1; i <= sz; i++) mfact[i] = mfact[i - 1] * i % MOD; minv_fact[sz] = inv(mfact[sz]); for (int i = sz - 1; i >= 0; i--) minv_fact[i] = minv_fact[i + 1] * (i + 1) % MOD; } // res = first * p^second std::pair fact_willson(Int n) const { Int e = 0; if (n <= sz) return std::make_pair(mfact[n], 0); Int res; std::tie(res, e) = fact_willson(n / MOD); e += n / MOD; if ((n / MOD) % 2 != 0) res = MOD - res; return std::make_pair(res * mfact[n % MOD] % MOD, e); } Int fact(Int n) const { return mfact[n]; } Int inv(Int n) const { return pow(n, MOD - 2); } Int pow(Int n, Int a) const { Int res = 1, exp = n; for(; a; a /= 2) { if (a & 1) res = res * exp % MOD; exp = exp * exp % MOD; } return res; } Int perm(Int n, Int r) { if (r < 0 || n < r) return 0; else return mfact[n] * minv_fact[n - r] % MOD; } Int binom_lucus(Int n, Int r) const { if (n < 0 || r < 0 || n < r) return 0; assert(n <= sz); if (n >= MOD) return binom(n % MOD, r % MOD) * binom(n / MOD, r / MOD); else return r > n ? 0 : mfact[n] * minv_fact[n - r] * minv_fact[r]; } Int binom(Int n, Int r) const { if (n < 0 || r < 0 || n < r) return 0; return mfact[n] * minv_fact[r] % MOD * minv_fact[n - r] % MOD; } Int homo(Int n, Int r) const { if (n < 0 || r < 0) return 0; return r == 0 ? 1 : binom(n + r - 1, r); } }; using comb = comb_util; #include int main(){ int T; comb cm(2000010); scanf("%d\n", &T); for(int i=0;i