結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0