最大公约数--------------------------------数论(线性欧拉筛)

给定整数N，求1<=x,y<=N且GCD(x,y)为素数的数对(x,y)有多少对。

GCD(x,y)即求x，y的最大公约数。

输入格式
输入一个整数N

输出格式
输出一个整数，表示满足条件的数对数量。

数据范围
1≤N≤107
输入样例：
4
输出样例：
4

解析:

gcd(x,y)=p p是质数
所以我们可以求出1~N所有质数p，枚举p
因为 gcd(x,y)=p;
所以 gcd(x/p,y/p)=1
所以问题就转换成
求 gcd(x',y')=1的个数
1<=x',y'<=N/p;

所以第一步线性欧拉筛处理出1~N所有的质数和1~N内互质的个数
然后枚举质数p,
对答案的贡献就是 (phi[1]+phi[2]+.....+phi[N/p] )*2+1
前面一部分可以前缀和处理。为什么+1,因为还存在(1,1)这种情况


本题phi[1]=1不用写 ，因为1号点是没有的

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e7+100;
int prime[N],n,cnt;
int phi[N]; 
bool st[N];
ll s[N];
void init(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
		{
			phi[i]=i-1;
			prime[++cnt]=i;
		 } 
		 for(int j=1;j<=cnt&&prime[j]<=n/i;j++)
		 {
		 	st[i*prime[j]]=true;
		 	if(i%prime[j]==0) 
		 	{
		 		phi[i*prime[j]]=phi[i]*prime[j];
		 		break;
			 }
			phi[i*prime[j]]=phi[i]*(prime[j]-1);
		 }
	} 
	for(int i=1;i<=N;i++) s[i]=s[i-1]+phi[i];
}
int main()
{
	cin>>n;
	ll res=1;
	init(n);
	for(int i=1;i<=cnt;i++)
	{
		int p=prime[i];
		res+=s[n/p]*2+1;
		
	}
	cout<<res<<endl;
}