#include using std::cin; using std::cout; using std::endl; using std::vector; const int mod=1e9+7; int main(){ int N; cin>>N; int M=sqrt(N); vector dp(N); vector> sum(M+1,vector(M+1)); dp[0]=1; vector A(N); for(int i=0;i>A[i]; if(i>0&&A[i-1]>1){ dp[i]+=dp[i-1]; dp[i]%=mod; } for(int j=1;j<=M;j++){ 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(int j=i+A[i];j