結果
問題 | No.109 N! mod M |
ユーザー | legosuke |
提出日時 | 2019-10-28 18:59:51 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,696 bytes |
コンパイル時間 | 1,585 ms |
コンパイル使用メモリ | 167,200 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-09-14 21:16:21 |
合計ジャッジ時間 | 8,464 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge6 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 126 ms
5,248 KB |
testcase_02 | AC | 118 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | TLE | - |
testcase_05 | -- | - |
testcase_06 | -- | - |
testcase_07 | -- | - |
testcase_08 | -- | - |
ソースコード
#include <bits/stdc++.h> #define endl '\n' #define int long long #define lint long long #define pii pair<int,int> #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() #define SZ(v) ((int)v.size()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) #define MINF(a) memset(a,0x3f,sizeof(a)) #define POW(n) (1LL<<(n)) #define POPCNT(n) (__builtin_popcount(n)) #define IN(i,a,b) (a <= i && i <= b) using namespace std; template <typename T> inline bool CHMIN(T& a,T b) { if(a>b) { a=b; return 1; } return 0; } template <typename T> inline bool CHMAX(T& a,T b) { if(a<b) { a=b; return 1; } return 0; } template <typename T> inline void SORT(T& a) { sort(ALL(a)); } template <typename T> inline void REV(T& a) { reverse(ALL(a)); } template <typename T> inline void UNI(T& a) { sort(ALL(a)); a.erase(unique(ALL(a)),a.end()); } template <typename T> inline T LB(vector<T>& v, T a) { return *lower_bound(ALL(v),a); } template <typename T> inline int LBP(vector<T>& v, T a) { return lower_bound(ALL(v),a) - v.begin(); } template <typename T> inline T UB(vector<T>& v, T a) { return *upper_bound(ALL(v),a); } template <typename T> inline int UBP(vector<T>& v, T a) { return upper_bound(ALL(v),a) - v.begin(); } template <typename T1, typename T2> ostream& operator<< (ostream& os, const pair<T1,T2>& p) { os << p.first << " " << p.second; return os; } template <typename T1, typename T2> istream& operator>> (istream& is, pair<T1,T2>& p) { is >> p.first >> p.second; return is; } template <typename T> ostream& operator<< (ostream& os, const vector<T>& v) { REP(i,v.size()) { if (i) os << " "; os << v[i]; } return os; } template <typename T> istream& operator>> (istream& is, vector<T>& v) { for(T& in : v) is >> in; return is; } template <typename T = int> vector<T> make_v(size_t a) { return vector<T>(a); } template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) { return vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...)); } template <typename T, typename V> typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v) { t = v; } template <typename T, typename V> typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v) { for(auto &e : t) fill_v(e,v); } const lint MOD = 1000000007; const lint INF = 0x3f3f3f3f3f3f3f3f; const double EPS = 1e-10; bool is_prime(int n) { if (n < 2) return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0) return false; } return true; } // x = a^(-1) // a と m が互いに素 long long mod_inv(long long a, long long m) { long long b = m, x = 1, y = 0, t; while (b > 0) { t = a / b; swap(a -= t * b, b); swap(x -= t * y, y); } x %= m; if (x < 0) x += m; return x; } void solve(int n, int m) { if (m == 1 || n >= m) { cout << 0 << endl; } else if (is_prime(m)) { int ans = m - 1; for (int i = m - 1; i >= n + 1; --i) { (ans *= mod_inv(i, m)) %= m; } cout << ans << endl; } else { int ans = 1; FOR(i, 1, n + 1) { (ans *= i) %= m; } cout << ans << endl; } } void _main() { int T; cin >> T; REP(i, T) { int N, M; cin >> N >> M; solve(N, M); } } signed main(signed argc, char **argv) { if (argc > 1) { if (strchr(argv[1], 'i')) freopen("input.txt", "r", stdin); if (strchr(argv[1], 'o')) freopen("output.txt", "w", stdout); } cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(10); _main(); return 0; }