def main1(n,a):
  dp=[0]*n
  mod=10**9+7
  dp[0]=1
  d={}
  base=0
  c=int(n**0.5)+1
  for i in range(n):
    # マスiに到達するケースを加算
    for j in range(2,min(i,c)):
      if (j,i%j) in d:
        dp[i]+=d[(j,i%j)]
        dp[i]%=mod
    dp[i]+=base
    dp[i]%=mod
    if i==n-1:break
    k=a[i]
    if k==1:
      base+=dp[i]
      base%=mod
    elif k<c:
      dp[i+1]+=dp[i]
      dp[i+1]%=mod
      if (k,i%k) in d:
        d[(k,i%k)]+=dp[i]
        d[(k,i%k)]%=mod
      else:
        d[(k,i%k)]=dp[i]
    else:
      dp[i+1]+=dp[i]
      dp[i+1]%=mod
      ii=i+k
      while ii<n:
        dp[ii]+=dp[i]
        dp[ii]%=mod
        ii+=k
  return dp[-1]

if __name__=='__main__':
  n=int(input())
  a=list(map(int,input().split()))
  ret1=main1(n,a)
  print(ret1)