N = int(input())
A = list(map(int,input().split()))

assert 1<=N<=10**5
assert len(A)==N
assert all(1<=A[i]<=N for i in range(N))

A = [0] + A
mod = 10**9 + 7

dp = [0 for i in range(N+1)]
data = [[0 for i in range(320)] for j in range(N)]
dp[-1] = 1
for j in range(1,min(N,320)):
    data[N-j][j] = 1

for i in range(N-1,0,-1):
    if A[i]**2 > N:
        for j in range(A[i],N+1-i,A[i]):
            dp[i] += dp[i+j]
            dp[i] %= mod
    else:
        dp[i] = data[i][A[i]]

    if A[i]!=1:
        dp[i] += dp[i+1]
        dp[i] %= mod

    for j in range(1,min(320,i)):
        data[i - j][j] = data[i][j] + dp[i]
        data[i - j][j] %= mod

print(dp[1])