HDU2588 数论 欧拉函数

原题:http://acm.hdu.edu.cn/showproblem.php?pid=2588

题解:求1<=X<=N and (X,N)>=M.的个数。题目中的n很大我么可以考虑枚举gcd。

显然对于每个gcd都应为N的因数。

gcd(a,N)=gg|N

变形下的:

gcd(\frac{a}{g},\frac{N}{g})=1

对于gcd为g的a的个数显然是\varphi(\frac{N}{g})

 

#include<cstdio>
#include<cstring>
#ifdef WIN32
#define LLD "%I64d"
#else
#define LLD "%lld"
#endif
using namespace std;

int n,m,cas,ans; 
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("gcd.in","r",stdin);
	scanf("%d",&cas);
	while(cas--){
		scanf("%d%d",&n,&m); ans=0;
		for(int i=1;1ll*i*i<=1ll*n;i++){
			if(n%i!=0) continue;
			if(i>=m) ans+=phi(n/i);
			if(n/i>=m && i*i!=n) ans+=phi(i);
		}
		printf("%d\n",ans);	
	}
	return 0;
}

 

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

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值