結果

問題 No.2120 場合の数の下８桁
コンテスト
ユーザー 259_Momone
提出日時 2022-11-04 22:28:25
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 180 ms / 2,000 ms
コード長 1,361 bytes
コンパイル時間 2,281 ms
コンパイル使用メモリ 246,860 KB
最終ジャッジ日時 2025-02-08 18:03:30
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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;
}
0