【POJ】3133 Manhattan Wiring

http://poj.org/problem?id=3133

题意：n×m的网格，有2个2，2个3，他们不会重合。还有障碍1。现在求2到2的路径和3到3的路径互不相交的最短长度-2。（2<=n,m<=9）

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

struct H {
	static const int M=1000007;
	struct E { int next, to; }e[M<<1];
	int mn[M], hash[M], cnt, ihead;
	H() { cnt=ihead=0; memset(mn, 0x3f, sizeof mn); memset(hash, -1, sizeof hash); }
	bool find(int x, int &pos) {
		pos=x%M;
		while(1) { if(hash[pos]==x) return 0; if(hash[pos]==-1) break; ++pos; if(pos==M) pos=0; }
		hash[pos]=x;
		return 1;
	}
	void ins(int a, int b) {
		int pos; if(!find(a, pos)) { mn[pos]=min(mn[pos], b); return; }
		mn[pos]=b; e[++cnt].next=ihead; ihead=cnt; e[cnt].to=pos;
	}
	void clr() { for(int i=ihead; i; i=e[i].next) hash[e[i].to]=-1, mn[e[i].to]=0x3f3f3f3f; ihead=cnt=0; }
}T[2];
int n, m, mp[12][12], all;
#define BIT(a,b) ((a)<<((b)<<1))
#define CLR(a,b) (a^=(a&BIT(3,b)))
#define GET(a,b) (3&((a)>>((b)<<1)))

bool next(int s, int row, int col, bool U, bool D, bool L, bool R, int id, int &t) {
	if((!row&&U) || (!col&&L) || (row==n-1&&D) || (col==m-1&&R)) return 0;
	if((D&&mp[row+1][col]==1) || (R&&mp[row][col+1]==1)) return 0;
	int l=GET(s, col), u=GET(s, col+1), d=0, r=0;
	if((L&&!l) || (!L&&l) || (U&&!u) || (!U&&u)) return 0;
	if((l&&l!=id) || (u&&u!=id)) return 0;
	t=s;
	if(!l && !u) { if(D) d=id; if(R) r=id; }
	else if(l && u) { if(l!=u || mp[row][col]>=1) return 0; }
	else if(l && !u) { if(D) d=l; if(R) r=l; }
	else if(!l && u) { if(D) d=u; if(R) r=u; }
	CLR(t, col); CLR(t, col+1);
	t|=BIT(d, col); t|=BIT(r, col+1);
	if(col==m-1) t<<=2; t&=all;
	return 1;
}

int bfs() {
	int now=0, last=1;
	T[now].clr();
	T[now].ins(0, 0);
	for(int row=0; row<n; ++row) for(int col=0; col<m; ++col) {
		T[last].clr();
		for(int i=T[now].ihead; i; i=T[now].e[i].next) {
			int t, s=T[now].hash[T[now].e[i].to], d=T[now].mn[T[now].e[i].to]+1, id=mp[row][col]-1;
			if(id>=1) {
				if(next(s, row, col, 1, 0, 0, 0, id, t)) T[last].ins(t, d);
				if(next(s, row, col, 0, 1, 0, 0, id, t)) T[last].ins(t, d);
				if(next(s, row, col, 0, 0, 1, 0, id, t)) T[last].ins(t, d);
				if(next(s, row, col, 0, 0, 0, 1, id, t)) T[last].ins(t, d);
			}
			else {
				if(next(s, row, col, 0, 0, 0, 0, 0, t)) T[last].ins(t, d-1);
				if(id==-1) {
					for(int k=1; k<=2; ++k) {
						if(next(s, row, col, 1, 1, 0, 0, k, t)) T[last].ins(t, d);
						if(next(s, row, col, 1, 0, 1, 0, k, t)) T[last].ins(t, d);
						if(next(s, row, col, 1, 0, 0, 1, k, t)) T[last].ins(t, d);
						if(next(s, row, col, 0, 1, 1, 0, k, t)) T[last].ins(t, d);
						if(next(s, row, col, 0, 1, 0, 1, k, t)) T[last].ins(t, d);
						if(next(s, row, col, 0, 0, 1, 1, k, t)) T[last].ins(t, d);
					}
				}
			}
		}
		swap(now, last);
	}
	int ret=0x3f3f3f3f;
	for(int i=T[now].ihead; i; i=T[now].e[i].next) ret=min(ret, T[now].mn[T[now].e[i].to]);
	return ret>=0x3f3f3f3f?0:ret-2;
}

int main() {
	while(scanf("%d%d", &n, &m), n|m) {
		for(int i=0; i<n; ++i) for(int j=0; j<m; ++j) scanf("%d", &mp[i][j]);
		all=BIT(1, m+1)-1;
		printf("%d
", bfs());
	}
	return 0;
}

　　

---

吐槽：

1、改着改着自己本来特判的东西然后改掉了QAQ，然后造数据对拍才发现....

2、为什么我写的插头dp常数那么大啊！！！

本题比简单回路简单多了吧？

自己很容易可以画出左插头和上插头的情况，然后在位中维护的是当前插头所连接的标号然后就行啦...