[hdu6134]Battlestation Operational

#include<bits/stdc++.h>
using namespace std;
#define mod 1000000007
#define N 1000005
long long n,ans,mu[N],f[N],t[N],vis[N],p[N];
void xxs(int n){
    mu[1]=f[1]=1;
    for(int i=2;i<=n;i++){
        if (!vis[i]){
            p[++p[0]]=i;
            mu[i]=-1;
            t[i]=1;
            f[i]=2;
        }
        for(int j=1;(j<=p[0])&&(i*p[j]<=n);j++){
            vis[i*p[j]]=1;
            if (i%p[j]==0){
                t[i*p[j]]=t[i]+1;
                f[i*p[j]]=f[i]/(t[i]+1)*(t[i]+2);
                break;
            }
            t[i*p[j]]=1;
            f[i*p[j]]=f[i]*2;
            mu[i*p[j]]=-mu[i];
        }
    }
    t[1]=1;
    for(int i=2;i<=n;i++){
        t[i]=(f[i-1]+t[i-1]+1)%mod;
        mu[i]+=mu[i-1];
    }
    for(int i=2;i<=n;i++)f[i]=(t[i]+f[i-1])%mod;
}
int main(){
    xxs(N-5);
    while (scanf("%lld",&n)!=EOF){
        ans=0;
        for(int i=1,j;i<=n;i=j+1){
            j=n/(n/i);
            ans=(ans+(mu[j]-mu[i-1]+mod)%mod*f[n/i])%mod;
        }
        printf("%lld\n",ans);
    }
}