結果
問題 | No.2120 場合の数の下8桁 |
ユーザー | NAMIDAIRO |
提出日時 | 2022-11-04 23:40:20 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 953 ms / 2,000 ms |
コード長 | 1,550 bytes |
コンパイル時間 | 940 ms |
コンパイル使用メモリ | 103,184 KB |
実行使用メモリ | 81,536 KB |
最終ジャッジ日時 | 2024-07-18 21:25:44 |
合計ジャッジ時間 | 4,195 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 1 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 145 ms
81,528 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 3 ms
5,376 KB |
testcase_10 | AC | 51 ms
41,700 KB |
testcase_11 | AC | 140 ms
81,536 KB |
testcase_12 | AC | 2 ms
5,376 KB |
testcase_13 | AC | 7 ms
5,376 KB |
testcase_14 | AC | 12 ms
5,376 KB |
testcase_15 | AC | 363 ms
41,116 KB |
testcase_16 | AC | 131 ms
77,568 KB |
testcase_17 | AC | 920 ms
80,256 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 953 ms
81,328 KB |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <map> #include <queue> #include <set> #include <random> #include <iomanip> using namespace std; typedef long long ll; #define rep(i, n) for(int i = 0; i < (n); i++) template<class T> using vi = vector<T>; template<class T> using vii = vector<vi<T>>; template<class T> using viii = vector<vii<T>>; void chmax(ll & x, ll y) { x = max(x, y); } void pf(int n, vi<int>& prime, vi<int>& cnt, bool plus = true) { if (n == 1) return; if (plus) cnt[prime[n]]++; else cnt[prime[n]]--; pf(n / prime[n], prime, cnt, plus); return; } ll modpow(ll x, ll n, ll mod) { if (n == 0) return 1LL; ll res = modpow(x, n / 2, mod); res *= res; res %= mod; if (n % 2) res *= x; res %= mod; return res; } int main() { int n, m; cin >> n >> m; if(n < m){ printf("%08d\n", 0); return 0; } vi<int> prime(n + 1); prime[1] = 1; for (int i = 2; i <= n; i++) { if (prime[i]) continue; prime[i] = i; for (ll j = (ll)i * i; j <= n; j += i) { if (prime[j]) continue; prime[j] = i; } } prime[1] = 1; vi<int> cnt(n + 1); m = min(n - m, m); rep(i, m) { pf(n - i, prime, cnt); pf(m - i, prime, cnt, false); } ll ans = 1; const ll mod = 1e8; for (int i = 2; i <= n; i++) { if (cnt[i] == 0) continue; ans *= modpow(i, cnt[i], mod); ans %= mod; } printf("%08d\n", ans); return 0; }