codeforces449D&&51nod 1407 与与与与 [dp][容斥原理]

中文题面

题解

1/10 UPD:经dalao提醒,原来我的代码里对f[x]有取模操作,这在后面的计算中显然会造成影响,而因为本题f[x]的上限不会达到被取模的规模,跪谢chesterking

设$f[x]$为&x==x的数的个数,$g[x]$为x转化为二进制以后1的个数

转化一下,求等于0的组数即 所有组合情况-含有位为1的组合数

容斥求一下含有位为1的组合数

$$\sum^{2^{20}}_{x=1}=(-1)^{g[x]-1}*(2^{f[x]}-1)$$

2的幂减一是因为不能一个都不选,最后用所有情况($2^n$)减一下就好了

答案就是

$$\sum^{2^{20}}_{x=1}=(-1)^{g[x]}*(2^{f[x]}-1)$$

至于求$f[x]$可以做一下51nod1406

1406题解

#include<cstdio>
#include<algorithm>
#include<cstring>
#define P 1000000007
#define LL long long
using namespace std;

int n,g[1<<20],ans;
LL f[1<<20],pw[21];

inline int read(){
    int a=0;char f=1,c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){a=a*10+c-'0';c=getchar();}
    return a*f;
}

LL ksm(LL a,LL k){
    LL tmp=1;
    while(k){
        if(k&1) tmp=tmp*a%P;
        a=a*a%P;k>>=1;
    }
    return tmp;
}

int main(){
    pw[0]=1;
    for(int i=0;i<=20;++i){
        if(i) pw[i]=pw[i-1]<<1;
        for(int j=0;j<(1<<20);++j)
            if(j&pw[i]) ++g[j];
    }
    while(scanf("%d",&n)!=EOF){
        memset(f,0,sizeof(f));
        for(int i=1,x;i<=n;++i){
            x=read();
            ++f[x];
        }
        for(int i=20;i>=0;--i)
            for(int j=0;j<(1<<20);++j)
                if(j&pw[i]) f[j^pw[i]]=f[j^pw[i]]+f[j];
        ans=0;
        for(int i=0;i<(1<<20);++i){
            int opt=(g[i]&1)?-1:1;
            ans=(ans+(opt*((ksm(2,f[i])-1+P)%P)%P)+P)%P;
        }
        printf("%lld\n",ans);
    }
    return 0;
}