[BZOJ 2818] Gcd 线性筛+欧拉函数前缀和

题目传送门：【BZOJ 2818】

---

题目大意：输入整数 N，求 1 ≤ x , y ≤ N 且 gcd ( x , y ) 为质数的数对 ( x , y ) 的数目。其中 1 ≤ N ≤ $10^7$。

---

题目分析：

由题，欲求 gcd ( x , y ) = p（1 ≤ x , y ≤ N，p 为质数），根据 gcd 的性质，
我们可将其转化一下，变成：求 gcd ( x , y ) = 1（1 ≤ x , y ≤ $\lfloor \frac{N}{p} \rfloor$）。

令 $\lfloor \frac{N}{p} \rfloor$ = d，设此时的答案为 ans$_{pi}$；在 1 ≤ x , y ≤ $\lfloor \frac{N}{p} \rfloor$ 这个范围内，根据欧拉函数 φ 的定义（与这个数互质的数的数量），可知：ans$_d$ = 2 * ( φ(d) + φ(d-1) + φ(d-2) + …… + φ(2) + φ(1) ) - 1。

对于上式，由于在数对 ( x , y ) 中，x 和 y 可以交换，因此答案要乘以 2；又因为对于 ( 1 , 1 ) 只有一种情况，所以最后的答案要减 1。

之后，我们对每个不大于 N 的质数进行类似的处理，最后得到的答案 ans = $\sum ans_{pi}$（所有的 $p_i$ 不大于 N）。

因此，我们只需要预处理出 $10^7$以内的所有质数，所有数的 φ 值，以及 φ 值的前缀和，枚举每个质数 $p_i$ 并求出对应的 $ans_{pi}$，最后求和即可（前缀和&最后的答案要开 long long）。

---

下面附上代码：

[cpp] view plain copy

print ?

1. #include<cstdio>  
2. const int MX=10000005;  
3.   
4. int n,ptot=0,prime[MX],phi[MX];                 //ptot:质数总数   
5. bool isnot[MX];   
6. long long ans=0,sum_phi[MX];                    //sum_phi: φ的前缀和   
7.   
8. void sieve(int n){  
9.     isnot[1]=true;  
10.     phi[1]=1;  
11.     for (register int i=2;i<=n;i++){  
12.         if (!isnot[i]){  
13.             prime[++ptot]=i;  
14.             phi[i]=i-1;  
15.         }  
16.         for (int t=1;t<=ptot;t++){  
17.             int j=prime[t]*i;  
18.             if (j>n) break;  
19.             isnot[j]=true;  
20.             phi[j]=phi[prime[t]]*phi[i];  
21.             if (i%prime[t]==0){  
22.                 phi[j]=prime[t]*phi[i];  
23.                 break;  
24.             }  
25.         }  
26.     }  
27. }  
28.   
29. int main(){  
30.     scanf(”%d”,&n);  
31.     sieve(n);  
32.     for (register int i=1;i<=n;i++)              //预处理 φ的前缀和   
33.         sum_phi[i]=sum_phi[i-1]+phi[i];  
34.           
35.     for (register int i=1;i<=ptot;i++){          //枚举每个质数 pi并求出此时的 ans   
36.         int d=n/prime[i];  
37.         ans+=sum_phi[d]*2-1;  
38.     }  
39.     printf(”%lld”,ans);  
40.     return 0;  
41. }