結果
問題 | No.1255 ハイレーツ・オブ・ボリビアン |
ユーザー |
👑 ![]() |
提出日時 | 2020-08-10 04:06:49 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,002 bytes |
コンパイル時間 | 168 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 52,480 KB |
最終ジャッジ日時 | 2024-06-28 05:39:37 |
合計ジャッジ時間 | 1,453 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | WA * 1 |
other | WA * 15 |
ソースコード
#素因数分解 def Prime_Factorization(N): if N<0: R=[[-1,1]] else: R=[] N=abs(N) k=2 while k*k<=N: if N%k==0: C=0 while N%k==0: C+=1 N//=k R.append([k,C]) k+=1 if N!=1: R.append([N,1]) if not R: R.append([N,1]) return R #Euler's Totient関数 def Euler_Totient(N): N=abs(N) if N==1: return 1 H=Prime_Factorization(N) R=1 for (p,e) in H: R*=p**(e-1)*(p-1) return R #約数全部 def Divisors(N): N=abs(N) L,U=[],[] k=1 while k*k <=N: if N%k== 0: L.append(k) if k!=N//k: U.append(N//k) k+=1 return L+U[::-1] #================================================ N=int(input()) M=2*N-1 D=Euler_Totient(M) L=Divisors(D) T=2 x=1 for a in L: U=pow(2,a-x,M) T*=U T%=M x=a if (T-1)%M==0: print(a) break