hdu 5119 Happy Matt Friends 2014ACM/ICPC亚洲区北京站-重现赛

这一题求得是方案总数，所以和平时写的 f[i][sum]=max(f[i-1][sum],f[i-1][k[i]^sum]+1);的背包不太一样……吧？

这里面的状态转移方程是f[i][sum]=f[i-1][sum]+f[i-1][pre_sum]; f[i][sum]表示取了前i个数，异或和是sum。

因为10^6二进制是20bit，所以异或和最大是1<<20

f[i][sum]的状态由之前的两个状态转移过来，一个是i-1没取k[i]，一个是i-1取了k[i]，pre_sum^k[i]=sum，So that pre_sum=k[i]^sum

注意这里的方案数最多会有2^40个，所以f[i][sum]要用long long，另外sum=0也是一种方案（不取任何一个数），对应的f[i][sum]=1；

#include<iostream>
#include<stdio.h>
#include<cstdio>
#include<stdlib.h>
#include<vector>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<stack>
#include<queue>
#include<ctype.h>
#include<map>
#include<time.h>
#include<bitset>
using namespace std;
//hdu 5119

//1<<20=1048575
const int maxn=50;
int T;
int N;
int M;
int k[maxn];
long long f[maxn][(1<<20)+10];
int main()
{
    freopen("input.txt","r",stdin);
    //freopen("data.txt","r",stdin);
    //freopen("out1.txt","w",stdout);
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++)
    {
        scanf("%d %d",&N,&M);
        for(int i=1;i<=N;i++)
        {
            scanf("%d",&k[i]);
        }
        memset(f,0,sizeof(f));
        //f[0][0]=1;
        //选0个也是一种方案
        f[1][0]=1;
        f[1][k[1]]=1;
        for(int i=2;i<=N;i++)
        {
            for(int sum=0;sum<1048576;sum++)
            {
                f[i][sum]=f[i-1][sum]+f[i-1][k[i]^sum];
//                if(f[i-1][k[i]^sum]!=0)//排除不存在的情况
//                {
//                    f[i][sum]=max(f[i-1][sum],f[i-1][k[i]^sum]+1);
//
//                }
//                else
//                {
//                    f[i][sum]=f[i-1][sum];
//                }
            }
        }
      //  for(int i=0;i<100;i++) cout<<f[1][i]<<" ";
        long long ans=0;
        for(int i=M;i<1048576;i++)
        {
            ans+=f[N][i];
        }
        printf("Case #%d: %I64d\n",ca,ans);
    }
    return 0;
}