HDU1964 Pipes 插头DP

　　题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1964

　　空白区域之间有权值，求经过所有空白区域的哈密顿回路的最小权值。简单的插头DP，空白区域特殊处理即可。

//STATUS:C++_AC_203MS_664KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<map>
using namespace std;
#define LL long long
#define pii pair<int,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N=21,INF=0x3f3f3f3f,MOD=40001,STA=8000010;
const double DNF=1e13;

int g[N][N],code[N],ma[N];
int T,n,m,ex,ey;

struct Hash{     //Hash表,MOD为表长,STA为表大小
    int first[MOD],next[STA],size;
    int f[STA];
    LL sta[STA];
    void init(){
        size=0;
        mem(first,-1);
    }
    void add(LL st,int ans){
        int i,u=st%MOD;
        for(i=first[u];i!=-1;i=next[i]){
            if(sta[i]==st){
                f[i]=Min(f[i],ans);
                return;
            }
        }
        sta[size]=st;
        f[size]=ans;
        next[size]=first[u];
        first[u]=size++;
    }
}hs[2];

void shift(int p)    //换行移位
{
    int k;
    LL sta;
    for(k=0;k<hs[!p].size;k++){
        sta=hs[!p].sta[k]<<3;
        hs[p].add(sta,hs[!p].f[k]);
    }
}

LL getsta()   //最小表示法
{
    LL i,cnt=1,sta=0;
    mem(ma,-1);
    ma[0]=0;
    for(i=0;i<=m;i++){
        if(ma[code[i]]==-1)ma[code[i]]=cnt++;
        code[i]=ma[code[i]];
        sta|=(LL)code[i]<<(3*i);
    }
    return sta;
}

void getcode(LL sta)
{
    int i;
    for(i=0;i<=m;i++){
        code[i]=sta&7;
        sta>>=3;
    }
}

void unblock(int i,int j,int p)
{
    int k,t;
    LL cnt,x,y;
    for(k=0;k<hs[!p].size;k++){
        getcode(hs[!p].sta[k]);
        x=code[j],y=code[j+1];
        cnt=hs[!p].f[k];
        if(x && y){     //合并连通分量
            code[j]=code[j+1]=0;
            if(x!=y){
                for(t=0;t<=m;t++)
                    if(code[t]==y)code[t]=x;
                hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));
            }
            else if(i==ex && j==ey){   //最后一个点特殊处理
                hs[p].add(getsta(),cnt);
            }
        }
        else if(x&&!y || !x&&y){   //延续连通分量
            t=x?x:y;
            if(g[i+1][j]>=0){
                code[j]=t;code[j+1]=0;
                hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));
            }
            if(g[i][j+1]>=0){
                code[j]=0;code[j+1]=t;
                hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));
            }
        }
        else {  //创建新连通分量
            if(g[i][j]>=0 && g[i][j]<=9){
                hs[p].add(getsta(),cnt);
            }
            if(g[i+1][j]>=0 && g[i][j+1]>=0){
                code[j]=code[j+1]=8;
                hs[p].add(getsta(),cnt);
            }
        }
    }
}

void block(LL j,int p)
{
    int k;
    for(k=0;k<hs[!p].size;k++){
        getcode(hs[!p].sta[k]);
        if(code[j]==0 && code[j+1]==0)
            hs[p].add(getsta(),hs[!p].f[k]);
    }
}

int slove()
{
    int i,j,p;
    hs[0].init();
    hs[p=1].init();
    hs[0].add(0,0);
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            if(g[i][j]>=0)unblock(i,j,p);
            else p=!p;
            hs[p=!p].init();
        }
        shift(p);   //换行移位
        hs[p=!p].init();
    }
    for(i=0;i<hs[!p].size;i++)
        if(hs[!p].sta[i]==0)return hs[!p].f[i];
    return 0;
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j,ans;
    char s[N];
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        mem(g,-1);
        n=n*2-1;
        m=m*2-1;
        scanf("%s",s);
        getchar();
        for(i=0;i<n;i++){
            gets(s);
            for(j=1;j<=m;j++){
                if(s[j]==' ')g[i][j-1]=10;
                else if(s[j]=='#')g[i][j-1]=-1;
                else g[i][j-1]=s[j]-'0';
            }
        }
        scanf("%s",s);
        ex=n-1,ey=m-1;

        ans=slove();

        printf("%d\n",ans);
    }
    return 0;
}