結果
| 問題 |
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())))
Kazun