容斥理论

输入N和M,M表示接下来会输入M个数,求从1到N中不能被这M个数整出的数的个数。

思路:先算出从1到N中能被这M个数整除(即是这M个数中任意一个数的倍数)的数的个数,在用N减去就是答案了。

#include <iostream>
#include <string.h>
using namespace std;
long long num[20];

long long gcd(long long a,long long b) //求a和b的最大公约数
{
    return b ? gcd(b,a%b) : a;//b为0时返回a,b大于0时返回gcd(b,a%b)
}

int main()
{
    long long n;
    int m;
    while(cin>>n>>m)
    {
        for(int i=0; i<m; i++)
        cin>>num[i];
        int ans=0;
        for(int i=1;i<(1<<m);i++)
        {
            int b=0;
            long long w=1;
            for(int j=0;j<m;j++)
            if(i&(1<<j))
               {
                  b++;
                  w=w/gcd(num[j],w)*num[j];//表示w和num[j]的最小公倍数
               }
            if(b&1)//表示如果b是奇数,就让最后结果加上所得个数
            ans+=n/w;
            else ans-=n/w;
        }
        cout<<n-ans<<endl;
    }
    return 0;
}

posted @ 2015-08-02 17:37  茶飘香~  阅读(172)  评论(0编辑  收藏  举报