容斥原理

容斥原理 ：重要应用 求出一个数n在区间[1,m]里面有多少个数与它互质。假设数据不超过int型。

 

实现过程分为两步：

1, 求出m的质因子 并保存在数组里面；

2, 求出区间[1,n]里面有多少个数与m不互质。

#include <cstdio>
#include <cmath>
int p[10];//保存质因子  int型n不会超过10个 
int k;//记录质因子个数
void getp(int n)//求出n的质因子 
{
	int i;
	k = 0;//初始化
	for(i = 2; i*i <= n; i++)
	{
		if(n % i == 0)
		{
			p[k++] = i;//保存质因子 
			while(n % i == 0)
			n /= i;
		}
	}
	if(n > 1) p[k++] = n;//本身是质数 
} 
int nop(int m)//求出区间[1,m]里面有多少个数与n不互质 
{
	int top = 0;//队列顶点 
	int que[10100];
	int i, j, t;
	que[top++] = -1;//队列数组保存n所有质因子任意不相同组合的乘积 
	for(i = 0; i < k; i++)
	{
		t = top;//利于下面计算  
		for(j = 0; j < t; j++)
		{
			que[top++] = que[j] * p[i] * (-1);//奇加偶减 
		}
	}
	int sum = 0;//统计个数
	for(i = 1; i < top; i++)
	sum += m / que[i];
	return sum;
}
int main()
{
	int n, m;
	while(scanf("%d%d", &n, &m), n||m)//求区间[1,m]内有多少个数与n互质 
	{
		getp(n);
		printf("%d\n", m-nop(m));
	}
	return 0;
}

dfs版本

#include <cstdio>
#include <cmath>
int p[10];//保存质因子  int型n不会超过10个 
int k;//记录质因子个数
void getp(int n)//求出n的质因子 
{
	int i;
	k = 0;//初始化
	for(i = 2; i*i <= n; i++)
	{
		if(n % i == 0)
		{
			p[k++] = i;//保存质因子 
			while(n % i == 0)
			n /= i;
		}
	}
	if(n > 1) p[k++] = n;//本身是质数 
} 
int nop(int m, int t)//求出区间[1,m]里面有多少个数与n不互质 
{
	int i, sum = 0;
	for(i = t; i < k; i++)
	sum += m / p[i] - nop(m/p[i],i+1);
	return sum;
}
int main()
{
	int n, m;
	while(scanf("%d%d", &n, &m), n||m)//求区间[1,m]内有多少个数与n互质 
	{
		getp(n);
		printf("%d\n", m-nop(m, 0));
	}
	return 0;
}