結果
問題 | No.2120 場合の数の下8桁 |
ユーザー | boatmuscles |
提出日時 | 2022-11-04 23:06:04 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 155 ms / 2,000 ms |
コード長 | 1,948 bytes |
コンパイル時間 | 1,355 ms |
コンパイル使用メモリ | 143,088 KB |
実行使用メモリ | 46,648 KB |
最終ジャッジ日時 | 2024-07-18 20:59:45 |
合計ジャッジ時間 | 3,152 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 43 ms
8,332 KB |
testcase_01 | AC | 43 ms
8,464 KB |
testcase_02 | AC | 44 ms
8,332 KB |
testcase_03 | AC | 43 ms
8,332 KB |
testcase_04 | AC | 42 ms
8,208 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 1 ms
5,376 KB |
testcase_07 | AC | 1 ms
5,376 KB |
testcase_08 | AC | 1 ms
5,376 KB |
testcase_09 | AC | 44 ms
8,204 KB |
testcase_10 | AC | 88 ms
28,088 KB |
testcase_11 | AC | 155 ms
46,520 KB |
testcase_12 | AC | 48 ms
8,204 KB |
testcase_13 | AC | 46 ms
8,332 KB |
testcase_14 | AC | 48 ms
8,328 KB |
testcase_15 | AC | 93 ms
28,084 KB |
testcase_16 | AC | 145 ms
46,648 KB |
testcase_17 | AC | 147 ms
46,648 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 145 ms
46,524 KB |
ソースコード
#include<cstdio> #include<vector> #include<map> #include<bitset> #include<atcoder/modint> #include<ext/pb_ds/assoc_container.hpp> #include<ext/pb_ds/tree_policy.hpp> using namespace __gnu_pbds; using namespace std; using namespace atcoder; using mint = modint; #define rep(i, n) for (int i = 0; i < (int)(n); ++i) #define rrep(i, n) for (int i = (int)(n)-1; i >= 0; --i) #define rep2(i, a, b) for (int i = (int)a; i < (int)(b); ++i) #define rrep2(i, a, b) for (int i = (int)(b)-1; i >= (int)(a); --i) constexpr int MAX_N = 10000000; bitset<MAX_N+1> is_composite;//prime_divisor[i]はiを割り切る素数の1つ vector<int> primes = {2}; int main(){ int M, N; scanf("%d\n%d", &M, &N); if(M < N){ printf("00000000\n"); return 0; }else if(M == N){ printf("00000001\n"); return 0; } for (int i = 4; i <= MAX_N; i += 2){ is_composite.set(i); } int _i = 3; for (int twoi; _i*_i <= MAX_N; _i += 2){ if(!is_composite[_i]){ primes.push_back(_i); twoi = _i<<1; for (int j = _i*_i; j <= MAX_N; j += twoi){ is_composite.set(j); } } } for (; _i <= MAX_N; _i += 2) if(!is_composite[_i]) primes.push_back(_i); gp_hash_table<int, int> prime_factorization;//注目している数を素因数分解したときの<素因数, 割り切る回数> auto solve = [&](int x, bool add){ int __x = x; for(auto &i : primes){ if(__x < i) break; int cnt = 0, _x = __x; while((_x /= i)){ cnt += _x; } prime_factorization[i] += add ? cnt : -cnt; } }; solve(M, true); solve(N, false); solve(M-N, false); mint::set_mod(100000000); mint answer{1}; for(auto [num, pow] : prime_factorization) answer *= mint::raw(num).pow(pow); printf("%08u\n", answer.val()); return 0; }