容斥原理应用

hdu 1796

题目意思：求出在1到N之间能背集合M中的元素整除的个数

这是一道简单的容斥原理应用题。先找出1...n内能被集合中任意一个元素整除的个数，再减去能被集合中任意两个整除的个数，即能被它们两只的最小公倍数整除的个数，因为这部分被计算了两次，然后又加上三个时候的个数，然后又减去四个时候的倍数...。用dfs得到集合M中所有元素的组合，从0到n-1枚举当前数的所有组合情况，组合数为奇的加，为偶的减。题目要注意除去M集合中为0的数，它对最后的结果没影响，但会对gcd造成影响。

#include <stdio.h>

__int64 n,m,num,t,a[15];

__int64 gcd(__int64 a,__int64 b)
{
    return b==0 ? a : gcd(b,a%b);
}

void dfs(__int64 beg,__int64 flag,__int64 k)
{
    k = k*a[beg]/gcd(k,a[beg]);
    if(flag==1) num+=n/k;
    else num-=n/k;
    for(__int64 j = beg+1; j < t; j++)
    {
        dfs(j,!flag,k);
    }
}

int main()
{
    __int64 i;
    while(~scanf("%I64d%I64d",&n,&m))
    {
        num = 0;
        __int64 nu;
        t = 0;
        for(i = 0; i < m; i++)
        {
            scanf("%I64d",&nu);
            if(nu!=0) a[t++] = nu;
        }
        n--;
        for(i = 0; i < t; i++)
        dfs(i,1,a[i]);
        printf("%I64d\n",num);
    }
    return 0;
}