n=int(input())
a=list(map(int,input().split()))
dp=[0]*(n+1)
mod=10**9+7
dp[1]=1
d={}
for i in range(1,n+1):
  # マスiに到達するケースを加算
  for j in range(1,n+1):
    if (j,i%j) in d:
      dp[i]+=d[(j,i%j)]
      dp[i]%=mod
  if i==n:break
  x=a[i-1]
  if x>1:dp[i+1]+=dp[i]
  if (x,i%x) in d:
    d[(x,i%x)]+=dp[i]
    d[(x,i%x)]%=mod
  else:
    d[(x,i%x)]=dp[i]
print(dp[-1])