結果

問題 No.2318 Phys Bone Maker
ユーザー FromBooskaFromBooska
提出日時 2023-05-27 10:31:26
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,562 bytes
コンパイル時間 308 ms
コンパイル使用メモリ 86,908 KB
実行使用メモリ 80,828 KB
最終ジャッジ日時 2023-08-26 19:11:22
合計ジャッジ時間 5,075 ms
ジャッジサーバーID
(参考情報)
judge15 / judge13
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 69 ms
75,800 KB
testcase_01 AC 70 ms
71,288 KB
testcase_02 TLE -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
testcase_46 -- -
testcase_47 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

# Nを素因数分解
# たとえば2**3が入っていれば、最後は2**3項が必要
# それまでは2**1, 2**2という2つのチョイスがあるので2**2で4通りある
# 逆から考えられないか
# Nを1以外の約数で割って、1までもっていく方法の数
# メモ化dfsなら実装できるが間に合わないか
# dpでできれば間に合うだろう
# 10**12でも素因数の数は忘れたけど少ない、たしか数十個ぐらい
# まずメモ化dfsでやってみるか
# いや、Nの約数にしか止まらないからそこでのdpとしたらどうだ

N = int(input())

def divisors(n):
    lower_divisors , upper_divisors = [], []
    i = 1
    while i*i <= n:
        if n % i == 0:
            lower_divisors.append(i)
            if i != n // i:
                upper_divisors.append(n//i)
        i += 1
    return lower_divisors + upper_divisors[::-1]

N_divs = divisors(N)
N_divs_set = set(N_divs)
L = len(N_divs)
dic = {}
for i in range(L):
    dic[N_divs[i]] = i


mod = 998244353

# Nの約数を大きい方から1まで並べるのはうまくいかないので小さい方からやる
# dp[i]パターン数

dp = [0]*L
dp[0] = 1 

for i in range(L-1):
    num = N_divs[i]
    #print(num)
    for d in N_divs[1:L-1]:
        #print('i', i, 'num', num, 'd', d, 'num*d', num*d, num*d in N_divs)
        if num%d != 0 and num*d in N_divs:
            dp[dic[num*d]] += dp[i]
            dp[dic[num*d]] %= mod
    dp[L-1] += dp[i]
    dp[L-1] %= mod
    #print(num, dp)

ans = dp[L-1]%mod
print(ans)

0