POJ 1739 Tony's Tour

题意：给定n*m的带障碍格子图，求从左下角走到右下角的哈密顿路的总数。

题解：更进一步的插头DP，主要解法可以参照我的上篇日志：http://blog.csdn.net/tmeteorj/article/details/8635090。

改动之处在于单向插头，上篇中用01和10分别代表左右括号，这道题则需要加11作为单向括号（因为它在括号序列状态中是没有与之匹配的括号的）

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int MAX=10007;typedef long long LL;struct HT{    int head,hash[MAX],nxt[MAX];    LL cnt[MAX];    void init()    {        head=-1;        memset(hash,-1,sizeof(hash));        memset(nxt,-1,sizeof(nxt));        memset(cnt,0,sizeof(cnt));    }    void insert(int st,LL ct)    {        LL pos=st%MAX;        LL ad=1;        while(hash[pos]!=st&&hash[pos]!=-1)        {            pos+=ad*ad;            ad++;            if(pos>=MAX)pos%=MAX;        }        if(hash[pos]==-1)        {            hash[pos]=st;            cnt[pos]=ct;            nxt[pos]=head;            head=pos;        }        else            cnt[pos]+=ct;    }} dp[2];int n,m;char mp[14][14];int findright(int st,int pos){    int tp,p=pos;    for(int dep=1; dep!=0; p+=1)    {        tp=((st>>((m-p+1)*2))&3)%4;        switch(tp)        {        case 0:            break;        case 1:            dep++;            break;        case 2:            dep--;            break;        }        if(dep==0)            break;    }    return 3<<((m-p+1)*2);}int findleft(int st,int pos){    int tp,p=pos;    for(int dep=1; dep!=0; p-=1)    {        tp=((st>>((m-p+1)*2))&3)%4;        switch(tp)        {        case 0:            break;        case 1:            dep--;            break;        case 2:            dep++;            break;        }        if(dep==0)            break;    }    return 3<<((m-p+1)*2);}LL getans(int fg){    LL pos=0,ad=1;    while(dp[fg].hash[pos]!=0&&dp[fg].hash[pos]!=-1)    {        pos+=ad*ad;        ad++;        if(pos>=MAX)pos%=MAX;    }    if(dp[fg].hash[pos]==-1)    {        return 0ll;    }    else        return dp[fg].cnt[pos];}int main(){    while(scanf("%d%d",&n,&m),n||m)    {        for(int i=1; i<=n; i++)scanf(" %s",mp[i]+1);        int fg1=1,fg2,ex,ey;        fg1=0,fg2=1;        dp[fg2].init();        dp[fg2].insert(0,1ll);        for(int i=1; i<=n; i++)        {            fg1=!fg1;            fg2=!fg2;            dp[fg2].init();            for(int tp=dp[fg1].head; tp!=-1; tp=dp[fg1].nxt[tp])            {                if(dp[fg1].hash[tp]&3)                    continue;                dp[fg2].insert(dp[fg1].hash[tp]>>2,dp[fg1].cnt[tp]);            }            for(int j=0; j<m; j++)            {                fg1=!fg1;                fg2=!fg2;                dp[fg2].init();                int a=3<<(2*(m-j)),b=3<<(2*(m-j-1));                for(int pos=dp[fg1].head; pos!=-1; pos=dp[fg1].nxt[pos])                {                    int st=dp[fg1].hash[pos];                    LL ct=dp[fg1].cnt[pos];                    int x=((st&a)>>(2*(m-j)))%4,y=((st&b)>>(2*(m-j-1)))%4;                    if(i==n)                    {                        if(j==0)                        {                            if(y)                                dp[fg2].insert(st&(~b)|findright(st,j+3),ct);                            else                                dp[fg2].insert(st|b,ct);                            continue;                        }                        else if(j==m-1)                        {                            if((x==3||y==3)&(x!=3||y!=3))                                dp[fg2].insert(st&(~a)&(~b),ct);                            continue;                        }                    }                    if(mp[i][j+1]=='#')                    {                        if(x||y)                            continue;                        dp[fg2].insert(st,ct);                    }                    else if(x)                    {                        if(y)                        {                            if(x==1&&y==1)                                dp[fg2].insert((st&(~a)&(~b))^findright(st,j+3),ct);                            else if(x==2&&y==1)                                dp[fg2].insert(st&(~a)&(~b),ct);                            else if(x==2&&y==2)                                dp[fg2].insert((st&(~a)&(~b))^findleft(st,j),ct);                            else if(x==3)                                dp[fg2].insert((st&(~a)&(~b))|findright(st,j+3),ct);                            else if(y==3)                                dp[fg2].insert((st&(~a)&(~b))|findleft(st,j),ct);                        }                        else                        {                            dp[fg2].insert(st,ct);                            dp[fg2].insert(st&(~a)|(x<<(2*(m-j-1))),ct);                        }                    }                    else if(y)                    {                        dp[fg2].insert(st,ct);                        dp[fg2].insert(st&(~b)|(y<<(2*(m-j))),ct);                    }                    else                        dp[fg2].insert(st|(1<<(2*(m-j)))|(2<<(2*(m-j-1))),ct);                }            }        }        printf("%lld\n",getans(fg2));    }    return 0;}