Acwing 220. 最大公约数

题目描述

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

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

输入格式
输入一个整数 N。

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

数据范围
$1\le N\le 10^7$

思路

结论：$1-N$的正整数中，互质二元组的个数为 $-1+2\ast {\sum }_{i=1}^{N}phi[i]$.

所以本题，首先筛 $1-n$ 的质数和欧拉函数，维护欧拉函数前缀和，

答案即为 ${\sum }_{k=1}^{cnt}(-1+2\ast {\sum }_{i=1}^{\lfloor \frac{n}{primes[k]}\rfloor }phi[i])$

代码

#include<iostream>

using namespace std;
typedef long long ll;

const ll N=1e7+10;
ll n,cnt;

ll primes[N],phi[N];
bool st[N];

void m_init()
{
	ll i,j;
	
	phi[1]=1;
	
	for(i=2;i<=n;++i)
	{
		if(!st[i])	primes[cnt++]=i,phi[i]=i-1;
		for(j=0;primes[j]<=n/i;++j)
		{
			st[primes[j]*i]=true;
			
			if(i%primes[j]==0)
			{
				phi[primes[j]*i]=phi[i]*primes[j];
				break;	
			}	
			phi[primes[j]*i]=phi[i]*(primes[j]-1);
		}	
	}	
	
	for(i=1;i<=n;++i)
		phi[i]+=phi[i-1];
}

int main()
{
	ll i,ans=0;
	
	scanf("%d",&n);
	m_init();
	for(i=0;i<cnt;++i)
		ans+=2*phi[n/primes[i]]-1;
	
	printf("%lld",ans);
	
	return 0;
}