poj 1091

最新推荐文章于 2019-10-28 13:43:35 发布

weixin_30871293

最新推荐文章于 2019-10-28 13:43:35 发布

阅读量149

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/zhaozhe/archive/2011/04/19/2020811.html

本文详细介绍了容斥原理，并通过一个具体的编程实现案例来帮助理解该原理的应用场景。通过分解数学问题为编程任务，文章展示了如何利用容斥原理解决特定类型的问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

容斥原理。

代码：

#include<iostream>
#include<fstream>

using namespace std;

long long n,m;
long long num[60],cnt[60];
long long top,per;

long long f(long long ,long long);

void get(long long s,long long t,long long p){
	long long i;
	if(t==p)
	{
		long long j=m;
		for(i=0;i<t;i++)
			j/=cnt[i];
		per+=f(j,n);
	}
	else
	{
		for(i=s;i<top;i++)
		{
			cnt[t]=num[i];
			get(i+1,t+1,p);
		}
	}
}

long long f(long long s,long long t){
	long long i,j=1;
	for(i=1;i<=t;i++)
		j*=s;
	return j;
}

void read(){
//	ifstream cin("in.txt");
	long long i,j,k;
	cin>>n>>m;
	long long ans=0;
	long long x=m;
	for(i=2;i*i<=x;i++)
		if(x%i==0)
		{
			while(x%i==0)
			{
				x/=i;
			}
			num[top++]=i;
		}
	if(x>1)
		num[top++]=x;

	ans=f(m,n);
	for(i=1;i<=top;i++)
	{
		per=0;
		get(0,0,i);
		if(i%2)
			ans-=per;
		else
			ans+=per;
	}
	cout<<ans<<endl;
}



int main(){
	read();
	return 0;
}

转载于:https://www.cnblogs.com/zhaozhe/archive/2011/04/19/2020811.html