P-2568（水题2）

第二道大佬们口中的数论水题

题目大意：给一个数n，求有多少对（x，y) (x，y均小于等于m），gcd（x，y）等于一个素数。注意：（x，y）和（y，x）不一样。

这个题和我们的水题1差不多的道理，就是方法不一样（但我还是不会）

题解：假设我们先找最大公约数为2的（x，y）对，那么我们正常的思路就是枚举每一个y%2==0的y然后将这些和y/2互素的数（欧拉数）加起来，这样的话y就是取2,4,6,8…那么y/2就是1,2,3,4…而我们的任务是吧所有的素数都求一遍，所以我们不如把欧拉数组的前缀和算出来，这样就不用一个一个算了。

代码：

//#pragma GCC optimize(3,"Ofast","inline")
//#include<unordered_map>
//#include<unordered_set>
#include<cstdio>
#include<iostream>
#include<cmath>
#include<functional>
#include<cstring>
#include<string>
#include<cstdlib>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#include<stack>
#include<vector>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e7+10;
const ll mod=1e9+7;
int prime[maxn];
int isprime[maxn];
int phi[maxn];
int num;
ll sum[maxn];
int n;
void getprime()
{
    memset(prime,0,sizeof(prime));
    memset(isprime,1,sizeof(isprime));
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(isprime[i])
        {
            phi[i]=i-1;
            prime[num++]=i;
        }
        for(int j=0;j<num&&i*prime[j]<=n;j++)
        {
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}
            else phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
}
int main()
{
    cin>>n;
    getprime();//欧拉筛。。。。
    ll ans=0;
    for(int i=1;i<=n;i++)sum[i]=sum[i-1]+phi[i];
    for(int i=0;i<num;i++)ans+=sum[n/prime[i]]*2-1;//为什么乘二？因为xy和yx不同，为什么减一？因为phi[1]=1,两个（1，1）要减去一个，假设当前素数为2，也就是不能有两个（2，2）
    cout<<ans<<endl;
    return 0;
}