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+1)
m=[1]*(L+1)
for i in range(2,L):
  if P[i]:
    m[i]*=-1
    for j in range(i+i,L,i):
      P[j]=0
      if j%(i*i)==0:
        m[j]=0
      else:
        m[j]*=-1

M=998244353
g=[0]*(L+1)
q=[0]*n
for i in range(n):
  f=1+g[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)