51nod(1237最大公约数之和 V3 反演+杜教筛)

题目

这里需要反演一波。。。

最后一个式子就可以整除分块+杜教筛。。。

#include<cstdio>
#include<map>
using namespace std;
typedef long long ll;
const int N=5e6+5,mod=1e9+7,inv=500000004;
int is[N+5],prime[N+5];
int phi[N+5];
void init()
{
    phi[1]=1;int tot=0;
    for(int i=2;i<=N;++i)
    {
        if(!is[i]) prime[++tot]=i,phi[i]=i-1;
        for(int j=1;j<=tot&&(ll)i*prime[j]<=N;++j)
        {
            is[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j]%mod;
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1)%mod;
        }
    }
    for(int i=2;i<=N;++i) (phi[i]+=phi[i-1])%=mod;
}
map<int,int>ans_phi;
inline int f_phi(ll n){
    if(n<=N) return phi[n];
    if(ans_phi[n]) return ans_phi[n];
    ll res=(n%mod)*((n+1)%mod)%mod*inv%mod;
    for(ll l=2,r;l<=n;l=r+1)
        r=n/(n/l),res-=(r-l+1)*f_phi(n/l)%mod;
    return ans_phi[n]=(res%mod+mod)%mod;
}
int main(){
    init();
    ll n,ans=0;scanf("%lld",&n);
    for(ll l=1,r;l<=n;l=r+1)
        r=n/(n/l),(ans+=(n/l%mod)*((n/l)%mod)%mod*((f_phi(r)-f_phi(l-1))%mod))%=mod;
    printf("%lld\n",(ans+mod)%mod);
}