import math
N=int(input())
mod=(int)(1e9+7)
A=list(map(int,input().split()))
dp=[0]*N
dp[0]=1
M=int(math.sqrt(N))
Sum=[[0]*(M+1) for i in range(M+1)]
for i in range(N):
    if i>0 and A[i-1]!=1:
        dp[i]+=dp[i-1];
        dp[i]%=mod
    for j in range(M+1):
        if j>0:
            dp[i]+=Sum[j][i%j]
            dp[i]%=mod
    if A[i]<=M:
        Sum[A[i]][i%A[i]]+=dp[i]
        Sum[A[i]][i%A[i]]%=mod
    else:
        for j in range(N):
            X=j*A[i]+A[i]+i
            if X>=N:
                break
            dp[X]+=dp[i]
            dp[X]%=mod
print(dp[N-1])