Hackerearth SPECIAL PAIRS DP(SOS).

题意:长度为n的序列a,选两个数(a[i],a[j]) 问a[i] & a[j] == 0的方案数.
n<=1e5,a[i]<=1e6.

假如x=101001 若x&y=0 则y:0?0??0 即y为010110的某个子集.
设F[mask] 序列a中有多少个数为mask的子集. 

dp[mask][i] 类似01背包DP:mask,i位数字的子集个数,第i位和mask相同或者不同

dp[mask[i]=dp[mask][i-1]+dp[mask^(1<<i)][i-1]  O(M*2^M) M位二进制的最大位数.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e5+20,M=21;
int n,a[N];
ll f[N*20];
int main()
{
    int T,x;
    cin>>T;
    while(T--)
    {
        memset(f,0,sizeof(f));
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]),f[a[i]]++;
        for(int i=0;i<=M;i++)
        {
            //for(int mask=0;mask<(1<<M);mask--)
            for(int mask=(1<<M)-1;mask>=0;mask--)
                if(mask&(1<<i))
                    f[mask]+=f[mask^(1<<i)];

        }
        ll ans=0;
        for(int i=1;i<=n;i++)
        {
            //101001,y=0?0??0->010110的某个子集
            int x=a[i],y=(~x)&((1<<(M))-1);
            ans+=f[y];
        }
        printf("%I64d\n",ans);

    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值