hdu 3929 Big Coefficients(Lucas定理+容斥原理)

题意：求F(x) = (1+x)^a1 + (1+x)^a2 + ... + (1+x)^am.中系数为奇数的项的个数。

思路：由Lucas定理可知，(1+x)^n中奇数项的个数等于2^(n的二进制表示中1的个数)。Lucas定理相关证明可以看这里。然后就是利用容斥原理，把重复的减去，两项相交的结果中奇数项的个数为2^(a[i]&a[j]的二进制表示中1的个数)。需要注意的是，由于偶数个奇数项在一起并不是奇数项，因此这里要减去。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=20;
ll a[maxn],ans;
int n;
int cal(ll x)
{
    int cnt=0;
    while(x)
    {
        x=x&(x-1);
        cnt++;
    }
    return cnt;
}
void dfs(int p,ll msk,ll v)
{
    ans+=(1LL<<cal(msk))*v;
    for(int i=p+1;i<n;++i)
        dfs(i,msk&a[i],-2*v);
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int t,tcase=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;++i)
            scanf("%I64d",&a[i]);
        ans=0;
        for(int i=0;i<n;++i)
            dfs(i,a[i],1);
        printf("Case #%d: %I64d\n",++tcase,ans);
    }
    return 0;
}