結果

問題 No.1781 LCM
ユーザー 👑 NachiaNachia
提出日時 2021-12-10 00:23:04
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,647 bytes
コンパイル時間 1,451 ms
コンパイル使用メモリ 82,484 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-17 17:31:46
合計ジャッジ時間 1,932 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 AC 1 ms
6,944 KB
testcase_08 AC 2 ms
6,940 KB
testcase_09 AC 16 ms
6,944 KB
testcase_10 AC 15 ms
6,944 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 2 ms
6,940 KB
testcase_13 AC 15 ms
6,944 KB
testcase_14 AC 15 ms
6,944 KB
testcase_15 AC 7 ms
6,944 KB
testcase_16 AC 6 ms
6,940 KB
testcase_17 AC 10 ms
6,944 KB
testcase_18 AC 8 ms
6,944 KB
testcase_19 AC 2 ms
6,940 KB
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 AC 2 ms
6,940 KB
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 AC 2 ms
6,940 KB
testcase_32 AC 1 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <atcoder/modint>
using namespace std;
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;
#define rep(i,n) for(int i=0; i<(n); i++)


using m32 = atcoder::modint998244353;


int main() {
    i64 N,M; cin >> N >> M;
    if(200'000 < N) return 0;
    if(200'000 < M) return 0;

    vector<int> is_prime(M+1, 1); is_prime[0] = is_prime[1] = 0;
    vector<int> largest_prime_factor(M+1, 1);
    vector<int> largest_prime_factor_index(M+1, 0);
    for(i64 i=2; i<=M; i++) if(is_prime[i]) {
        for(i64 j=i*i; j<=M; j+=i) is_prime[j] = 0;
        for(i64 j=i; j<=M; j+=i) largest_prime_factor[j] = i;
    }
    for(i64 i=2; i<=M; i++){
        i64 d = largest_prime_factor[i];
        if(i/d%d == 0) largest_prime_factor_index[i] = largest_prime_factor_index[i/d] + 1;
        else largest_prime_factor_index[i] = 1;
    }
    
    vector<m32> dp(M+1, 1); dp[0] = 0;
    vector<m32> powN(100, 1); rep(i,100) powN[i] = m32(i).pow(N);
    for(int i=2; i<=M; i++){
        int pre = i;
        rep(t, largest_prime_factor_index[i]) pre /= largest_prime_factor[i];
        dp[i] = dp[pre] * powN[1 + largest_prime_factor_index[i]];
    }

    //for(auto a : dp) cout << a.val() << " "; cout << endl;

    for(int i=1; i<=M; i++) if(is_prime[i]) for(int j=M/i; j>=1; j--) dp[j*i] -= dp[j];
    m32 ans = 0;
    for(int i=1; i<=M; i++) ans += dp[i] * (M/i);
    cout << ans.val() << endl;
    return 0;
}


struct ios_do_not_sync{
    ios_do_not_sync(){
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
    }
} ios_do_not_sync_instance;
0