hdu 1964 Plan (插头dp)

题意：

给出一个迷宫，并且迷宫可走的格子量量之间有间隔墙，打破强需要一些话费，现在问如何走完所有的格子并且让总的花费最小。

题解：

这题和之前的计数类问题相似只不过改动dp方程一些地方而已，其实这个和之前的状态压缩题目对比就很容易理解这样转移。轻松1a，kuangbin的模板就是牛逼！赞！

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<map>
using namespace std;
typedef __int64 lld;
#define oo 0x3f3f3f3f
#define OO 0x3f3f3f3f3f3f3f3f
#define HASH 10007
#define STATE 1000010
#define MAXD 15
int N,M,ex,ey;
int code[MAXD],maze[MAXD][MAXD];
int ch[MAXD];
struct NODE
{
    int right,down;
}node[MAXD][MAXD];

struct HASHMAP
{
    int head[HASH],next[STATE],sizes;
    int dp[STATE];
    lld state[STATE];
    void init()
    {
        sizes=0;
        memset(head,-1,sizeof head);
    }
    void push(lld st,int ans)
    {
        int h=st%HASH;
        for(int i=head[h];i!=-1;i=next[i])
        {
            if(st==state[i])
            {
                dp[i]=min(dp[i],ans);
                return ;
            }
        }
        dp[sizes]=ans;
        state[sizes]=st;
        next[sizes]=head[h];
        head[h]=sizes++;
    }
}hm[2];

void decode(int code[],int m,lld st)
{
    for(int i=m;i>=0;i--)
    {
        code[i]=st&7;
        st>>=3;
    }
}

lld encode(int code[],int m)///最小表示法
{
    lld st=0;
    int cnt=0;
    memset(ch,-1,sizeof ch);
    ch[0]=0;
    for(int i=0;i<=m;i++)
    {
        if(ch[code[i]]==-1) ch[code[i]]=++cnt;
        code[i]=ch[code[i]];
        st<<=3;
        st|=code[i];
    }
    return st;
}

void shift(int code[],int m)///换行 移位
{
    for(int i=m;i>0;i--)
        code[i]=code[i-1];
    code[0]=0;
}

void dpblank(int i,int j,int cur)
{
    int left,up;
    for(int k=0;k<hm[cur].sizes;k++)
    {
        decode(code,M,hm[cur].state[k]);
        left=code[j-1];
        up=code[j];
        if(left&&up)///11 -> 00       有上插头和左插头，这种情况下相当于合并两个连通分量
        {
            if(left==up)///如果在同一个联通分量
            {
                if(i==ex&&j==ey)
                {
                    code[j-1]=code[j]=0;
                    if(j==M)shift(code,M);
                    hm[cur^1].push(encode(code,M),hm[cur].dp[k]);
                }
            }
            else///不再同一个联通分量里面可以进行合并
            {
                code[j-1]=code[j]=0;
                for(int t=0;t<=M;t++)
                    if(code[t]==up)
                        code[t]=left;
                if(j==M)shift(code,M);
                hm[cur^1].push(encode(code,M),hm[cur].dp[k]);
            }
        }
        else if(left||up)///01 || 10  上插头和左插头恰好有一个，这种情况相当于延续原来的连通分量
        {
            int temp;
            if(left) temp=left;
            else temp=up;
            if(maze[i][j+1])
            {
                code[j-1]=0;
                code[j]=temp;
                hm[cur^1].push(encode(code,M),hm[cur].dp[k]+node[i][j].right);
            }
            if(maze[i+1][j])
            {
                code[j-1]=temp;
                code[j]=0;
                if(j==M)shift(code,M);///切记不可忘记，换行要shift
                hm[cur^1].push(encode(code,M),hm[cur].dp[k]+node[i][j].down);
            }
        }
        else///没有上插头和左插头，有下插头和右插头，相当于构成一个新的连通块
        {
            if(maze[i][j+1]&&maze[i+1][j])
            {
                code[j]=code[j-1]=13;
                hm[cur^1].push(encode(code,M),hm[cur].dp[k]+node[i][j].down+node[i][j].right);
            }
        }
    }
}

void init()
{
    memset(maze,0,sizeof maze);
    char c,str[2*MAXD+1][2*MAXD+1];
    scanf("%d %d",&N,&M);
    getchar();
    for(int i=1;i<=2*N+1;i++)
        gets(str[i]+1);
    for(int i=1;i<=N;i++)
    {
        for(int j=1;j<=M;j++)
        {
            maze[i][j]=1;
            node[i][j].right=str[2*i][2*j+1]-'0';
            node[i][j].down=str[2*i+1][2*j]-'0';
        }
    }
    ex=N;
    ey=M;
}

void solve()
{
    int cur=0;
    lld ans=0;
    hm[cur].init();
    hm[cur].push(0,0);
    for(int i=1;i<=N;i++)
        for(int j=1;j<=M;j++)
        {
            hm[cur^1].init();
            dpblank(i,j,cur);
            cur^=1;
        }
    for(int i=0;i<hm[cur].sizes;i++)
            ans+=hm[cur].dp[i];
    printf("%I64d\n",ans);
}

void Debug()
{
    for(int i=1;i<=N;i++)
    {
        for(int j=1;j<=M;j++)
            printf("%d %d: right=%d down=%d\n",i,j,node[i][j].right,node[i][j].down);
    }
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        init();
        //Debug();
        solve();
    }
    return 0;
}
/**
2
2 2
OO
O*
4 4
***O
XO**
**O*
XX**

*/