hdu 1964 插头dp

题目链接

每个格子间的墙之间都有花费。要使一条回路遍历所有格子，并且总花费最小。

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
#define mes(a,x,s)  memset(a,x,(s)*sizeof a[0])
#define mem(a,x)  memset(a,x,sizeof a)
#define ysk(x)  (1<<(x))
typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
const int maxn=10    ;
const int N=300000;
const int mod=19993    ;
char mp[30][30];
int cur,n,m;
struct COST
{
   int le,up;
}cost[maxn+3][maxn+3];
struct pli
{
  ll first;int second;
};
void update(int &x,int val)  {if(val<x)  x=val; }
struct Hashmap
{
    pli e[N];
    int fir[mod+3],nex[N],nedge;
    void init()
    {
        nedge=0;
        mem(fir,-1);
    }
    void insert(ll state,int val)
    {
        int p=state%mod;
        for(int i=fir[p];~i;i=nex[i])
        {
            if(state==e[i].first)
            {
                update(e[i].second,val);
                return;
            }
        }
        e[nedge].first=state,e[nedge].second=val;
        nex[nedge]=fir[p];
        fir[p]=nedge++;
    }

}ha[2];


struct Code
{
   int tot,bit[maxn+3],ch[maxn+3];
   ll s;
   void decode(int x)
   {
       tot=0;
       for(int i=0;i<=m;i++)
       {
           bit[i]=x&7;
           tot=max(tot,bit[i]);
           x>>=3;
       }
       tot++;
    }
    ll normalize()
    {
        mem(ch,-1);ch[0]=0;tot=1;
        s=0;
        for(int i=m;i>=0;i--)
        {
            if(ch[bit[i]]==-1) ch[bit[i]]=tot++;
            bit[i]=ch[bit[i]];
            s=(s<<3)|bit[i];
        }
        return s;
    }
    void shift()
    {
        for(int i=m;i>0;i--)  bit[i]=bit[i-1];
        bit[0]=0;
    }
    void color(int a,int b)//a to b
    {
        for0(i,m+1) if(bit[i]==a)
        {
            bit[i]=b;
        }
    }
}code;


void getCost(int i,int j)
{
   int x=2*i,y=2*j;
   if( mp[x][y-1]!='#') cost[i][j].le=mp[x][y-1]-'0';
   if( mp[x-1][y]!='#') cost[i][j].up=mp[x-1][y]-'0';
}


void DP(int x,int y,ll state, int val)
{
    const int p=code.bit[y-1],q=code.bit[y];
    if(p&&q)
    {
        if(p==q)
        {
            if(x==n&&y==m)
            {
                 code.bit[y-1]=code.bit[y]=0;
                 code.shift();
                 ha[cur].insert(code.normalize(),val+cost[x][y].le+cost[x][y].up);
            }
        }
        else
        {
           code.color(p,q);
           code.bit[y-1]=code.bit[y]=0;
           if(y==m) code.shift();
           ha[cur].insert(code.normalize(),val+cost[x][y].le+cost[x][y].up);
        }
        return;
    }

    if(!p&&!q)
    {
        if(y+1<=m&&x+1<=n)
        {
            code.bit[y-1]=code.bit[y]=code.tot;
            ha[cur].insert(code.normalize(),val);
        }
        return;
    }
    int tc=p?cost[x][y].le:cost[x][y].up;
    int t=p+q;
    if(y+1<=m)
    {
        code.bit[y-1]=0,code.bit[y]=t;
        ha[cur].insert(code.normalize(),val+tc);
    }
    if(x+1<=n)
    {
        code.bit[y-1]=t,code.bit[y]=0;
        if(y==m) code.shift();
        ha[cur].insert(code.normalize(),val+tc);
    }
}

void solve()
{
    cur=0;ha[cur].init();
    ha[cur].insert(0,0);
    for1(i,n)  for1(j,m)
    {
        cur^=1;ha[cur].init();
        for0(k,mod)  for(int l=ha[1^cur].fir[k];~l;l=ha[1^cur].nex[l])
        {
            pli now=ha[cur^1].e[l];
            ll state=now.first;
            int val=now.second;
            code.decode(state);
            DP(i,j,state,val);
        }
    }
    int ans=INF;
    for0(k,mod)  for(int l=ha[cur].fir[k];~l;l=ha[cur].nex[l])
    {
        update(ans,ha[cur].e[l].second);
    }
    cout<<ans<<endl;

}
int main()
{
   std::ios::sync_with_stdio(false);
   int T;cin>>T;
   while(T--)
   {
       cin>>n>>m;cin.get();
       for1(i,2*n+1) for1(j,2*m+2) cin.get(mp[i][j]);
       for1(i,n) for1(j,m) {mem(&cost[i][j],0);getCost(i,j);}
       solve();
   }
   return 0;
}