ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Co-prime(容斥原理求互素数个数)(模板)

2021-09-19 21:33:33  阅读:176  来源: 互联网

标签:prime Co int 个数 ll 容斥 long include define


HUD4135

题意

求a~b之间与n互素的个数

思路

求出 1~b 之间与n互素的个数减去 1~a-1之间与n互素的的个数即为所求。

#include<cstdio>
#include<queue>
#include<set>
#include<cstdlib>
#include<string.h>
#include<string>
#include<iostream>
#include<cmath>
#include<unordered_map>
#include<map>
#include<algorithm>
#define endl "\n"
#define IOS ios::sync_with_stdio(0), cin.tie(0),cout.tie(0)
#define ft first
#define sd second
#define pll pair<ll, ll>
#define pii pair<int, int>
#define ll long long int
#define ull unsigned long long int
#define mt(a,b) memset(a, b, sizeof a)
//#define int long long
const double PI = acos(-1.0);
const int inf = 0x3f3f3f3f;
const int INF = 0x7fffffff;
using namespace std;
const int N = 1e5 + 7, M = 1e6 + 10;
ll fac[N], m;

//对n分解质因数,质因数个数为m
void div(ll x)
{
	m = 0;
	int q = sqrt(x);
	for (int i = 2; i <= q; i++)
		if (x % i == 0)
		{
			fac[++m] = i;
			while (x % i == 0) x /= i;
		}

	if (x > 1) fac[++m] = x;
}

//应用容斥原理求出1~b之间含有n的质因子的数的个数,也即与n不互素的数的个数
ll co_prime(ll n)
{
	ll ans = 0;
	for (ll i = 0; i < 1ll << m; i++)//i == 0的时候,t = 1, ans = n,也即是b可以含有任何质因子
	{
		int k = 0;
		ll t = 1;
		for (int j = 0; j < m; j++)
			if (i & 1ll << j) t *= fac[j + 1], k++;
		if (k & 1) ans -= n / t;
		else ans += n / t;
	}

	return ans;
}

int main()
{
	IOS;
	int cas = 1;
	int T; cin >> T;
	while (T--)
	{
		ll a, b, n; cin >> a >> b >> n;
		div(n);

		cout << "Case #" << cas++ << ": ";
		cout << co_prime(b) - co_prime(a - 1) << endl;
	}


	return 0;
}

参考链接

标签:prime,Co,int,个数,ll,容斥,long,include,define
来源: https://blog.csdn.net/QQ2530063577/article/details/120384894

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有