結果
問題 | No.2249 GCDistance |
ユーザー | 👑 Kazun |
提出日時 | 2023-03-17 22:25:14 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,622 ms / 5,000 ms |
コード長 | 967 bytes |
コンパイル時間 | 213 ms |
コンパイル使用メモリ | 81,700 KB |
実行使用メモリ | 243,184 KB |
最終ジャッジ日時 | 2023-10-18 15:22:47 |
合計ジャッジ時間 | 20,340 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1,468 ms
231,452 KB |
testcase_01 | AC | 1,536 ms
243,180 KB |
testcase_02 | AC | 1,531 ms
243,180 KB |
testcase_03 | AC | 1,522 ms
243,180 KB |
testcase_04 | AC | 1,496 ms
235,772 KB |
testcase_05 | AC | 1,537 ms
243,184 KB |
testcase_06 | AC | 1,622 ms
243,180 KB |
testcase_07 | AC | 1,525 ms
243,180 KB |
testcase_08 | AC | 1,476 ms
231,880 KB |
testcase_09 | AC | 1,529 ms
243,184 KB |
testcase_10 | AC | 1,541 ms
241,536 KB |
ソースコード
#Euler's Totient関数 def Euler_Totient_List(N): """k=0,1,...,N に対して, 1以上k以下の整数のうち, kと互いに素な整数の個数 φ(k) を求める. N:正の整数 """ assert N>=0,"Nが非負整数ではない." phi=list(range(N+1)) for p in range(2,N+1): if phi[p]==p: for j in range(p,N+1,p): phi[j]=phi[j]//p*(p-1) return phi #================================================== def solve(): T=int(input()) N=[0]*T for t in range(T): N[t]=int(input()) N_max=max(N) P=Euler_Totient_List(N_max) Pc=[0]*(N_max+1); Pc[0]=-1 for n in range(1,N_max+1): Pc[n]=Pc[n-1]+P[n] Ans=[0]*T for t in range(T): K=N[t]*(N[t]-1)//2 Ans[t]=Pc[N[t]]+2*(K-Pc[N[t]]) return Ans #================================================== import sys input=sys.stdin.readline write=sys.stdout.write write("\n".join(map(str, solve())))