BZoj 2705: [SDOI2012]Longge的问题【数论】
来源:互联网 发布:怎么在淘宝上搜115 编辑:程序博客网 时间:2024/06/05 13:22
Longge的问题
Time Limit:3000MS Memory Limit:131072KB 64bit IO Format:%lld & %lluDescription
Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。
Input
一个整数,为N。
Output
一个整数,为所求的答案。
Sample Input
6
Sample Output
15
Hint
【数据范围】
对于60%的数据,0<N<=2^16。
对于100%的数据,0<N<=2^32。
题解:
最直观的想法肯定是要枚举n 的所有因子,然后想办法统计出来所有的和 n 的最大公约数正好是这个因子的数的个数,之后乘上这个因子,最后累加,就是需要的答案了,但是问题来了...
问题就在于:如何统计每个因子(设为x)对应的满足 gcd(n,m)= x 的m 的个数,下面是数学分析:
需要求出满足 gcd(n,m)= x 的m 的个数,不难发现,n 和m 的最大公约数为x,那么就可以同时除去x
相当于求满足 gcd(n/x,m/x)=1 的数的个数,额...这不就是欧拉函数吗........
然后就是相乘再累加了.....
明白到这里之后,就能很简单的编程实现了....
另外,枚举n 的因子,只需要枚举到根号n,但是要特别注意 n 是完全平方数的时候,切忌重复统计!
ps: 数学不好的搞ACM真的伤不起啊,一见到数学题就卡壳了,搞了半天感觉搞懂了,但是下一次见到的时候让然还是不会做..............
#include<stdio.h>#include<math.h>typedef long long ll;ll euler(ll x)//欧拉函数{ll ans=x,tp=sqrt(x*1.0);for(ll i=2;i<=tp;++i){if(x%i==0){ans=ans-ans/i;while(x%i==0){x/=i;}}}if(x>1){ans=ans-ans/x;}return ans;}int main(){ll n;while(~scanf("%d",&n)){ll ans=0,tp=sqrt(n*1.0);for(ll i=1;i<=tp;++i){if(n%i==0){ans+=i*euler(n/i)+n/i*euler(i);}}if(tp*tp==n) //注意这个{ans-=tp*euler(tp);}printf("%lld\n",ans);}}
1 0
- BZoj 2705: [SDOI2012]Longge的问题【数论】
- BZOJ 2705 [SDOI2012]Longge的问题 数论
- BZOJ 2705 [SDOI2012]Longge的问题 题解
- 【BZOJ 2705】 [SDOI2012]Longge的问题
- ❤BZOJ 2705: [SDOI2012]Longge的问题
- bzoj 2705: [SDOI2012]Longge的问题
- Longge的问题 [Bzoj 2705,SDOI2012]
- bzoj 2705: [SDOI2012]Longge的问题
- BZOJ 2705 [SDOI2012]Longge的问题
- BZOJ 2705 [SDOI2012]Longge的问题
- bzoj 2705 [SDOI2012]Longge的问题
- 【bzoj 2705】[SDOI2012]Longge的问题
- 新视野OJ 2705 [SDOI2012]Longge的问题 (数论)
- 【SDOI2012】【BZOJ2705】【Longge的问题】【题解】【数论】
- [BZOJ2705][SDOI2012]Longge的问题(数论)
- [杂题 数论] BZOJ2705: [SDOI2012]Longge的问题
- bzoj2705: [SDOI2012]Longge的问题(数论)
- BZOJ 2705: [SDOI2012]Longge的问题(euler函数)
- 支付宝接口
- C++ 和java的不同点
- Maven环境的搭建
- Android 事件分发机制
- 关于tp框架中ueditor生成乱码问题
- BZoj 2705: [SDOI2012]Longge的问题【数论】
- Spring Security(02)——关于登录
- MySQL建外键简直吓死人error150
- caffe代码阅读7:LayerRegistry的实现细节-2016.3.18
- 判断字符串中是否是数字
- 二维码使用之二维码扫描
- LEFT JOIN 与 RIGHT JOIN 的使用区别
- Crawling
- Android Studio 之 gradle 理解