[Poj3133]Manhattan Wiring （插头DP）

Description

题目大意：给你个N x M（1≤N, M≤9）的矩阵，0表示空地，1表示墙壁，2和3表示两对关键点。现在要求在两对关键点之间建立两条路径，其中两条路径不可相交或者自交（就是重复经过同一格子），并且不能经过墙壁，路径只能从一个格子走到相邻的下一格子。求两条路径最少需要经过的格子数减二。如果不存在解，输出0。

Solution

插头DP，用三进制来表示轮廓线，0表示没有插头，1表示“2”的的插头，2表示“3”的插头

1. 对于0的格子：00 -> 00，11，22；01/10 -> 01/10；02/20 -> 02/20；11/22 -> 00。
2. 对于1的格子：00 -> 00。
3. 对于2的格子：01/10 -> 00，00 -> 01/10。
4. 对于3的格子：02/20 -> 00，00 -> 02/20。

但是if太多会超时，所以有必要减少if语句的使用

Code

#include <cstdio>

#include <cstring>

#include <algorithm>

#define N 11

using namespace std;

const int A[]={1,3,9,27,81,243,729,2187,6561,19683,59049};

int n,m,g[N][N],dp[N][N][59049];

inline int read(){

	int x=0,f=1;char ch=getchar();

	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}

	return x*f;

}

int main(){

	while(~scanf("%d%d",&n,&m)&&n+m){

		memset(g,0,sizeof(g));

		memset(dp,0x3f,sizeof(dp));

		for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)g[i][j]=read();

		dp[0][m][0]=0;

		for(int i=1;i<=n;++i){

			for(int j=0;j<A[m];++j)//A[m]

				dp[i][0][j*3]=dp[i-1][m][j];

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

				for(int k=0;k<A[m+1];++k){//A[m+1]

					if (j==m&&k>=A[m]) break;

					int a=(k/A[j-1])%3,b=(k/A[j])%3;

					if(g[i][j]==1){

						if(a==0&&b==0) dp[i][j][k]=dp[i][j-1][k];

					}else if(g[i][j]==2){//min,+1

						if(a==1&&b==0) dp[i][j][k]=dp[i][j-1][k-A[j-1]]+1;

						else if(a==0&&b==1) dp[i][j][k]=dp[i][j-1][k-A[j]]+1;

						else if(a==0&&b==0) dp[i][j][k]=min(dp[i][j-1][k+A[j-1]],dp[i][j-1][k+A[j]])+1;

					}else if(g[i][j]==3){

						if(a==2&&b==0) dp[i][j][k]=dp[i][j-1][k-2*A[j-1]]+1;

						else if(a==0&&b==2) dp[i][j][k]=dp[i][j-1][k-2*A[j]]+1;

						else if(a==0&&b==0) dp[i][j][k]=min(dp[i][j-1][k+2*A[j-1]],dp[i][j-1][k+2*A[j]])+1;

					}

					else if(g[i][j]==0){

						if(a==0&&b==0) dp[i][j][k]=min(dp[i][j-1][k],min(dp[i][j-1][k+A[j]+A[j-1]]+1,dp[i][j-1][k+2*A[j]+2*A[j-1]]+1));

						else if(a==0&&b>0) dp[i][j][k]=min(dp[i][j-1][k],dp[i][j-1][k+b*A[j-1]-b*A[j]])+1;

						else if(a>0&&b==0) dp[i][j][k]=min(dp[i][j-1][k],dp[i][j-1][k+a*A[j]-a*A[j-1]])+1;

						else if(a&&b&&a==b) dp[i][j][k]=dp[i][j-1][k-a*A[j]-a*A[j-1]]+1;

					}

					if (i==n&&j==m) break;

				}

			}

		}

		if(dp[n][m][0]<0x3f3f3f3f) printf("%d\n",dp[n][m][0]-2);

		else printf("0\n");

	}

	return 0;

}