HDU 1796 How many integers can you find

一个月没做题了。。。

容斥原理，注意处理非互素的数，应用最小公倍数，注意0不用处理，若集合里有1，则直接输出n-1，同时要注意最小公倍数不要溢出，适当的时候应该跳出dfs。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <map>
#include <queue>
#define PB push_back
#define SZ(x) ((int)(x).size())
#define FOR(it,c) for ( __typeof((c).begin()) it=(c).begin(); it!=(c).end(); it++ )
using namespace std;

int n,m,a[20],res;

int gcd(long long a,int b){
    return a%b?gcd(b,a%b):b;
}

void dfs(int now,int k,int d,long long cnt){
    if(cnt>=n) return;
    if(k==d){
        res+=(n-1)/cnt;
        return;
    }
    for(int i=now;i<m;i++)
        if(a[i]&&a[i]<n)
            dfs(i+1,k+1,d,cnt*a[i]/gcd(cnt,a[i]));
}

int main(){
#ifndef ONLINE_JUDGE	          
	freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
#endif
    int ans;
    while(cin>>n>>m){
        ans=0;
        bool ok=false;
        for(int i=0;i<m;i++){
            cin>>a[i];
            if(a[i]==1) ok=true;
        }
        if(ok){
            cout<<n-1<<endl;
            continue;
        }
        int flag=1;
        for(int i=1;i<=m;i++){
            res=0;
            dfs(0,0,i,1);
            ans+=flag*res;
            flag*=-1;
        }
        cout<<ans<<endl;
    }
    return 0;
}