結果
問題 | No.2249 GCDistance |
ユーザー | 👑 Kazun |
提出日時 | 2023-03-17 22:25:14 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,567 ms / 5,000 ms |
コード長 | 967 bytes |
コンパイル時間 | 323 ms |
コンパイル使用メモリ | 82,440 KB |
実行使用メモリ | 243,648 KB |
最終ジャッジ日時 | 2024-09-18 11:34:09 |
合計ジャッジ時間 | 20,258 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 10 |
ソースコード
#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())))