結果
問題 | No.2120 場合の数の下8桁 |
ユーザー | 259_Momone |
提出日時 | 2022-11-04 22:28:25 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 205 ms / 2,000 ms |
コード長 | 1,361 bytes |
コンパイル時間 | 3,983 ms |
コンパイル使用メモリ | 259,132 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-18 20:24:32 |
合計ジャッジ時間 | 4,685 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <bits/extc++.h> int main() { using namespace std; constexpr unsigned long mod{100000000}; unsigned long M, N; cin >> M >> N; unsigned long value_p{1}, ord_2_p{}, ord_5_p{}; unsigned long value_m{1}, ord_2_m{}, ord_5_m{}; if(N > M){ printf("%08lu\n", 0UL); return 0; } const auto multiplies{[](unsigned long& value, unsigned long& ord_2, unsigned long& ord_5, unsigned long x){ while(x && x % 2 == 0){ ++ord_2; x /= 2; } while(x && x % 5 == 0){ ++ord_5; x /= 5; } (value *= x) %= mod; }}; for(unsigned long i{}; i < N; ++i){ multiplies(value_p, ord_2_p, ord_5_p, M - i); multiplies(value_m, ord_2_m, ord_5_m, i + 1); } assert(ord_2_p >= ord_2_m); assert(ord_5_p >= ord_5_m); auto ord_2{ord_2_p - ord_2_m}; auto ord_5{ord_5_p - ord_5_m}; auto value{[](unsigned long a, unsigned long n, unsigned long b){ while(n){ if(n & 1)(b *= a) %= mod; (a *= a) %= mod; n /= 2; } return b; }(value_m, 4999999, value_p)}; while(ord_2 > 0){ --ord_2; (value *= 2) %= mod; } while(ord_5 > 0){ --ord_5; (value *= 5) %= mod; } printf("%08lu\n", value); return 0; }