結果

問題 No.2075 GCD Subsequence
ユーザー sasa8uyauyasasa8uyauya
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,027 ms
137,220 KB
testcase_01 AC 996 ms
137,180 KB
testcase_02 AC 1,023 ms
137,256 KB
testcase_03 AC 993 ms
137,132 KB
testcase_04 AC 1,025 ms
137,084 KB
testcase_05 AC 988 ms
137,264 KB
testcase_06 AC 1,023 ms
137,124 KB
testcase_07 AC 992 ms
137,160 KB
testcase_08 AC 2,523 ms
164,476 KB
testcase_09 AC 3,427 ms
179,908 KB
testcase_10 AC 2,641 ms
165,576 KB
testcase_11 AC 3,107 ms
173,388 KB
testcase_12 AC 2,892 ms
169,048 KB
testcase_13 AC 2,301 ms
161,532 KB
testcase_14 AC 3,079 ms
172,756 KB
testcase_15 AC 2,429 ms
162,924 KB
testcase_16 AC 2,551 ms
164,588 KB
testcase_17 AC 3,353 ms
179,836 KB
testcase_18 AC 1,616 ms
146,644 KB
testcase_19 AC 1,629 ms
146,864 KB
testcase_20 AC 1,605 ms
146,780 KB
testcase_21 AC 1,638 ms
146,908 KB
testcase_22 AC 1,604 ms
146,744 KB
testcase_23 AC 1,607 ms
146,572 KB
testcase_24 AC 1,635 ms
146,788 KB
testcase_25 AC 1,599 ms
146,684 KB
testcase_26 AC 1,592 ms
146,536 KB
testcase_27 AC 1,599 ms
146,648 KB
testcase_28 AC 1,536 ms
144,612 KB
testcase_29 AC 1,021 ms
137,232 KB
testcase_30 AC 996 ms
137,324 KB
権限があれば一括ダウンロードができます

ソースコード

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