容斥原理概况:
用二进制表示状态的小技巧非常常用,后面的状态压缩DP也用到了这个技巧,因此一定要掌握
。
能被整除的数
给定一个整数 nn 和 mm 个不同的质数 p1,p2,…,pmp1,p2,…,pm。
请你求出 1∼n1∼n 中能被 p1,p2,…,pmp1,p2,…,pm 中的至少一个数整除的整数有多少个。
输入格式
第一行包含整数 nn 和 mm。
第二行包含 mm 个质数。
输出格式
输出一个整数,表示满足条件的整数的个数。
数据范围
1≤m≤161≤m≤16,
1≤n,pi≤1091≤n,pi≤109输入样例:
10 2 2 3
输出样例:
7
#include<iostream>
using namespace std;
typedef long long ll;
const int N = 20;
int p[N], n, m;
int main()
{
cin >> n >> m;
for(int i = 0; i < m; i++) cin >> p[i];
int res = 0;
枚举从1 到 1111....(m个 1) 的没一个状态集合,(至少选中一个集合)
for(int i = 1; i < 1 << m; i ++)
{
int t = 1; 选中集合对应质数的乘积
int s = 0; 选中的集合数量
枚举当前状态的每一位
for(int j = 0; j < m; j++)
{
if(i >> j & 1)
{
if((ll)t * p[j] > n) 乘积大于 n, 则 n / t = 0,跳出这次循环
{
t = -1;
break;
}
s ++ ; 有一个状态为 1 ,集合数量加 1
t *= p[j];
}
if( t == -1 ) continue;
if(s & 1) res += n / t; 容斥原理,奇数个相加,偶数个相减
else res -= n / t;
}
}
cout << res << endl;
return 0;
}