結果
| 問題 | No.2075 GCD Subsequence |
| コンテスト | |
| ユーザー |
sasa8uyauya
|
| 提出日時 | 2024-09-01 18:12:44 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 3,427 ms / 4,000 ms |
| コード長 | 778 bytes |
| 記録 | |
| コンパイル時間 | 421 ms |
| コンパイル使用メモリ | 82,500 KB |
| 実行使用メモリ | 179,908 KB |
| 最終ジャッジ日時 | 2024-09-01 18:13:44 |
| 合計ジャッジ時間 | 59,152 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 28 |
ソースコード
n=int(input())
a=list(map(int,input().split()))
L=1000000
p=[[] for i in range(L+1)]
p[1]=[1]
for i in set(a):
if i==1:
continue
v=i
j=1
while j*j<=v:
if v%j==0:
p[i]+=[j]
if (v//j)>j:
p[i]+=[v//j]
j+=1
P=[1]*L
for i in range(2,L):
if P[i]:
for j in range(i+i,L,i):
P[j]=0
P=[i for i in range(2,L) if P[i]]
m=[0]*(L+1)
for i in range(L+1):
if i==1:
continue
m[i]=1
v=i
for j in P:
if j*j>v:
break
if v%j==0:
m[i]*=-1
if v%(j*j)==0:
m[i]=0
while v%j==0:
v//=j
j+=1
if v>1:
m[i]*=-1
M=998244353
g=[0]*(L+1)
q=[0]*n
for i in range(n):
f=1
for v in p[a[i]]:
f+=(-m[v])*g[v]
q[i]=f%M
for v in p[a[i]]:
g[v]+=q[i]
g[v]%=M
print(sum(q)%M)
sasa8uyauya