BZOJ_P2818 GCD(数论GCD+欧拉函数线性筛)

BZOJ传送门
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 3081 Solved: 1369
[Submit][Status][Discuss]
Description
给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的
数对(x,y)有多少对.

Input
一个整数N

Output
如题

Sample Input
4

Sample Output
4

HINT
对于样例(2,2),(2,4),(3,3),(4,2)
1<=N<=10^7

Source
湖北省队互测

先丢个链接跑被丢的链接
自己写的题解:

phi用欧拉筛即可，随后求和，计算

#include<cstdio>
#include<iostream>
using namespace std;
#define N 10000005
int n,cnt;long long ans;
long long p[N],phi[N];bool b[N];
int main(){
    scanf("%d",&n);
    for(int i=2;i<=n;i++){
        if(!b[i]){phi[i]=i-1;p[++cnt]=i;}
        for(int j=1;j<=cnt&&i*p[j]<=n;j++){
            b[i*p[j]]=true;
            if(i%p[j]) phi[i*p[j]]=phi[i]*(p[j]-1);
            else{phi[i*p[j]]=phi[i]*p[j];break;}
        }
    }
    phi[1]=1;
    for(int i=2;i<=n;i++) phi[i]+=phi[i-1];
    for(int i=1;i<=cnt;i++) ans+=phi[n/p[i]]*2-1;
    printf("%lld",ans);
    return 0;
}