POJ-3133 Manhattan Wiring 插头DP

　　题目链接：http://poj.org/problem?id=3133

　　用两条线段分别连接2-2和3-3，使得他们长度的总和最小。

　　典型的插头DP题目，用0,2,3,种状态来表示每个格子的插头分别是没有插头，2号插头，3号插头，然后把所有可能的转移情况列出来，只是情况比较多，而且一定要考虑仔细，不然检查代码会很惨的。。话说A完这道题后感觉真爽啊O(∩_∩)O~

  1 //STATUS:C++_AC_360MS_900KB

  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=10,INF=0x3f3f3f3f,MOD=60007,STA=1000010;

 22 const LL LNF=0x3f3f3f3f3f3f3f3f;

 23 const double DNF=100000000000;

 24 

 25 int g[N][N];

 26 int n,m;

 27 

 28 struct Hash{

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

 30     int f[STA],sta[STA];

 31     void init(){

 32         size=0;

 33         mem(first,-1);

 34     }

 35     void add(int 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         f[size]=ans;

 44         sta[size]=st;

 45         next[size]=first[u];

 46         first[u]=size++;

 47     }

 48 }hs[2];

 49 

 50 void shift(int p)

 51 {

 52     int i,flag=(1<<((m+1)<<1))-1;

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

 54         hs[p].add((hs[!p].sta[i]<<2)&flag,hs[!p].f[i]);

 55     }

 56 }

 57 

 58 int slove()

 59 {

 60     int i,j,k,x,y,left,up,p,sta,f;

 61     hs[0].init();

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

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

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

 65         for(j=0;j<m;j++,hs[p=!p].init()){

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

 67                 sta=hs[!p].sta[k];

 68                 f=hs[!p].f[k]+1;

 69                 x=(3<<(j<<1))&sta;

 70                 y=(3<<((j+1)<<1))&sta;

 71                 left=x>>(j<<1);

 72                 up=y>>((j+1)<<1);

 73                 if(g[i][j]){

 74                     if(g[i][j]==2){  //独立插头

 75                         if(left==2 && up==0){

 76                             hs[p].add(sta^x,f);

 77                         }

 78                         else if(left==0 && up==2){

 79                             hs[p].add(sta^y,f);

 80                         }

 81                         else if(left==0 && up==0){

 82                             if(g[i+1][j]==1 || g[i+1][j]==2)hs[p].add(sta|(2<<(j<<1)),f);

 83                             if(g[i][j+1]==1 || g[i][j+1]==2)hs[p].add(sta|(2<<((j+1)<<1)),f);

 84                         }

 85                     }

 86                     else if(g[i][j]==3){  //独立插头

 87                         if(left==3 && up==0){

 88                             hs[p].add(sta^x,f);

 89                         }

 90                         else if(left==0 && up==3){

 91                             hs[p].add(sta^y,f);

 92                         }

 93                         else if(left==0 && up==0){

 94                             if(g[i+1][j]==1 || g[i+1][j]==3)hs[p].add(sta|(3<<(j<<1)),f);

 95                             if(g[i][j+1]==1 || g[i][j+1]==3)hs[p].add(sta|(3<<((j+1)<<1)),f);

 96                         }

 97                     }

 98                     else {  //g[i][j]==1

 99                         if(left==up && left){

100                             hs[p].add(sta^x^y,f);

101                         }

102                         else if((left==2 && up==0) || (left==0 && up==2)){

103                             if(g[i+1][j])hs[p].add(sta^y|(2<<(j<<1)),f);

104                             if(g[i][j+1])hs[p].add(sta^x|(2<<((j+1)<<1)),f);

105                         }

106                         else if((left==3 && up==0) || (left==0 && up==3)){

107                             if(g[i+1][j])hs[p].add(sta^y|(3<<(j<<1)),f);

108                             if(g[i][j+1])hs[p].add(sta^x|(3<<((j+1)<<1)),f);

109                         }

110                         else if(left==0 && up==0){   //这里一定要考虑仔细

111                             if(g[i+1][j] && g[i][j+1]){

112                                 hs[p].add(sta|(2<<(j<<1))|(2<<((j+1)<<1)),f);

113                                 hs[p].add(sta|(3<<(j<<1))|(3<<((j+1)<<1)),f);

114                             }

115                             hs[p].add(sta,f-1);

116                         }

117                     }

118                 }

119                 else if(left==0 && up==0){

120                     hs[p].add(sta,f-1);

121                 }

122             }

123         }

124         shift(p);

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

126     }

127     p=!p;

128     for(i=0;i<hs[p].size;i++){

129         if(hs[p].sta[i]==0)return hs[p].f[i]-2;

130     }

131     return 0;

132 }

133 

134 int main()

135 {

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

137     int i,j,ans;

138     while(~scanf("%d%d",&n,&m) && (n || m))

139     {

140         mem(g,0);

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

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

143                 scanf("%d",&g[i][j]);

144                 if(g[i][j]==0)g[i][j]=1;

145                 else if(g[i][j]==1)g[i][j]=0;

146             }

147         }

148 

149         ans=slove();

150 

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

152     }

153     return 0;

154 }