URAL 1519Formula 1

链接：http://ln1.vjudge.net/contest/3#overview

题意：给定一个矩阵，每个格子为空格或者有障碍物。求覆盖到所有的空白格的哈密顿回路的方案数。

分析：插头dp，换了一种更模块化的写法，手写hash，我用的是4进制的括号匹配表示。因为要求只能有一条回路，那么我们只要确定()括号合并只能在最后一个空格即可。这样就不会产生多余的回路啦。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=15;
const int M=1e5+7;
const int HASHSIZE=3e5+9;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int INF=1000000010;
const ll MAX=1ll<<55;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
const int hashsize=1e5+7,statesize=3e5+10;
char s[N];
int w[2*N],bit[2*N];
int n,m,ex,ey,pre,now,ma[N][N];
int tot[2],head[hashsize],nex[hashsize];
ll ans,state[2][statesize],dp[2][statesize];
void init() {
    memset(ma,0,sizeof(ma));
    ans=now=tot[0]=0;
    state[now][0]=0;
    dp[now][0]=tot[now]=1;
}
void Hash(ll s,ll num) {
    int i,pos=s%hashsize;
    for (i=head[pos];~i;i=nex[i])
    if (state[now][i]==s) { dp[now][i]+=num;return ; }
    state[now][tot[now]]=s;dp[now][tot[now]]=num;
    nex[tot[now]]=head[pos];head[pos]=tot[now]++;
}
void DP() {
    int i,j,k,t,b,p,q;ll s,num;
    for (i=1;i<=n;i++) {
        for (k=0;k<tot[now];k++) state[now][k]<<=2;
        for (j=1;j<=m;j++) {
            pre=now;now^=1;tot[now]=0;
            memset(head,-1,sizeof(head));
            for (k=0;k<tot[pre];k++) {
                s=state[pre][k];num=dp[pre][k];
                p=s>>w[j-1]&3;q=s>>w[j]&3;
                if (!ma[i][j]) {
                    if (!p&&!q) Hash(s,num);
                } else if (!p&&!q) {
                    if (ma[i][j+1]&&ma[i+1][j]) {
                        s^=bit[w[j-1]]^bit[w[j]+1];Hash(s,num);
                    }
                } else if (p&&!q) {
                    if (ma[i+1][j]) Hash(s,num);
                    if (ma[i][j+1]) {
                        s^=bit[w[j-1]+p-1]^bit[w[j]+p-1];Hash(s,num);
                    }
                } else if (!p&&q) {
                    if (ma[i][j+1]) Hash(s,num);
                    if (ma[i+1][j]) {
                        s^=bit[w[j-1]+q-1]^bit[w[j]+q-1];Hash(s,num);
                    }
                } else if (p==1&&q==1) {
                    for (b=1,t=j+1;t<=m;t++) {
                        if ((s>>w[t]&3)==2) b--;
                        if ((s>>w[t]&3)==1) b++;
                        if (!b) {
                            s^=bit[w[t]]^bit[w[t]+1];break ;
                        }
                    }
                    s^=bit[w[j-1]]^bit[w[j]];Hash(s,num);
                } else if (p==1&&q==2) {
                    if (i==ex&&j==ey) ans+=num;
                } else if (p==2&&q==1) {
                    s^=bit[w[j-1]+1]^bit[w[j]];Hash(s,num);
                } else {
                    for (b=1,t=j-2;t>=0;t--) {
                        if ((s>>w[t]&3)==2) b++;
                        if ((s>>w[t]&3)==1) b--;
                        if (!b) {
                            s^=bit[w[t]]^bit[w[t]+1];break ;
                        }
                    }
                    s^=bit[w[j-1]+1]^bit[w[j]+1];Hash(s,num);
                }
            }
        }
    }
}
void solve() {
    init();
    for (int i=1;i<=n;i++) {
        scanf("%s", s+1);
        for (int j=1;j<=m;j++) {
            ma[i][j]=s[j]=='.';
            if (ma[i][j]) ex=i,ey=j;
        }
    }
    DP();
    printf("%lld\n", ans);
}
int main()
{
    for (int i=0;i<2*N;i++) w[i]=i<<1,bit[i]=1<<i;
    while (scanf("%d%d", &n, &m)!=EOF) solve();
    return 0;
}