zju 1836 Number Puzzle

容斥原理的应用

#include<iostream>
#include <string>
using namespace std;

int a[12];
int b[12];
int n, m;
double total;
int gcd(int x, int y)
...{
    return y==0?x:gcd(y,x%y);
}
int gsm(int x, int y)
...{
    return (x*y)/gcd(x,y);
}
void dfs(int u, int sum, int cnt, int vn)
...{
    int i;
    if(cnt==vn)
    ...{
        total += m/sum;
        return;
    }
    for (i=u; i<n; ++i)
    ...{
        dfs(i+1, gsm(sum,b[i]), cnt+1, vn);
    }
}
int main()
...{
    int i, j, temp;
    double ans;
    while (cin >> n >> m)
    ...{
        memset(a, 0, sizeof(a));
        for (i=0; i<n; ++i)
        ...{
            scanf("%d", &temp);
            a[temp] = 1;
        }
        n = 0;
        for (i=1; i<=10; ++i)
        ...{
            if(a[i] && i<=m)
                b[n++] = i;
        }
        ans = 0;
        for (i=1; i<=n; ++i)
        ...{
            total = 0;
            dfs(0,1,0,i);
            if(i%2!=0)
                ans += total;
            else
                ans -= total;
        }
        printf("%.0lf ",ans);
    }
    return 0;
}