洛谷P2568 GCD 题解
题目链接:P2568 GCD
题意:
给定正整数 n n n,求 1 ≤ x , y ≤ n 1\le x,y\le n 1≤x,y≤n 且 gcd ( x , y ) \gcd(x,y) gcd(x,y)为素数的数对 ( x , y ) (x,y) (x,y) 有多少对。
1 ≤ n ≤ 1 0 7 1 \le n\le 10^7 1≤n≤107
正着去想不太好搞,那就反过来
考虑一个素数对答案的贡献
显然所有满足条件的 ( x , y ) (x,y) (x,y) 都可以写作 ( a p i , b p i ) (ap_i,bp_i) (api,bpi) 的形式,其中 gcd ( x , y ) = p i \gcd(x,y)=p_i gcd(x,y)=pi
注意到 gcd ( a , b ) = 1 \gcd(a,b)=1 gcd(a,b)=1
则对于每个素数
p
(
p
≤
n
)
p\ (p\le n)
p (p≤n) ,它对答案的贡献就是
2
∑
i
=
1
k
φ
(
i
)
−
1
2\sum\limits_{i=1}^{k}\varphi(i)-1
2i=1∑kφ(i)−1
其中,
k
=
⌊
n
p
⌋
k=\left\lfloor\dfrac{n}{p}\right\rfloor
k=⌊pn⌋,
φ
(
n
)
\varphi(n)
φ(n) 为欧拉函数
× 2 \times 2 ×2 是因为 ( x , y ) (x,y) (x,y) 为有序点对, − 1 -1 −1 是因为 ( 1 , 1 ) (1,1) (1,1) 会被多算一次
直接线性筛 φ ( n ) \varphi(n) φ(n) 然后前缀和一下就好了
时间复杂度 O ( n ) O(n) O(n) ,记得开 long long \text{long long} long long 哦
代码如下
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define N (int)(1e7+15)
int prime[N],pcnt,phi[N],sum[N];
bool ck[N];
void Euler(int n)
{
phi[1]=1;ck[1]=1;
for(int i=2; i<=n; i++)
{
if(!ck[i])
{
prime[++pcnt]=i;
phi[i]=i-1;
}
for(int j=1; j<=pcnt&&i*prime[j]<=n; j++)
{
int pos=i*prime[j];
ck[pos]=1;
if(i%prime[j])
{
phi[pos]=phi[i]*phi[prime[j]];
}else
{
phi[pos]=phi[i]*prime[j];
break;
}
}
}
}
signed main()
{
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
// freopen("check.in","r",stdin);
// freopen("check.out","w",stdout);
int n,res=0;
cin >> n;
Euler(n);
for(int i=1; i<=n; i++)
sum[i]+=sum[i-1]+phi[i];
for(int i=1; i<=pcnt&&prime[i]<=n; i++)
res+=sum[n/prime[i]]*2-1;
cout << res << endl;
return 0;
}
转载请说明出处