HDU-1964 Pipes 插头DP

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

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

  1 //STATUS:C++_AC_203MS_664KB

  2 #include<stdio.h>

  3 #include<stdlib.h>

  4 #include<string.h>

  5 #include<math.h>

  6 #include<iostream>

  7 #include<string>

  8 #include<algorithm>

  9 #include<vector>

 10 #include<queue>

 11 #include<stack>

 12 #include<map>

 13 using namespace std;

 14 #define LL long long

 15 #define pii pair<int,int>

 16 #define Max(a,b) ((a)>(b)?(a):(b))

 17 #define Min(a,b) ((a)<(b)?(a):(b))

 18 #define mem(a,b) memset(a,b,sizeof(a))

 19 #define lson l,mid,rt<<1

 20 #define rson mid+1,r,rt<<1|1

 21 const int N=21,INF=0x3f3f3f3f,MOD=40001,STA=8000010;

 22 const double DNF=1e13;

 23 

 24 int g[N][N],code[N],ma[N];

 25 int T,n,m,ex,ey;

 26 

 27 struct Hash{     //Hash表,MOD为表长,STA为表大小

 28     int first[MOD],next[STA],size;

 29     int f[STA];

 30     LL sta[STA];

 31     void init(){

 32         size=0;

 33         mem(first,-1);

 34     }

 35     void add(LL st,int ans){

 36         int i,u=st%MOD;

 37         for(i=first[u];i!=-1;i=next[i]){

 38             if(sta[i]==st){

 39                 f[i]=Min(f[i],ans);

 40                 return;

 41             }

 42         }

 43         sta[size]=st;

 44         f[size]=ans;

 45         next[size]=first[u];

 46         first[u]=size++;

 47     }

 48 }hs[2];

 49 

 50 void shift(int p)    //换行移位

 51 {

 52     int k;

 53     LL sta;

 54     for(k=0;k<hs[!p].size;k++){

 55         sta=hs[!p].sta[k]<<3;

 56         hs[p].add(sta,hs[!p].f[k]);

 57     }

 58 }

 59 

 60 LL getsta()   //最小表示法

 61 {

 62     LL i,cnt=1,sta=0;

 63     mem(ma,-1);

 64     ma[0]=0;

 65     for(i=0;i<=m;i++){

 66         if(ma[code[i]]==-1)ma[code[i]]=cnt++;

 67         code[i]=ma[code[i]];

 68         sta|=(LL)code[i]<<(3*i);

 69     }

 70     return sta;

 71 }

 72 

 73 void getcode(LL sta)

 74 {

 75     int i;

 76     for(i=0;i<=m;i++){

 77         code[i]=sta&7;

 78         sta>>=3;

 79     }

 80 }

 81 

 82 void unblock(int i,int j,int p)

 83 {

 84     int k,t;

 85     LL cnt,x,y;

 86     for(k=0;k<hs[!p].size;k++){

 87         getcode(hs[!p].sta[k]);

 88         x=code[j],y=code[j+1];

 89         cnt=hs[!p].f[k];

 90         if(x && y){     //合并连通分量

 91             code[j]=code[j+1]=0;

 92             if(x!=y){

 93                 for(t=0;t<=m;t++)

 94                     if(code[t]==y)code[t]=x;

 95                 hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));

 96             }

 97             else if(i==ex && j==ey){   //最后一个点特殊处理

 98                 hs[p].add(getsta(),cnt);

 99             }

100         }

101         else if(x&&!y || !x&&y){   //延续连通分量

102             t=x?x:y;

103             if(g[i+1][j]>=0){

104                 code[j]=t;code[j+1]=0;

105                 hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));

106             }

107             if(g[i][j+1]>=0){

108                 code[j]=0;code[j+1]=t;

109                 hs[p].add(getsta(),cnt+(g[i][j]<10?g[i][j]:0));

110             }

111         }

112         else {  //创建新连通分量

113             if(g[i][j]>=0 && g[i][j]<=9){

114                 hs[p].add(getsta(),cnt);

115             }

116             if(g[i+1][j]>=0 && g[i][j+1]>=0){

117                 code[j]=code[j+1]=8;

118                 hs[p].add(getsta(),cnt);

119             }

120         }

121     }

122 }

123 

124 void block(LL j,int p)

125 {

126     int k;

127     for(k=0;k<hs[!p].size;k++){

128         getcode(hs[!p].sta[k]);

129         if(code[j]==0 && code[j+1]==0)

130             hs[p].add(getsta(),hs[!p].f[k]);

131     }

132 }

133 

134 int slove()

135 {

136     int i,j,p;

137     hs[0].init();

138     hs[p=1].init();

139     hs[0].add(0,0);

140     for(i=0;i<n;i++){

141         for(j=0;j<m;j++){

142             if(g[i][j]>=0)unblock(i,j,p);

143             else p=!p;

144             hs[p=!p].init();

145         }

146         shift(p);   //换行移位

147         hs[p=!p].init();

148     }

149     for(i=0;i<hs[!p].size;i++)

150         if(hs[!p].sta[i]==0)return hs[!p].f[i];

151     return 0;

152 }

153 

154 int main()

155 {

156  //   freopen("in.txt","r",stdin);

157     int i,j,ans;

158     char s[N];

159     scanf("%d",&T);

160     while(T--)

161     {

162         scanf("%d%d",&n,&m);

163         mem(g,-1);

164         n=n*2-1;

165         m=m*2-1;

166         scanf("%s",s);

167         getchar();

168         for(i=0;i<n;i++){

169             gets(s);

170             for(j=1;j<=m;j++){

171                 if(s[j]==' ')g[i][j-1]=10;

172                 else if(s[j]=='#')g[i][j-1]=-1;

173                 else g[i][j-1]=s[j]-'0';

174             }

175         }

176         scanf("%s",s);

177         ex=n-1,ey=m-1;

178 

179         ans=slove();

180 

181         printf("%d\n",ans);

182     }

183     return 0;

184 }