mod=10**9+7
n=int(input())
a=[0]+list(map(int,input().split()))
dp=[0]*(n+3)
dp[1]=1
B=400
sdp=[[0]*(B+3) for i in range(B+3)]
for i in range(1,n+1):
    dp[i]+=dp[i-1]
    for j in range(1,B+1):
        dp[i]+=sdp[j][i%j]
    dp[i]%=mod
    if a[i-1]==1:dp[i]-=dp[i-1]
    dp[i]%=mod
    if a[i]<=B:
        sdp[a[i]][i%a[i]]+=dp[i]
        sdp[a[i]][i%a[i]]%=mod
    else:
        for j in range(i+a[i],n+1,a[i]):dp[j]+=dp[i]
print(dp[n]%mod)