记录编号 600722 评测结果 AAAAAAAAAA
题目名称 [河南省队2012] 最大公约数和 最终得分 100
用户昵称 Gravatarxxz 是否通过 通过
代码语言 C++ 运行时间 0.030 s
提交时间 2025-05-12 20:06:29 内存使用 3.84 MiB
显示代码纯文本
#include<bits/stdc++.h>
#define ll long long

ll n,t,ans,prime[200000];
bool flag[1000000];
long long ou(ll k){
	ll i,s=k,tmp=k;
	for(i=1;i<=t&&prime[i]*prime[i]<=k;i++){
		if(k%prime[i]==0)s=s/prime[i]*(prime[i]-1);
		while(k%prime[i]==0)k/=prime[i];
	}
	if(k>1)s=s/k*(k-1);
	return s;
}

int main(){
	freopen("gcdsum.in","r",stdin);freopen("gcdsum.out","w",stdout);
	scanf("%lld",&n);
    t=0;
	for(ll i=2;i*i<=n;i++){
		if(!flag[i]){
			prime[++t]=i;
			for(ll j=i+i;j*j<=n;j+=i){
				flag[j]=1;
			}
		}
	}
	for(ll i=1;i*i<=n;i++){
		if(n%i==0){
			if(i*i==n){
                ans+=ou(i)*i;
                break;
            }else{
                ans+=ou(i)*n/i+ou(n/i)*i;
            }
		}
	}
	printf("%lld\n",ans);
	return 0;
}