POJ 2411 插头DP

//插头DP，算是广义路径的吧。
/*
我是这样想的，定义填数的为0，未填的为1.然后，初始自然是(0,0).我还定义了整个棋盘的状态，不知是否多此一举。
这样，把轮廓线上的格子状态记录。当（I,J）上方的格子为空，必定要填一个竖的。当左边格子为空，当前可填一个横的，也可不填。
当左边格子不为空，当前格子必为空。。。AC.
*/

#include <iostream>
#include <cstdio>
using namespace std;
int n,m;
int code[15];
const int MAXL=6000;
const int MAXM=1000000;
struct HASHMAP{
    int hash[MAXL],next[MAXM];
    __int64 f[MAXM]; int state[MAXM],alst[MAXM];
    int size;
    void init(){
        size=0;
        memset(hash,-1,sizeof(hash));
    }
    void push(int st,__int64 ans,int als){
        int h=st%MAXL;
        for(int i=hash[h];i!=-1;i=next[i]){
            if(state[i]==st){
                f[i]+=ans;
                return ;
            }
        }
        f[size]=ans;
        state[size]=st;
        alst[size]=als;
        next[size]=hash[h];
        hash[h]=size++;
    }
}hm[2];

void decode(int st){
    for(int i=1;i<=m;i++){
        code[i]=(st&1);
        st=(st>>1);
    }
}

int encode(){
    int st=0; int i;
    for(i=m;i>1;i--){
        st=(st|code[i]);
        st=(st<<1);
    }
    st=(st|code[i]);
    return st;
}

void dp(int i,int j,int cur){
    for(int e=0;e<hm[cur].size;e++){
        decode(hm[cur].state[e]);
        int als=hm[cur].alst[e];
        if(code[j]==1){
            code[j]=0; als--;
            hm[cur^1].push(encode(),hm[cur].f[e],als);
        }
        else{
            if(code[j-1]==1){
                code[j]=1;
                hm[cur^1].push(encode(),hm[cur].f[e],als+1);
                code[j]=code[j-1]=0;
                hm[cur^1].push(encode(),hm[cur].f[e],als-1);
            }
            else{
                code[j]=1;
                hm[cur^1].push(encode(),hm[cur].f[e],als+1);
            }
        }
    }
}


int main(){
    int i,j;
    while(scanf("%d%d",&n,&m)!=EOF){
        if(!n&&!m) break;
        code[0]=0;
        int cur=0;
        hm[cur].init();
        hm[cur].push(0,1,0);
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++){
                hm[cur^1].init();
                dp(i,j,cur);
                cur=cur^1;
            }
        __int64 ans=0;
        for(i=0;i<hm[cur].size;i++){
            if(hm[cur].alst[i]==0){
                ans+=hm[cur].f[i];
            }
        }
        printf("%I64d
",ans);
    }
    return 0;
}