比赛 2025.5.4 评测结果 AAAAAAAAAA
题目名称 GCD 最终得分 100
用户昵称 彭欣越 运行时间 0.688 s
代码语言 C++ 内存使用 53.76 MiB
提交时间 2025-05-04 10:55:20
显示代码纯文本
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e7+5;
ll n,f[N],p[N],phi[N],s[N],cnt;
ll ans;
void judge (int n) {
	phi[1]=1;
	for (int i=2;i<=n;i++) {
		if (!f[i]) p[++cnt]=i,phi[i]=i-1;
	    for (int j=1;j<=cnt&&i*p[j]<=n;j++) {
	    	f[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++) s[i]=s[i-1]+phi[i];
} 
int main () {
	freopen("gcd_prime.in","r",stdin);
	freopen("gcd_prime.out","w",stdout);
	cin >> n;
	judge(n);
	for (int i=1;i<=cnt;i++) {
		ans+=2*s[n/p[i]]-1;
	}
	cout << ans <<endl;
	return 0;
}