hdu 4135 求[l,r]之间与n互质的数的个数(容斥原理)

先把n素因子分解,再用位运算,枚举倍数,加上容斥原理

#include<stdio.h>
#include<vector>
using namespace std;
#define ll long long
ll work(ll r, ll n){
	vector<ll> p;
	for(ll i = 2; i * i <= n; i++){
		if(n % i==0){
			p.push_back(i);
			while(n%i==0) n /= i;
		}
	}
	if(n>1) p.push_back(n);
	ll ans = 0;
	for(ll m = 1; m < (1<<p.size()); m++){
		ll mul = 1, bits = 0;
		for(ll i = 0; i < p.size(); i++){
			if(m & (1<<i)){	
				++bits;
				mul *= p[i];
			}
		}
		ll cur = r / mul;
		if(bits & 1) ans += cur;
		else ans -= cur;
	}
	return r - ans;
}
int main()
{
	//freopen("in.txt","r",stdin);
	int t;
	ll n, l, r;
	scanf("%d", &t);
	int ca = 1;
	while(t--)
	{
		scanf("%lld%lld%lld", &l, &r, &n);
		printf("Case #%d: %lld\n", ca++, work(r, n) - work(l-1, n));
	}
	return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值