GCD2[caioj1281]\[BZOJ2818]莫反

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

caioj
BZOJ

思路

这是一道不用莫函的莫反(伪)题啊。

我们观察一下,

$$1\le x,y\le n,\gcd (x,y)=p,p是质数$$

可以转化为:

$$1\le x,y\le n/p,\gcd (x,y)=1，也就是x,y互质$$

那么我们先考虑$y\ge x$的情况

当$y=x$时，只有在$y=1,x=1$时，$\gcd (x,y)=1$成立。

当$y>x$时，恰好$\forall y\in [2,n/p],\forall x\in [1,y)$，与y互质的$x$，共有$\phi (y)$个。

通俗地说，就是对于每一个$y$,都会有$\phi (y)$个数与它互质。

对于每一个质数$p$，当$y\ge x,y\in [1,n/p]$时，$ans=({\sum }_{i=2}^{n/p}\phi (i))+1={\sum }_{i=1}^{n/p}\phi (i)$

那么$x,y$没有大小关系时，$ans=({\sum }_{i=1}^{n/p}\phi (i)\ast 2)-1$，首先$\ast 2$是因为，$(2,3),(3,2)$是两对不同的二元组，$-1$是因为$y=x=1$，即$(1,1)$的情况被选了两遍，要减去$1$。

那么对于$\forall 质数p\in [1,n]$，都有$({\sum }_{i=1}^{n/p}\phi (i)\ast 2)-1$。

则最终答案为:

$$\sum _{p是质数}^n(\sum _{i=1}^{n/p}(\phi (i)\ast 2)-1)$$

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ll long long
using namespace std;
const int N=1e7+10;
const int M=1e6+10;
int prime[M],m,n;ll phi[N];bool v[N];
inline void g_p()
{
    m=0;phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!v[i])prime[++m]=i,phi[i]=i-1;
        for(int j=1;j<=m&&i*prime[j]<=n;j++)
        {
            v[i*prime[j]]=1;
            if(i%prime[j]==0){phi[i*prime[j]]=prime[j]*phi[i];break;}
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
    for(int i=2;i<=n/2;i++)phi[i]+=phi[i-1];
}
int main()
{
    scanf("%d",&n);
    g_p();
    ll ans=0;
    for(int i=1;i<=m;i++)
        ans+=phi[n/prime[i]]*2-1;//1,1
    printf("%lld\n",ans);
    return 0;
}