結果
| 問題 |
No.1262 グラフを作ろう!
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-03-12 01:59:29 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 570 ms / 3,000 ms |
| コード長 | 910 bytes |
| コンパイル時間 | 1,468 ms |
| コンパイル使用メモリ | 81,920 KB |
| 実行使用メモリ | 269,480 KB |
| 最終ジャッジ日時 | 2024-09-16 06:22:38 |
| 合計ジャッジ時間 | 50,069 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge6 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 96 |
ソースコード
def ind(b,n):
res=0
while n%b==0:
res+=1
n//=b
return res
def eratosthenes(n):
spf = [0 for i in range(n+1)]
for i in range(2,n+1):
if not spf[i]:
for j in range(1,n//i+1):
spf[i*j] = i
return spf
#spf = eratosthenes(10**6)
def euler_phi(n):
phi = [1 for i in range(n+1)]
for i in range(2,n+1):
if phi[i] == 1:
for j in range(1,n//i+1):
t = 1
tmp = j
while tmp%i==0:
t *= i
tmp //= i
phi[i*j] *= t*(i-1)
return phi
phi = euler_phi(10**6)
N,M = map(int,input().split())
A = list(map(int,input().split()))
data = [0 for i in range(10**6+1)]
for a in A:
data[a] += 1
res = - sum(A)
for i in range(1,10**6+1):
for j in range(1,10**6//i+1):
res += (i)*phi[j]*data[i*j]
print(res)