容斥原理

容斥原理：所有情况相加-重复的情况。

hdu1796
推荐指数：****
这是最最基础的模板
题意：给n和m，m个数，在小于n的数中，有多少为m个数的倍数。
思路：m个数的倍数-重复计入的部分，主要是理解容斥原理。

//二进制 
#include<cstdio>
#include<iostream>
using namespace std;

int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);
}
int arr[25];

int main()
{
	int n,m;
	while(scanf("%d %d",&n,&m)!=EOF)
	{
		int num=0;
		int ans=0;
		for(int i=1;i<=m;i++)
		{
			int ai;
			scanf("%d",&ai);
			if(ai!=0)//数据好坑，还可能为0
				arr[num++]=ai;
		}
		//num种情况，对应二进制num位（0~num-1）
		for(int i=1;i<(1<<num);i++)
		{
			int lcm=1,sum=0;
			for(int j=0;j<num;j++)
			{
				if(i&(1<<j))
				{
					sum++;//位上为一的个数
					lcm=lcm/gcd(lcm,arr[j])*arr[j];//最小公倍数
				}
			}
			if(sum%2==0)ans-=(n-1)/lcm;//偶减
			else ans+=(n-1)/lcm;//奇加
		}
		printf("%d\n",ans);
	}
	
	return 0;
}

//dfs
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int arr[25];
int n,m,num,ans;
int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);
}
void dfs(int a,int lcm,int id)
{
	lcm=arr[a]/gcd(lcm,arr[a])*lcm;
	if(id%2==0)
		ans-=(n-1)/lcm;
	else
		ans+=(n-1)/lcm;
	//深搜,eg:A->AB->ABC,B->BC,C	
	for(int i=a+1;i<num;i++)
		dfs(i,lcm,id+1); 
}

int main()
{
	while(scanf("%d %d",&n,&m)!=EOF)
	{
		 num=0,ans=0;
		 for(int i=1;i<=m;i++)
		 {
		 	int ai;//注意0这种情况 
		 	cin>>ai;
			if(ai!=0)
			arr[num++]=ai;
		 }
		 
		 for(int i=0;i<num;i++)
		 	dfs(i,arr[i],1);
		 	
		 cout<<ans<<endl;
	}
	
	return 0;
}

hdu 4135
推荐指数***

题意：给定区间[a,b]，n，求在[a,b]中有多少个和n互质的数；
思路：[a,b]总数-与n不互质的数=ans；求出n的质因子，利用容斥原理，求出质因子在[a,b]的倍数个数-重复部分；

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;

ll num;
ll arr[20];
void init(ll n)//分解质因子+二进制化
{

   	num=0;
   	for(ll i=2;i*i<=n;i++)
   	if(n%i==0)
   	{
   		while(n%i==0)n/=i;
   		arr[num++]=i;
   	}
   	if(n>1)arr[num++]=n;
}
ll solve(ll n)//
{
   	ll ans=0;
   	for(ll i=1;i<(1ll<<num);i++)
   	{
   		ll l=1,s=0;
   		for(ll j=0;j<num;j++)
   		if(i>>j&1)
   		{
			s++;//该为为1
   			l*=arr[j];//已是质数，不用gcd
   		}
   		if(s&1)ans+=n/l;
   		else ans-=n/l;
   	}
   	return ans;
}
int main()
{
   	int t;
   	ll a,b,n;
   	scanf("%d",&t);
   	int cas=0;
   	while(t--)
   	{
   		scanf("%lld%lld%lld",&a,&b,&n);
   		init(n);
   		printf("Case #%d: %lld\n",++cas,b-a+1-solve(b)+solve(a-1));
   	}
   	return 0;
}

hdu6286
推荐指数***

题意：给两个区间[a,b] [c,d]，求有多少个（x，y）存在(a<=x<=b,c<=y<=d),使xy=k*2018（k是整数）；
思路：2018的因子只有1 2 1009 2018；只需判断在[a,b]区间其因子个数乘对应另一因子在[c,d]的个数，再减去重复计算的。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>

using namespace std;
typedef long long ll;

 //在当前[l,r]区间x因子的个数
ll fac(int l,int r,int x)
{
	return r/x-(l-1)/x;
}
int main()
{
	int l1,r1,l2,r2;
	while(~scanf("%d%d%d%d",&l1,&r1,&l2,&r2))
	{
		ll ans=0;
		ans+=(fac(l1,r1,2018)*(r2-l2+1));
		ans+=(fac(l2,r2,2018)*(r1-l1+1));
		ans-=(fac(l1,r1,2018)*fac(l2,r2,2018));
		ans+=((fac(l1,r1,1009)-fac(l1,r1,2018))*(fac(l2,r2,2)-fac(l2,r2,2018)));
		ans+=((fac(l2,r2,1009)-fac(l2,r2,2018))*(fac(l1,r1,2)-fac(l1,r1,2018)));
		printf("%lld\n",ans);
	}
	return 0;
}

hdu2204
推荐指数：*****
题意：给你一个数n，求[1,n]之间有多少数可以表示为m^k，eg：8=2³；64=8²=2⁶;可见求出质数的次方包含了合数的情况。
下班了下班了 明天写
上线了理解借鉴

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
int prime[30];
int main()
{
	//求质数：10^18约等于2^64，求100以内即可
	int num=0;
	prime[0]=2;
	for(int i=3;i<100;i++)
	{
		int flag=1;
		for(int j=2;j<i;j++)
		{
			if(i%j==0)
			{
				flag=0;
				break;
			}
		}
		if(flag)prime[++num]=i;
	}
	//2,3,5次方相当于容斥定理的三个集合；
	ll n;
	while(cin>>n)
	{
		ll ans=0;
		ll temp;
		for(int i=0;;i++)
		{
			temp=pow(n,1.0/prime[i]);
			if(temp<2)break;
			ans+=temp-1;
			//eg：10,10^(1/2)=3;所以1,2，3的平方都小于10，减去1；（1单独考虑）
			
			for(int j=i+1;;j++)
			{
				temp=pow(n,1.0/prime[i]/prime[j]);
				if(temp<2)break;
				ans-=temp-1;
				
				for(int k=j+1;;k++)
				{
					temp=pow(n,1.0/prime[i]/prime[j]/prime[k]);
					if(temp<2)break;
					ans+=temp-1;
				}
			}
		}//终于懂了。第一轮:开2次、2*3次、2*3*5次；第二轮:开3次、3*5次、3*5*7次；所以素数存到>64；
		cout<<ans+1<<'\n'; 
	}
	return 0;
}