hdu 5135 Co-prime(求m以内与n互质的个数)

文章目录

题目链接:

http://acm.hdu.edu.cn/showproblem.php?pid=4135
如果是求某一个n互质的个数的话,那确实是这样容斥比较好
原理也感觉比较简单,要找与n互质那就是要找,那就减去不互质的个数不就行了,于是就把n的质因子找出来,那么质因子的倍数就肯定与n不互质,然后就容斥。

容斥有dfs,和位运算枚举出来,而今天学到了一个直接用队列数组容斥出来的,哇感觉很爽~
参考这位大佬的博客:https://www.cnblogs.com/jiangjing/archive/2013/06/03/3115470.html

#include"bits/stdc++.h"
#define out(x) cout<<#x<<"="<<x
#define C(n,m) (m>n?0:(long long)fac[(n)]*invf[(m)]%MOD*invf[(n)-(m)]%MOD)
using namespace std;
typedef long long LL;
const int maxn=1e5+5;
const int MOD=1e9+7;
LL L,R,N;
int t;
int prime[maxn];
int que[maxn];
LL solve(LL n)//求n以内与N不互质的个数 
{
	LL ans=0;
	for(int i=2; i<=t; i++)
	{
		ans+=n/que[i];
	}
	return ans;
}
int main()
{
	int T;
	cin>>T;
	for(int Case=1; Case<=T; Case++)
	{
		cin>>L>>R>>N;
		LL n=N;
		int cnt=0;
		for(LL i=2; i*i<=n; i++)
		{
			if(n%i==0)
			{
				while(n%i==0)n/=i;
				prime[++cnt]=i;
			}
		}
		if(n>1)prime[++cnt]=n;
		t=0;
		que[++t]=-1;
		for(int i=1; i<=cnt; i++)//用队列数组进行容斥 
		{
			int len=t;
			for(int k=1; k<=len; k++)que[++t]=que[k]*prime[i]*(-1);
		}
		cout<<"Case #"<<Case<<": "<<R-L+1-(solve(R)-solve(L-1))<<endl;

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值