結果

問題 No.2318 Phys Bone Maker
ユーザー aaaaaaaaaa2230
提出日時 2023-05-27 10:20:54
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 921 ms / 3,000 ms
コード長 974 bytes
コンパイル時間 221 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 78,208 KB
最終ジャッジ日時 2024-12-25 20:18:09
合計ジャッジ時間 11,516 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 45
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
n = int(input())
mod = 998244353

divs = set()
for i in range(1,int(n**0.5)+1):
    if n%i:
        continue
    divs.add(i)
    divs.add(n//i)

divs = sorted(divs)
le = len(divs)
primes = []


for i in range(2,int(n**0.5)+1):
    if n%i:
        continue
    primes.append(i)
    while n % i == 0:
        n //= i

if n != 1:
    primes.append(n)

div_count = [[0]*len(primes) for i in range(le)]

for i,div in enumerate(divs):
    now = div
    for j,prime in enumerate(primes):
        while now%prime == 0:
            now //= prime
            div_count[i][j] += 1

dp = [0]*le
dp[0] = 1
lp = len(primes)

for i in range(le-1):
    divi = div_count[i]
    for j in range(i+1,le):
        if divs[j]%divs[i]:
            continue
        divj = div_count[j]
        pat = 1
        for x in range(lp):
            if divi[x] == divj[x]:
                pat *= divi[x]+1
                pat %= mod
        dp[j] += dp[i]*pat
        dp[j] %= mod
print(dp[-1])
0