#include typedef long long ll; typedef unsigned long long ull; #define FOR(i,a,b) for(int (i)=(a);i<(b);i++) #define REP(i,n) FOR(i,0,n) #define RANGE(vec) (vec).begin(),(vec).end() using namespace std; std::vector primeFactor(int n) { std::vector res; for (int i = 2; i*i <= n; ++i) { if (n%i == 0) { res.push_back(i); while (n%i==0) n /= i; } } if (n!= 1) res.push_back(n); return std::move(res); } template T extgcd(T a, T b, T &x, T &y) { T d = a; if (b != 0) { d = extgcd(b, a%b, y, x); y -= (a/b)*x; } else { x=1; y=0; } return d; } // mod M の逆元を求める ll mod_inverse(ll a, int M) { ll x, y; ll d = extgcd(a, (ll)M, x, y); assert(d == 1); return (M + x%M) % M; } class NFactorialModM { public: void solve(void) { int T; cin>>T; // // N >= M なら N*(N-1)*...*(M+1)*M*(M-1)*...*1 = 0 mod M // N < M としてよい ~~ // // M-N <= 10^5 より // (M-1)!/N! = (M-1)*...*(N+1) は O(10^5) で計算可能 // // N! = N!/(M-1)! * (M-1)! = (M-1)! * ((M-1)!/N!)^(-1) であることに留意する。 // // * M が素数のときウィルソンの定理より // (M-1)! == (M-1) mod M // // * M が合成数のとき M = (M/p) * p (pは素数) とかける // * M/p <= N ならば N! で M が現れるので 0 // * M/p > N ならば M > p*N >= 2*N なので 10^5 > M-N > N (N! 計算すればよい) // REP(t,T) { int N,M; cin>>N>>M; if (N >= M || M == 1) { cout<<0< M) // M is prime { ll a = 1; for (int i = 1; i < M-N; ++i) (a *= (N+i)) %= M; cout<<(mod_inverse(a,M)*(M-1))%M<solve(); delete obj; return 0; } #endif