結果

問題 No.2318 Phys Bone Maker
ユーザー tatyamtatyam
提出日時 2023-03-15 03:14:12
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,460 bytes
コンパイル時間 991 ms
コンパイル使用メモリ 100,284 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-18 12:59:08
合計ジャッジ時間 6,003 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 TLE -
testcase_03 AC 5 ms
6,944 KB
testcase_04 AC 7 ms
6,944 KB
testcase_05 AC 7 ms
6,948 KB
testcase_06 AC 6 ms
6,940 KB
testcase_07 AC 3 ms
6,944 KB
testcase_08 AC 7 ms
6,944 KB
testcase_09 AC 6 ms
6,944 KB
testcase_10 AC 10 ms
6,948 KB
testcase_11 AC 7 ms
6,944 KB
testcase_12 AC 10 ms
6,944 KB
testcase_13 AC 10 ms
6,940 KB
testcase_14 AC 7 ms
6,944 KB
testcase_15 AC 7 ms
6,944 KB
testcase_16 AC 7 ms
6,940 KB
testcase_17 AC 9 ms
6,944 KB
testcase_18 AC 10 ms
6,940 KB
testcase_19 AC 4 ms
6,944 KB
testcase_20 AC 8 ms
6,944 KB
testcase_21 AC 7 ms
6,940 KB
testcase_22 AC 6 ms
6,944 KB
testcase_23 AC 6 ms
6,944 KB
testcase_24 AC 9 ms
6,944 KB
testcase_25 AC 8 ms
6,940 KB
testcase_26 AC 8 ms
6,944 KB
testcase_27 AC 9 ms
6,944 KB
testcase_28 AC 5 ms
6,944 KB
testcase_29 AC 5 ms
6,944 KB
testcase_30 AC 4 ms
6,944 KB
testcase_31 AC 10 ms
6,940 KB
testcase_32 AC 8 ms
6,940 KB
testcase_33 AC 1 ms
6,948 KB
testcase_34 AC 8 ms
6,940 KB
testcase_35 AC 355 ms
6,940 KB
testcase_36 AC 1,885 ms
6,944 KB
testcase_37 AC 1,893 ms
6,944 KB
testcase_38 AC 1,803 ms
6,944 KB
testcase_39 AC 2,966 ms
6,944 KB
testcase_40 AC 2,902 ms
6,944 KB
testcase_41 TLE -
testcase_42 TLE -
testcase_43 AC 12 ms
6,944 KB
testcase_44 AC 255 ms
6,940 KB
testcase_45 AC 238 ms
6,944 KB
testcase_46 TLE -
testcase_47 AC 11 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdint>
#include <iostream>
#include <vector>
#include <atcoder/modint>
using namespace std;
using u64 = uint64_t;
using Modint = atcoder::modint998244353;


u64 binary_gcd(u64 x, u64 y) {
    if (!x or !y) return x + y;
    int n = __builtin_ctzll(x), m = __builtin_ctzll(y);
    x >>= n; y >>= m;
    while (x != y) {
        if (x > y) {
            x = (x - y) >> __builtin_ctzll(x - y);
        } else {
            y = (y - x) >> __builtin_ctzll(y - x);
        }
    }
    return x << (n > m ? m : n);
}
int main() {
    u64 N;
    cin >> N;
    
    auto divisor = [&]{
        vector<u64> divisor, divisor2;
        for(u64 d = 1; d * d <= N; d++) if(N % d == 0) {
            divisor.push_back(d);
            divisor2.push_back(N / d);
        }
        if(divisor.back() == divisor2.back()) divisor2.pop_back();
        divisor.insert(divisor.end(), divisor2.rbegin(), divisor2.rend());
        return divisor;
    }();
    
    auto press = [&](u64 x) {
        return lower_bound(divisor.begin(), divisor.end(), x) - divisor.begin();
    };
    
    vector<Modint> dp(divisor.size());
    dp[0] = 1;
    
    for(int i = 0; i < divisor.size(); i++) {
        const u64 A = divisor[i];
        for(u64 B : divisor) {
            const u64 GCD = binary_gcd(A, B);
            if(GCD == B) continue;
            const u64 LCM = A / GCD * B;
            dp[press(LCM)] += dp[i];
        }
    }
    
    cout << dp.back().val() << endl;
}
0