比赛 |
2025.5.4 |
评测结果 |
AAAAAAAAAA |
题目名称 |
GCD |
最终得分 |
100 |
用户昵称 |
会挽弯弓满月 |
运行时间 |
0.820 s |
代码语言 |
C++ |
内存使用 |
39.61 MiB |
提交时间 |
2025-05-04 11:10:06 |
显示代码纯文本
#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e7+10;
ll n,ans,cnt,k;
ll pri[N],sum[N],phi[N];
bool st[N];
void prime(ll x){
ll t;
phi[1]=1;
for(ll i=2;i<=n;i++){
if(!st[i]){
pri[cnt++]=i;
phi[i]=i-1;
}
t=x/i;
for(ll j=0;pri[j]<=t;j++){
st[pri[j]*i]=1;
if(i%pri[j]==0){
phi[pri[j]*i]=phi[i]*pri[j];
break;
}
phi[pri[j]*i]=phi[i]*(pri[j]-1);
}
}
for(int i=1;i<=n;i++){
sum[i]=sum[i-1]+phi[i];
}
return;
}
int main(){
freopen("gcd_prime.in","r",stdin);
freopen("gcd_prime.out","w",stdout);
scanf("%lld",&n);
prime(n);
for(ll i=0;i<cnt;i++){
k=n/pri[i];
if(k>=1) ans+=2*sum[k]-1;
}
printf("%lld",ans);
return 0;
}