poj2480 数论

原题:poj.org/problem?id=2480

题目:求\sum gcd(i,N) 1<=i<=N。可以枚举gcd即:gcd(x,N)=g

变形一下:

gcd(\frac{x}{g},\frac{N}{g})=1,这样就转化为:\varphi(\frac{N}{g})个gcd为g的数

答案为:

\sum_{g|x} \varphi(\frac{N}{g})*g

#include<cstdio>
#include<cstring>
#include<algorithm>
#ifdef WIN32
#define LLD "%I64d"
#else
#define LLD "%lld"
#endif
#define ll long long
using namespace std;
ll n;
inline int phi(int x){
	int ans=x;
	for(int i=2;1ll*i*i<=1ll*x;i++){
		if(x%i==0) {
			ans=ans/i*(i-1);
			while(x%i==0) x/=i;
		}
	}
	if(x>1) ans=ans/x*(x-1);
	return ans;
}
int main(){
//	 freopen("poj2480.in","r",stdin);
	 while(scanf(LLD,&n)!=EOF){
	 	ll ans=0;
	 	for(int i=1;1ll*i*i<=1ll*n;i++){
			if(n%i==0){
				int x=phi(n/i);ans+=1ll*x*i;
				if(i*i!=n) {
					int y=phi(i);ans+=1ll*y*(n/i);
				} 	
			}
		 }
		printf(LLD"\n",ans);
	}
	return 0;
}

 

 

转载于:https://www.cnblogs.com/Exception2017/p/10252090.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值