from numpy import *
N=int(input())
A=int64(input().split()+[1])
s=ones((20,N+1),int64)
s[0]=A
for n in range(19):
  d=1<<n
  s[n+1,:-d]=gcd(s[n,:-d],s[n,d:])
L=arange(N)
R=copy(L)
G=zeros(N,int64)
for n in range(19, -1, -1):
  H=gcd(G,s[n][R])
  g=H!=1
  R[g]+=1<<n
  G[g]=H[g]
print(N*(N+1)//2-sum(R-L))