poj1739 Tony‘s Tour

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

普通通路问题，左下角走到右下角，那么只要在上面加两层壳，一层全不可走，一层全可走，那么总左下角走到右下角是一条路径且他和从外面这层壳的唯一路径走形成了一条回路，问题就转换成了直接求回路，就等同于模板题了

https://blog.csdn.net/liufengwei1/article/details/108114948

#include<cstdio>
using namespace std;
typedef long long ll;

const int maxl=15;
const int hs=299997;

int n,m,now,mt,cnt[2];ll ans;
int mi[maxl],mark[hs+3];
int q[2][hs+3],hd[hs+3],nxt[hs+3];
ll dp[2][hs+3];
int a[maxl][maxl];
char s[maxl];

inline void prework()
{
	for(int i=3;i<=n+2;i++)
	{
		scanf("%s",s+1);
		for(int j=1;j<=m;j++)
			a[i][j+2]=(s[j]=='.');
	}
	n+=2;m+=4;
	for(int j=1;j<=m;j++)
		a[1][j]=1;
	a[n][1]=a[n][2]=a[n][m]=a[n][m-1]=1;
	for(int i=1;i<=n;i++)
		a[i][1]=a[i][m]=1;
	for(int i=2;i<=n-1;i++)
		a[i][2]=a[i][m-1]=0;
	for(int j=2;j<=m-1;j++)
		a[2][j]=0;
	for(int i=0;i<=n;i++)
		a[i][0]=a[i][m+1]=0;
	for(int i=0;i<=m;i++)
		a[0][i]=a[n+1][i]=0;
}

inline void insert(int bit,ll val)
{
	int u=bit%hs+1;
	if(mark[u]!=mt)
		hd[u]=0,mark[u]=mt;
	for(int i=hd[u];i;i=nxt[i])
	if(q[now][i]==bit)
	{
		dp[now][i]+=val;
		return;
	}
	nxt[++cnt[now]]=hd[u];
	hd[u]=cnt[now];
	q[now][cnt[now]]=bit;
	dp[now][cnt[now]]=val;
}

inline void mainwork()
{
	int bit,b1,b2,top;ll val;ans=0;
	cnt[now]=1;q[now][1]=0;dp[now][1]=1;
	for(int i=1;i<=n;i++)
	{
		for(int k=1;k<=cnt[now];k++)
			q[now][k]<<=2;
		for(int j=1;j<=m;j++)
		{
			now^=1;cnt[now]=0;++mt;
			for(int k=1;k<=cnt[now^1];k++)
			{
				bit=q[now^1][k];val=dp[now^1][k];
				b1=(bit>>(j-1)*2)&3;
				b2=(bit>>(j*2))&3;
				if(!a[i][j])
				{
					if(b1==0 && b2==0)
						insert(bit,val);
					continue;
				}
				if(b1==0 && b2==0)
				{
					if(a[i+1][j] && a[i][j+1])
						insert(bit^mi[j-1]^(mi[j]*2),val);
				}
				else if(b1==0 && b2>0)
				{
					if(a[i+1][j])
						insert(bit^(b2*mi[j])^(b2*mi[j-1]),val);
					if(a[i][j+1])
						insert(bit,val);
				}
				else if(b1>0 && b2==0)
				{
					if(a[i+1][j])
						insert(bit,val);
					if(a[i][j+1])
						insert(bit^(b1*mi[j-1])^(b1*mi[j]),val);
				}
				else if(b1==2 && b2==1)
					insert(bit^(mi[j-1]*2)^mi[j],val);
				else if(b1==1 && b2==1)
				{
					top=1;
					for(int l=j+1;l<=m;l++)
					{
						if((bit>>(l*2))%4==1) top++;
						if((bit>>(l*2))%4==2) top--;
						if(top==0)
						{
							insert(bit^mi[j-1]^mi[j]^(3*mi[l]),val);
							break;	
						}
					}
				}
				else if(b1==2 && b2==2)
				{
					top=1;
					for(int l=j-2;l>=0;l--)
					{
						if((bit>>(l*2))%4==1) top--;
						if((bit>>(l*2))%4==2) top++;
						if(!top)
						{
							insert(bit^(2*mi[j-1])^(2*mi[j])^(3*mi[l]),val);
							break;
						}
					}
				}
				else if(b1==1 && b2==2)
				{
					if(i==n && j==m)
						ans+=val;
				}
			}
		}
	}
}

inline void print()
{
	printf("%lld\n",ans);
}

int main()
{
	mi[0]=1;
	for(int i=1;i<=13;i++)
		mi[i]=mi[i-1]<<2;
	while(~scanf("%d%d",&n,&m))
	{
		if(n==0 && m==0) break;
		prework();
		mainwork();
		print();
	}
	return 0;
}