Battlestation Operational（hdu6134 莫比乌斯反演）

题目：点击打开链接

题意：f(n) = sigma[1<=i<=n]sigma[1<=j<=i]ceil[i/j] (gcd(i,j)==1)

sigma代表求和符号，celi代表除完向上取整；

解题：个人觉得这个代码还挺容易看懂的

点击打开链接

#include <bits/stdc++.h>using namespace std;typedef long long ll;const ll mod=1e9+7;const ll maxn=1e6+5;int mu[maxn],prime[maxn],not_prime[maxn];ll f[maxn],g[maxn],h[maxn];void mobius(){    ll tot;    mu[1] = 1;    tot = 0;    for(int i = 2; i < maxn; i++){        if(!not_prime[i]){            prime[++tot] = i;            mu[i] = -1;        }        for(int j = 1; prime[j]*i<maxn; j++){            not_prime[prime[j]*i] = 1;            if(i%prime[j]==0){                mu[prime[j]*i] = 0;                break;            }            mu[prime[j]*i] = -mu[i];        }    }}//莫比乌斯函数的模板void solve(){    memset(g,0,sizeof(g));    memset(h,0,sizeof(h));    memset(f,0,sizeof(f));    for(int i=1;i<=1e6;i++)    {        g[i]++;        for(int j=i+1;j<=1e6;j+=i)            g[j]++;    }    for(int i=2;i<=1e6;i++)    g[i]=(g[i]+g[i-1])%mod;    for(int i=1;i<=1e6;i++)    {        for(int j=i;j<=1e6;j+=i)        {            h[j]=(h[j]+mu[i]*g[j/i]%mod+mod)%mod;        }    }   // cout<<h[1]<<endl;    for(int i=1;i<=1e6;i++)    {        f[i]=(f[i-1]+h[i])%mod;    }}int main(){    int n;    mobius();    solve();    while(cin>>n)    {        cout<<f[n]<<endl;    }    return 0;}

                                                                                  

1],