容斥原理的应用

三个集合
S = A + B + C - AB - AC - BC + ABC
四个集合
S = A + B + C + D - AB - AC - AD -BC - BD - CD + ABC + ABD + ACD + BCD - ABCD
则n个集合
S = 1个集合 - 2个集合交集 + 3个集合交集 + … + (-1) ^ (n+1) * n个集合交集。共有(2 ^ n) - 1项

证明：Cn0+Cn1+Cn2+...+Cnn=2^n 则Cn1+Cn2+...+Cnn=2^n-1

暴力解法 

#include<bits/stdc++.h>
using namespace std;
const int N=20;
int p[N];
int n,m;
int res=0;
int main(){
	cin>>n>>m;
	for(int i=0;i<m;i++) cin>>p[i];
	for(int i=1;i<=n;i++){
		for(int j=0;j<m;j++){
			if(i%p[j]==0){
				res++;
				break;
			}
		}
	}
	cout<<res;
	return 0;
}

解析点这！！！

#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const int N=20;
int p[N];//存质数的数组
int n,m,res=0;//res用来存答案

int main(){
	scanf("%d%d",&n,&m);
	for(int i=0;i<m;i++) scanf("%d",&p[i]);
	for(int i=1;i<(1<<m);i++){//把i转化为m位的二进制数  2的m次方可表示为1<<m  
	                          //由容斥原理可得共有2的m次方-1种选法 (m个数 每个数选或不选)
		int t=1,cnt=0;        //t表示所选数的乘积  cnt表示该二进制数中1的个数
		for(int j=0;j<m;j++){  //每个数共有m位
			if((i>>j)&1){     //这里不是每次移动一位 直接移动j位
				cnt++;        //最后一位是1则加1
				if((LL)t*p[j]>n){  //乘积的范围大于n n的范围内不存在该数倍数的数
					t=-1;
					break;//把该数跳出
				}
				t=t*p[j];
			}
		} 
		if(t!=-1){
			if(cnt%2) res+=n/t;//  n/t代表n内包含t的倍数的个数 1的个数为奇数则加
			else res-=n/t;
		}
	}
	cout<<res<<endl;
	return 0;
}