hdu 4317 Unfair Nim (状态压缩DP) 【2012 MultiUniversity Training Contest 2】

一道状态压缩DP，比赛时没时间仔细想了，赛后想了想转移那部分也有点糊里糊涂的，后来请教了一下yk同学，才肯定了一下那个状态表示，后来自己一边探索，一边敲码，一边调试，确认无误后得到1A。

看到题很自然的想到从低位开始往高位处理，基本处理加的思路是：在本位上只有可能+1（由1变成0，由0变成1）或者+2（由0变成0，由1变成1，但是存在进位），基本的情况是本位只可能主动+1，如果+2则可能是后一位进位过来了的。

基本状态是dp[i][j]：处理玩了第i位（第0～i位异或已经全为0），第i位有j堆需要想第i+1位进位。

状态转移：由第i+1位本身的值，加上第i位至第i+1位的进位，会得到一串新的值，但是如果不是偶数个1的话，第i+1位的异或结果是不会为0的，并且第i+1位也是可以要有自己本身的状态的，所以要再枚举2^n种变化，只对能够得到异或结果为0的状态进行更新。

更新的方法：

　　最后第i+1位的结果是由第i+1位本身的值，第i位进位的值，以及枚举的第i+1位是否要变化+1的值，三个值一起决定的，第i+1位本身的值 异或 第i位进位的值，得到的是目前第i+1位的值，可以先预处理出目前的情况与那些变化相结合是可以得到合法状态的，继而跟新。将3个值相异或即可得到第i+1位剩下的结果，是偶数个1就是合法的结果。

　　而也要知道第i+1位向第i+2位进位的结果，仔细观察后可以发现如果那三个值中某一位存在两个或两个以上的1时就会有对上一位的进位，而怎么顺利得到呢？这个就要仰慕一下severus大牛的指导，居然用到了《数电》里面的知识，F=ABC+AB¬C+A¬BC+¬ABC，化简之后就是得到的AB+(A^B)C，这样dp[i][]就可以由dp[i-1][]转化得到了。下面附代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 10
#define inf (1<<30)
#define see(x) cout<<#x<<":"<<x<<endl;
using namespace std;
const int maxn = 1<<N;
bool cha[maxn][maxn];
int dp[20][maxn];
int num[20];
int c[maxn];
inline int cal(int x){        //取出X里2进制1的个数
    if(c[x]!=-1) return c[x];
    int y = x;
    x=(x & 0x55555555) + ((x >>1 ) & 0x55555555);
    x=(x & 0x33333333) + ((x >>2 ) & 0x33333333);
    x=(x & 0x0F0F0F0F) + ((x >>4 ) & 0x0F0F0F0F);
    x=(x & 0x00FF00FF) + ((x >>8 ) & 0x00FF00FF);
    x=(x & 0x0000FFFF) + ((x >>16) & 0x0000FFFF);
    return c[y]=x;
}
int a[N];
int main(){
    int n, i, j, k, m, ans;
    memset(cha,0,sizeof(cha));
    memset(c,-1,sizeof(c));
    for(i=0;i<maxn;i++){
        for(j=0;j<maxn;j++){
            k = i^j;
            if((cal(k)&1)==0){
                cha[i][j] = 1;
            }
        }
    }
    while(~scanf("%d",&n)){
        m = 0;
        memset(num,0,sizeof(num));
        for(i=0;i<n;i++){
            scanf("%d",&a[i]);
            for(k=0;a[i]>=(1<<k);k++){
                if(a[i]&(1<<k)){
                    num[k] += (1<<i);
                }
            }
            m = max(m,k);
        }
        for(i=0;i<=m;i++) for(j=0;j<(1<<n);j++){
            dp[i][j] = inf;
        }
        for(k=0;k<(1<<n);k++){
            if(cha[num[0]][k]){
                dp[0][num[0]&k] = min(dp[0][num[0]&k],cal(k));
            }
        }
        for(i=1;i<=m;i++){
            for(j=0;j<(1<<n);j++){
                for(k=0;k<(1<<n);k++){
                    if(cha[num[i]^j][k]){
                        dp[i][(num[i]&j)|((num[i]^j)&k)] = min(dp[i][(num[i]&j)|((num[i]^j)&k)],dp[i-1][j] + cal(k)*(1<<i));
                    }
                }
            }
        }
        ans = dp[m][0];
        printf("%d\n",ans);
    }
    return 0;
}