結果

問題 No.1396 Giri
コンテスト
ユーザー ああ
提出日時 2026-05-17 15:04:26
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 894 ms / 2,000 ms
コード長 723 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 408 ms
コンパイル使用メモリ 85,504 KB
実行使用メモリ 194,816 KB
最終ジャッジ日時 2026-05-17 15:04:35
合計ジャッジ時間 8,059 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge3_0
純コード判定待ち
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

n=int(input())
x=[0]*(n+1)
y=[1]*(n+1)
z=[0]*(n+1);mod=998244353
for i in range(2,n+1):
    if x[i]==0:
        for j in range(i,n+1,i):
        	if x[j]:
        		continue
        	x[j]=i
    c=x[i]
    for j in range(2,int(i**0.5)+1):
        if i%j==0:
            c=j;break
    j=i
    while j%c==0:
        y[i]*=c;z[i]+=1;j//=c
s=[[0,0] for i in range(n+1)]
for i in range(2,n+1):
    c=i
    s[x[c]][1]=max(s[x[c]][1],z[c])
    if s[x[c]][1]>s[x[c]][0]:
        s[x[c]][0],s[x[c]][1]=s[x[c]][1],s[x[c]][0]
ans=1;p=1
k=[1]*(n+1)
for i in range(2,n+1):
    ans*=pow(i,s[i][0],mod);ans%=mod
    f=k[i//y[i]]
    if z[i]>s[x[i]][1]:
    	f*=pow(x[i],s[x[i]][0]-s[x[i]][1])
    p=max(p,f)
print(ans*pow(p,mod-2,mod)%mod)
0