求GCD(x,y) 为素数的数对 (x,y) 有多少对。
关键:将x y 拆成 a * d,b * d,(d为素数)、前缀和优化、phi[1]的另外处理.
题解转自秦大佬题解。
代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 10000010;
ll prime[N],phi[N],n,idx,sum[N],ans;
bool st[N];
void get_phi(int n)//欧拉函数之和
{
phi[1] = 1;
for(int i = 2; i <= n; i ++ )
{
if(!st[i])
{
prime[idx ++ ] = i;
phi[i] = i - 1;
}
for(int j = 0; i * prime[j] <= n; j ++ )
{
st[prime[j] * i] = true;
if(i % prime[j] == 0)
{
phi[i * prime[j]] = phi[i] * prime[j];
break;
}
phi[i * prime[j]] = phi[i] * (prime[j] - 1);
}
}
}
int main()
{
cin>>n;
get_phi(n);
for(int i = 2; i <= n; i ++ )//前缀和优化,1另外处理
sum[i] = sum[i - 1] + phi[i];
for(int i = 0; i < idx; i ++ )
{
int k = n / prime[i];
ans += sum[k] * 2 + 1;//1 x y没有差别不用乘2
}
cout<<ans;
return 0;
}