[题解]（gcd/欧拉函数）luogu_P2568_GCD

求gcd(x,y)=p等价于求gcd(x/p,y/p)=1，转化为了n/p内互质的个数

所以欧拉函数，因为有序所以乘2，再特判一下只有在1,1情况下才会重复计算，所以每次都减一

数组开小一时爽，提交wa火葬场！！！

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=;
int n;
int ck[maxn],prime[maxn],phi[maxn],tot;
long long sum[maxn];
void eular(int n){
    phi[]=;
    for(int i=;i<=n;i++){
        if(!ck[i]){
            ck[i]=i,prime[++tot]=i;
            phi[i]=i-;
        }
        for(int j=;j<=tot;j++){
            if(prime[j]>ck[i] || i*prime[j]>n)break;
            ck[prime[j]*i]=prime[j];
            phi[i*prime[j]]=phi[i]*(i%prime[j]?prime[j]-:prime[j]);
        }
    }
}
int main(){
    scanf("%d",&n);
    eular(n);
    for(int i=;i<=n;i++)sum[i]=sum[i-]+phi[i];
    long long ans=;
    for(int i=;i<=tot;i++){
        ans+=*sum[n/prime[i]]-;
    }
    printf("%lld\n",ans);
}