結果
問題 | No.2249 GCDistance |
ユーザー |
![]() |
提出日時 | 2023-03-17 22:40:26 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,187 bytes |
コンパイル時間 | 342 ms |
コンパイル使用メモリ | 81,820 KB |
実行使用メモリ | 232,940 KB |
最終ジャッジ日時 | 2024-09-18 11:47:32 |
合計ジャッジ時間 | 18,566 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 1 TLE * 9 |
ソースコード
import sysinput = sys.stdin.readline# エラトステネスの篩を用いた素因数分解・約数列挙MAX=10**7+10 # 使いたい最大値を指定# Sieve[i]で、iの最も小さい約数を返す。Sieve=[i for i in range(MAX)]for i in range(2,MAX):if Sieve[i]!=i:continuefor j in range(i,MAX,i):if Sieve[j]==j:Sieve[j]=i# 素因数分解def fact(x):D=dict()while x!=1:k=Sieve[x]if k in D:D[k]+=1else:D[k]=1x//=kreturn D# 約数列挙def faclist(x):LIST=[1]while x!=1:k=Sieve[x]count=0while x%k==0:count+=1x//=kLIST2=[]for l in LIST:for i in range(1,count+1):LIST2.append(l*k**i)LIST+=LIST2return LISTdef euler_totient(x):FACT=fact(x)ANS=xfor f in FACT:ANS=ANS*(f-1)//freturn ANSANS=[0]*(10**7+10)for i in range(2,10**7+1):k=euler_totient(i)ANS[i]=ANS[i-1]+i-1+(i-1-k)#print(ANS[:10])T=int(input())for tests in range(T):N=int(input())print(ANS[N])