HDU 1796 How many integers can you find 【容斥】

题目：点击打开链接

题意：给出n和集合m，求出所有满足的x，x<n且x能整除m集合中的任意值

分析：考虑到容斥关系，先求出能被任意一个数整除的所有数的数量，减去能被任意两个数整除的所有数的数量，再加上能被任意三个数整除的所有数的数量...即可。注意能被两数整除的数也就是能整除它们最小公倍数的数，集合m中的数可能为0。

 #pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define sqr(x) ((x)*(x))
#define ll long long
const ll INF = 3e18;
const int inf=0x3f3f3f3f;
const int M=100010;
const int N=100100;
const ll MOD=1000000007;
const double eps=1e-8;
const double pi=acos(-1.0);
using namespace std;
int n,m;
ll ans;
vector<ll> a;
ll gcd(ll x,ll y) {return y==0?x:gcd(y,x%y);}
ll lcm(ll x,ll y) {return x/gcd(x,y)*y;}
void dfs(int idx,int num,int flag)
{
    for(int i=idx;i<m;i++)
    {
        ll tmp=lcm(num,a[i]);
        if(flag%2) ans+=n/tmp;
        else ans-=n/tmp;
        dfs(i+1,tmp,flag+1);
    }
}
int main()
{
    ll tmp;
    int i,j;
    while(~scanf("%d%d",&n,&m))
    {
        a.clear();
        n--;ans=0;
        for(i=0;i<m;i++)
        {
            scanf("%I64d",&tmp);
            if(tmp) a.push_back(tmp);
        }
        m=(int)a.size();
        dfs(0,1,1);
        printf("%I64d\n",ans);
    }
}