結果

問題 No.2576 LCM Pattern
コンテスト
ユーザー LyricalMaestro
提出日時 2026-09-11 00:46:24
言語 PyPy3
(7.3.23 + ACL)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 38 ms / 2,000 ms
+ 509µs
コード長 1,059 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 90 ms
コンパイル使用メモリ 80,768 KB
実行使用メモリ 64,768 KB
最終ジャッジ日時 2026-09-11 00:46:27
合計ジャッジ時間 2,494 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

## https://yukicoder.me/problems/no/2576

import math

MOD = 998244353

def main():
    N, M = map(int, input().split())

    # 素因数分解
    m = int(math.sqrt(M))

    primes = {}
    for p in range(2, m + 1):
        if M % p == 0:
            primes[p] = 0
            while M % p == 0:
                primes[p] += 1
                M //= p
    if M > 1:
        primes[M] = 1

    prime_array = [(p, e) for p, e in primes.items()]
    p_len = len(prime_array)

    dp = [0] * (2 ** p_len)
    x = 1
    for _, e in prime_array:
        x *= e
    dp[0] = pow(x, N, MOD)

    for bit in range(1, 2 ** p_len):
        x = 1
        for k in range(p_len):
            _, e = prime_array[k]
            if bit & (1 << k) > 0:
                x *= (e + 1)
            else:
                x *= e
        base = pow(x, N, MOD)
        for pre_bit in range(bit):
            if bit & pre_bit == pre_bit:
                base -= dp[pre_bit]
                base %= MOD
        dp[bit] = base
    print(dp[-1])








if __name__ == "__main__":
    main()
0