51Nod 1237 最大公约数之和 V3 (杜教筛)

题目链接：哆啦A梦传送门

题意：

求：$\begin{array}{r}\sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)\end{array}$

题解：

$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)\\ =& \sum _{d=1}^{n}d\sum _{i=1}^n\sum _{j=1}^{n}gcd(i,j)=d\\ =& \sum _{d=1}^{n}d\sum _{i=1}^{\frac{n}{d}}\sum _{j=1}^{\frac{n}{d}}gcd(i,j)=1\\ =& \sum _{d=1}^{n}d((2\sum _{i=1}^{\frac{n}{d}}\sum _{j=1}^i[(i,j)=1])-1).......解释1\\ =& \sum _{d=1}^{n}d((2\sum _{i=1}^{\frac{n}{d}}\varphi (i)])-1)\end{array}$

我们再设：
$S(n)={\sum }_{i=1}^n\varphi (i)$

解释1见：
故只需杜教筛下欧拉函数前n项和就可以了。
杜教筛求欧拉函数前n项和

那么最后题目就变为：
$$S(n)=\sum _{i=1}^n-\sum _{d=2}^{n}S((\lfloor \frac{n}{d}\rfloor )$$
$$ans=\sum _{d=1}^{d}d\ast (2\ast S(\frac{n}{d})-1)$$

分块下解决就好了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<tr1/unordered_map>

using namespace std;

typedef long long LL;
const int N=6e6;
const LL mod=1e9+7;

bool vis[N+10];
int pri[N+10],phi[N+10],tot;
tr1::unordered_map<LL,LL> w;
LL sum[N+10];
LL inv2=500000004;
void init()
{
    phi[1]=1;

    for(int i=2;i<=N;i++)
    {
        if(!vis[i]){
            pri[++tot]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=tot&&i*pri[j]<=N;j++)
        {
            int x=pri[j];
            vis[i*x]=true;
            if(i%x==0){
                phi[i*x]=phi[i]*x;break;
            }
            phi[i*x]=phi[i]*(x-1);
        }
    }

    ///预处理前N项欧拉函数和
    for(int i=1;i<=N;i++)
        sum[i]=(sum[i-1]+phi[i])%mod;

}

LL s1(LL x) ///前x项和
{
    x=x%mod;
    return x*(x+1)%mod*inv2%mod;
}
///杜教筛求欧拉函数前n项和
LL djs(LL x)
{
    if(x<=N) return sum[x];
    if(w[x]) return w[x];

    LL ans=s1(x);

    for(LL l=2,r;l<=x;l=r+1)
    {
        r=x/(x/l);
        ans=(ans-(r-l+1)%mod*djs(x/l)%mod+mod)%mod;
    }
    return w[x]=ans;
}

int main()
{

    LL n;
    init();
    scanf("%lld",&n);

    LL ans=0;

    for(LL l=1,r;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans=(ans+(s1(r)-s1(l-1)+mod)%mod*(2*djs(n/l)%mod-1+mod)%mod)%mod;
    }
    printf("%lld\n",ans%mod);
    return 0;



}