poj 1739 Tony's Tour 插头dp模板题

题意：

给一个迷宫，求左下角到右下角的路径数。

分析：

插头dp的模板题，建议先看cdq的论文再看代码，这份代码在模板基础上略微有改动。论文地址http://wenku.baidu.com/view/ed2b3e23482fb4daa58d4b74.html

代码：

#include <iostream>
using namespace std;
const int maxD=16;
const int HASH=10007;
const int STATE=1000024;

int N,M;
int maze[maxD][maxD];
int code[maxD];
int ch[maxD];

struct HASHMAP
{
	int head[HASH],next[STATE],size;
	long long state[STATE],f[STATE];
	void init(){
		size=0;
		f[0]=0;
		memset(head,-1,sizeof(head));
	}
	void push(long long st,long long ans){
		int h=st%HASH;
		for(int i=head[h];i!=-1;i=next[i])
			if(state[i]==st){
				f[i]+=ans;
				return ;
			}
		state[size]=st;
		f[size]=ans;
		next[size]=head[h];
		head[h]=size++;
	}		
}hm[2];

void decode(int *code,int m,long long st)
{
	for(int i=m;i>=0;--i){
		code[i]=st&7;
		st>>=3;
	}	
}

long long encode(int *code,int m)
{
	int cnt=1;
	memset(ch,-1,sizeof(ch));
	long long st=0;	
	ch[0]=0;
	for(int i=0;i<=m;++i){
		if(ch[code[i]]==-1) ch[code[i]]=cnt++;
		code[i]=ch[code[i]];
		st<<=3;
		st|=code[i];
	}
	return st;
}

void shift(int *code,int m)
{
	for(int i=m;i>0;--i) code[i]=code[i-1];
	code[0]=0;	
}

void init()
{
	char str[20];
	memset(maze,0,sizeof(maze));
	for(int i=1;i<=N;++i){
		scanf("%s",str);
		for(int j=0;j<M;++j)
			if(str[j]=='.')
				maze[i][j+1]=1;
	}	
}

void dpblank(int i,int j,int cur)
{
	int k,left,up;
	for(k=0;k<hm[cur].size;++k){
		decode(code,M,hm[cur].state[k]);
		left=code[j-1];
		up=code[j];
		
		if((i==N&&j==1)||(i==N&&j==M)){
			if((left&&!up)||(!left&&up)){
				code[j]=code[j-1]=0;
				if(j==M) shift(code,M);
				hm[cur^1].push(encode(code,M),hm[cur].f[k]);
			}else if(left==0&&up==0){
				if(maze[i][j+1]){
					code[j-1]=0;
					code[j]=13;
					hm[cur^1].push(encode(code,M),hm[cur].f[k]);
				}
			}
			continue;		
		}
		if(left&&up){
			if(left!=up){
				code[j]=code[j-1]=0;
				for(int t=0;t<=M;++t)
					if(code[t]==left)
						code[t]=up;
				if(j==M) shift(code,M);
				hm[cur^1].push(encode(code,M),hm[cur].f[k]);
			}
		}
		else if((left&&!up)||(!left&&up)){
			int t;
			if(left) t=left;
			else t=up;
			if(maze[i][j+1]){
				code[j-1]=0;
				code[j]=t;
				hm[cur^1].push(encode(code,M),hm[cur].f[k]);
			}
			if(maze[i+1][j]){
				code[j-1]=t;
				code[j]=0;
				if(j==M) shift(code,M);
				hm[cur^1].push(encode(code,M),hm[cur].f[k]);
			}
		}
		else{
			if(maze[i][j+1]&&maze[i+1][j]){
				code[j-1]=code[j]=13;
				hm[cur^1].push(encode(code,M),hm[cur].f[k]);
			}
		}
	}	
}

void dpblock(int i,int j,int cur)
{
	for(int k=0;k<hm[cur].size;++k){
		decode(code,M,hm[cur].state[k]);
		code[j-1]=code[j]=0;
		if(j==M) shift(code,M);
		hm[cur^1].push(encode(code,M),hm[cur].f[k]);
	}		
}

void solve()
{
	int i,j,cur=0;
	hm[cur].init();
	hm[cur].push(0,1);
	for(int i=1;i<=N;++i)
		for(int j=1;j<=M;++j){
			hm[cur^1].init();
			if(maze[i][j]) dpblank(i,j,cur);
			else dpblock(i,j,cur);
			cur^=1;
		}
	printf("%I64d\n",hm[cur].f[0]);
}

int main()
{
	while(scanf("%d%d",&N,&M)==2){
		if(N==0&&M==0) break;
		init();
		solve();
	}
	return 0;	
}