結果
問題 | No.2798 Multiple Chain |
ユーザー | 👑 p-adic |
提出日時 | 2024-07-07 20:21:45 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 781 bytes |
コンパイル時間 | 152 ms |
コンパイル使用メモリ | 82,400 KB |
実行使用メモリ | 77,976 KB |
最終ジャッジ日時 | 2024-07-07 20:21:54 |
合計ジャッジ時間 | 6,001 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | RE | - |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | RE | - |
testcase_17 | RE | - |
testcase_18 | RE | - |
testcase_19 | RE | - |
testcase_20 | RE | - |
testcase_21 | RE | - |
testcase_22 | RE | - |
testcase_23 | RE | - |
testcase_24 | RE | - |
testcase_25 | RE | - |
testcase_26 | RE | - |
testcase_27 | RE | - |
testcase_28 | RE | - |
testcase_29 | RE | - |
testcase_30 | TLE | - |
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 | -- | - |
testcase_48 | -- | - |
testcase_49 | -- | - |
testcase_50 | -- | - |
testcase_51 | -- | - |
testcase_52 | -- | - |
testcase_53 | -- | - |
ソースコード
R=range n=int(input()) P,t=set(),n for i in R(2,min(n+1,31623)): if t%i<1:P.add(i) while t%i<1:t//=i D=[t]*(t>1) def G(a,b): while a:a,b=b%a,a return b def r(t): for j in R(len(D)): while D[j]%t<1:D[j]//=t P.add(t) while D: t=D.pop() c=int(t**0.25) if t<10**9:r(t) else: while c: x,d=c,1 y=(x*x+1)%t while (d<2)*c:x,y,d,c=(x*x+1)%t,(y**4+2*y*y+2)%t,G(abs(x-y),t),c-1 if d<t and c:D+=[t//d,d];break else:r(t) P=list(P) L=1 C=0 for i in R(len(P)): p,c=P[i],0 while n%p<1:n//=p;c+=1 P[i]=[p,c] L+=c C=max(c,C) assert(False) D=[[s<1]*(s+1)for s in R(C+1)] for l in R(L): t=[[0]for s in R(C+1)] for s in R(C+1): for d in R(s+1):t[s]+=[t[s][-1]+D[s][d]] D=[[t[s-d][min(s-d,d)+1]for d in R(s+1)]for s in R(C+1)] a=1 for[p,c]in P:a*=t[c][c+1] print(a)