版权声明:大佬您能赏脸,蒟蒻倍感荣幸,还请联系我让我好好膜拜。 https://blog.csdn.net/ShadyPi/article/details/82055566
暂无链接
第一题
【问题描述】
给定 n,求所有小于等于 n 的数,与 n 的最大公约数之和。
【输入格式】
输入包含一个非负整数 n。
【输出格式】
输出一个整数,为求和的结果。
【输入样例 1】
10
【输出样例 1】
27
【样例说明 1】
有四个数和 10 的最大公约数为 1,四个数为 2,一个数为 5,一个数为 10。
【输入样例 2】
123456
【输出样例 2】
1644800
【数据范围】
对于 40%的数据点, 。
对于 60%的数据点, 。
对于 100%的数据点, 。
题解
我特么为什么信了 的邪要去推式子。。。
考虑 肯定是 的约数,那么不妨枚举 的约数来统计每个约数的贡献。
因为 ,就有 ,即 出现了所有小于等于 且与 互质的数的个数次,所以 的贡献就是 。
注意特判 为完全平方数的情况。
代码
#include<bits/stdc++.h>
using namespace std;
const int M=1e5+5;
int mx,n;
long long ans;
int phi(int x)
{
if(x==1)return 1;
int sqt=sqrt(x),ans=1,flag=0;
for(int i=2;i<=sqt;++i){for(;x%i==0;x/=i,ans*=i,flag=1);if(flag)ans=ans/i*(i-1),flag=0;}
if(x>sqt)return ans*(x-1);
return ans;
}
void in(){scanf("%d",&n);}
void ac()
{
mx=sqrt(n);
for(int i=1;i<=mx;++i)if(n%i==0)ans+=1ll*i*phi(n/i)+1ll*(n/i)*phi(i);
if(mx*mx==n)ans-=1ll*mx*phi(n/mx);
printf("%lld",ans);
}
main(){in();ac();}