记录编号	600451	评测结果	AAAAAAAAAA
题目名称	[BZOJ 2818]GCD	最终得分	100
用户昵称	Ruyi	是否通过	通过
代码语言	C++	运行时间	0.480 s
提交时间	2025-05-04 15:02:58	内存使用	30.59 MiB

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,cnt,p[10000001],phi[10000001];
ll sum[10000001],ans;
bool flag[10000001];
void f(int n){
	phi[1]=1;
	for(int i=2;i<=n;++i){
		if(!flag[i]){
			p[++cnt]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=cnt&&i*p[j]<=n;++j){
			flag[i*p[j]]=1;
			if(i%p[j]==0){
				phi[i*p[j]]=phi[i]*p[j];
				break;
			}else{
				phi[i*p[j]]=phi[i]*phi[p[j]];
			}
		}
	}
	for(int i=1;i<=n;++i) sum[i]=sum[i-1]+phi[i];
}
int main(){
	freopen("gcd_prime.in","r",stdin);
	freopen("gcd_prime.out","w",stdout);
	cin>>n;
	f(n);
	for(int i=1;i<=cnt;++i) ans+=2*sum[n/p[i]]-1;
	cout<<ans<<endl;
	return 0;
}