容斥定理-能被整除的数

题意:给定一个整数 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;
}