hdu6134[莫比乌斯反演] 2017多校8

/*hdu6134[莫比乌斯反演] 2017多校8*/
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL MOD = ;
const int maxn =  + ;
LL f[maxn], g[maxn], d[maxn];
bool vis[maxn];
int prime[maxn], primes, mu[maxn];
void init_mu()
{
    memset(vis, , sizeof(vis));
    mu[] = ;
    primes = ;
    for (int i = ; i < maxn; i++)
    {
        if (!vis[i]) {
            prime[primes++] = i;
            mu[i] = -;
        }
        for (int j = ; j < primes && i * prime[j] < maxn; j++)
        {
            vis[i * prime[j]] = ;
            if (i % prime[j]) mu[i * prime[j]] = -mu[i];
            else { mu[i * prime[j]] = ; break;}
        }
    }
}
void init() {
    init_mu();
    for (int i = ; i < maxn; i++) {
        for (int j = i; j < maxn; j += i) {
            d[j]++;
            if (d[j] > MOD) d[j] -= MOD;
        }
    }
    g[] = ;
    for (int i = ; i < maxn; i++) {
        g[i] = (g[i - ] + d[i - ] + ) % MOD;
    }
    for (int i = ; i < maxn; i++) {
        for (int j = i; j < maxn; j += i) {
            f[j] = (f[j] + mu[i] * g[j / i] % MOD) % MOD;
        }
    }
    for (int i = ; i < maxn; i++) {
        f[i] = (f[i - ] + f[i]) % MOD;
    }
}
int main() {
    init();
    int n;
    while (~scanf("%d", &n)) {
        printf("%lld\n", f[n]);
    }
    return ;
}