题意:给定一个整数 n 和 m 个不同的质数 p1,p2,…,pm。
1≤m≤16,
1≤n,pi≤109
请你求出 1∼n 中能被 p1,p2,…,pm 中的至少一个数整除的整数有多少个。
能被整除的数
题目分析:
先考虑朴素做法:
1.对于每个1-n每个数试除法判断,复杂度O(nm)一定会超时
2.打表枚举p1-pm的倍数,开个布尔数组记录是否存在,再历遍1-n求得个数,bool数组开不到那么大(可恶).
接着引入今天的主题-容斥定理
用人话说就是,当选取集合交集的个数为奇数时符号为正,选取集合符号为负.我们举个例子当n=10 p={2,3}
我们的答案等于n/2+n/3-n/(23)(下取整),也就是S1+S2-S12,Si的含义是1-n中能被质数pi整除的个数,由于元素互质
推广成3个集合即为S1+S2+S3-S12-S13-S23+S123
证明:
求三圆面积:S表示圆的面积,S12表示其重合面积,很容易得到
S总=S1+S2+S3-S12-S13-S23+S123
我们接着讲题对于每个pi我们可以知道其在1-n有多少个倍数,那么这里的集合S为pi的所有倍数,我们想知道S的交集,那么我们枚举所有组合的情况求和即可,这里给上代码配合讲解
#include<bits/stdc++.h>
using namespace std;
int p[20];
int main()
{ int n,m;
cin>>n>>m;
for(int i=0;i<m;i++)
cin>>p[i];
int ans=0;
for(int i=1;i<1<<m;i++)//用二进制数来表示方案选取情况,0为不选第j位,1为选取第j位
{ int t=1,cnt=0;
for(int j=0;j<m;j++)
{
if(i>>j&1)//第j位上是否为1
{
cnt++;
if((long long)t*p[j]>n){t=-1;break;}
t*=p[j];//由于互质求并集就是将他们相乘
}
}
if(t!=-1)
{
if(cnt%2)ans+=n/t;
else ans-=n/t;
}
}
cout<<ans<<endl;
return 0;
}