bzoj2818 Gcd

题目链接：bzoj2818
题目大意：
给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对.

题解：
欧拉函数或莫比乌斯反演
设d∈[1,n]是一个质数。
原式即

$$\sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{gcd(i,j)=d}1$$
套路一下： $$\sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{d}\sum_{gcd(\frac{i}{d},\frac{j}{d})=1}1$$
额第三个求和是在枚举范围内的质数。

欧拉函数的做法就是化成：

$$\sum_{d}\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(\frac{i}{d},\frac{j}{d})=1]$$
然后看式子的后半部分，就是在求[1,n]中的与 $\frac{i}{d}$互质的数的个数，就是 $\phi(i)$。所以可以预处理欧拉函数及其前缀和，枚举质数 $d$然后直接累加就好了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
#define maxn 10001000
#define N 701000

bool ispri[maxn];LL sum[maxn];
int n,cnt,phi[maxn],pri[N];
void Eular(int lim)
{
    cnt=0;phi[1]=1;
    for (int i=2;i<=lim;i++)
    {
        if (!ispri[i]) {pri[++cnt]=i;phi[i]=i-1;}
        for (int j=1;j<=cnt && pri[j]*i<=lim;j++)
        {
            ispri[i*pri[j]]=true;
            if (i%pri[j]==0)
            {
                phi[i*pri[j]]=pri[j]*phi[i];
                // mu[i*pri[j]]=0;
                break;
            }
            phi[i*pri[j]]=(pri[j]-1)*phi[i];
            // mu[i*pri[j]]=-mu[i];
        }
    }sum[0]=0;
    for (int i=1;i<=lim;i++) sum[i]=sum[i-1]+phi[i];
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int n,i;LL ans=0;
    scanf("%d",&n);
    Eular(n);
    for (i=1;i<=cnt;i++)
    {
        int lim=n/pri[i];
        ans+=sum[lim]*2-1;
    }
    printf("%lld\n",ans);
    return 0;
}

莫比乌斯反演的做法我现在还不是很懂怎么化出来..先挖个坑
来填坑
莫比乌斯反演的做法则是化成：

$$\sum_{d}\sum_{i=1}^{n}\sum_{j=1}^{n}\epsilon((\frac{i}{d},\frac{j}{d}))$$
把 $\frac{i}{d}$、 $\frac{j}{d}$什么的换成 $i$、$j$。
反演一下： $$\sum_{d}\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{n}{d}}\sum_{t|i,t|j}\mu(t)$$
即 $$\sum_{d}\sum_{t=1}^{n}\sum_{i=1,t|i}^{\frac{n}{d}}\sum_{j=1,j|t}^{\frac{n}{d}}\mu(t)$$
实际上后半部分就是在求有多少对 $i$,$j$是 $t$的倍数，然后再乘上$\mu(t)$。
所以式子就是： $$\sum_{d}\sum_{t=1}^{n}\frac{\frac{n}{d}}{t}\times \frac{\frac{n}{d}}{t}\times \mu(t)$$
预处理 $mobius$函数及其前缀和，枚举质数 $d$，然后后面的$\sqrt{n}$分块求就好了。
懒得再打了。。
看这位 dalao的吧，谢谢dalao的解疑。