結果

問題 No.3505 Sum of Prod of Root
コンテスト
ユーザー kitatai
提出日時 2026-04-18 17:45:58
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
TLE  
実行時間 -
コード長 1,159 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 199 ms
コンパイル使用メモリ 84,908 KB
実行使用メモリ 653,824 KB
最終ジャッジ日時 2026-04-18 17:46:28
合計ジャッジ時間 6,682 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 1
other AC * 2 TLE * 1 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

N = int(input())
M = 998244353

next_u = [1] * 60
next_ui = [1] * 60
for k in range(2, 60):
    next_ui[k] += 1
    next_u[k] = next_ui[k] ** k
u_time = [1]
u_change = []
v = 1
while v <= N:
    min_u = next_u[2]
    min_uk = 2
    for k in range(2, 60):
        if min_u > next_u[k]:
            min_u = next_u[k]
            min_uk = k
    ks = []
    for k in range(2, 60):
        if next_u[k] == min_u:
            ks.append(k)
    for k in ks:
        next_ui[k] += 1
        next_u[k] = next_ui[k] ** k
    if u_time[-1] != min_u:
        u_time.append(min_u)
        u_change.append(ks)
    v = min_u
ans = 0
f = False
max_c = int(N ** (1 / 2)) + 10
inv = [0] * (max_c + 1)
inv[1] = 1
for i in range(2, max_c + 1):
    inv[i] = M - (M // i) * inv[M % i] % M
cu = [1] * 60
p = 1
for i in range(len(u_time) - 1):
    u_i = u_time[i]
    u_i1 = u_time[i + 1]
    if u_i1 > N:
        u_i1 = N + 1
        f = True
    ans = (ans + (p * ((u_i1 - u_i) * (u_i + u_i1 - 1) // 2))) % M
    if f:
        break
    for k in u_change[i]:
        old = cu[k]
        cu[k] += 1
        new = cu[k]
        p = p * inv[old] % M
        p = p * new % M
print(ans)
0