bzoj1814 Ural 1519 Formula 1（插头DP）

对插头DP的理解还不是很透彻。

先说一下肤浅的理解吧。

插头DP使用范围：指数级复杂度，且适用于解决网格图连通性问题，如哈密顿回路等问题。插头一般指每相邻2个网格的接口。

题目难度：一般不可做。

使用三进制状态，用0表示没有插头，1表示左括号插头，2表示右括号插头，而由于位运算常数特别小，可以采用四进制+手工hash的做法。处理好状态与编号的对应关系后进行dp，复杂度就只与合法状态数有关了，时间复杂度O(snm^2)，其中s为合法状态数（不超过42000），而第二维的m显然也是不满的，因此可以通过本题。

然后递推就是大力讨论情况，详细见代码。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=3e5+,mod=;
int n,m,D,ex,ey,tot[],mp[][],pw[],hd[N],nxt[N],a[][N];
ll ans,f[][N];
char str[];
void add(int x,ll y)
{
    int u=x%mod+;
    for(int i=hd[u];i;i=nxt[i])if(a[D][i]==x){f[D][i]+=y;return;}
    nxt[++tot[D]]=hd[u],hd[u]=tot[D];
    a[D][tot[D]]=x,f[D][tot[D]]=y;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)
    {
        scanf("%s",str+);
        for(int j=;j<=m;j++)
        if(str[j]=='.'){mp[i][j]=,ex=i,ey=j;}
    }
    pw[]=;
    for(int i=;i<=max(n,m);i++)pw[i]=pw[i-]*;
    a[][]=,tot[]=f[][]=;
    for(int i=;i<=n;i++)
    {
        for(int j=;j<=tot[D];j++)a[D][j]*=;
        for(int j=;j<=m;j++)
        {
            memset(hd,,sizeof hd);
            tot[D^=]=;
            for(int k=;k<=tot[D^];k++)
            {
                int S=a[D^][k],b1=(S>>(j*-))%,b2=(S>>(j*))%;
                ll p=f[D^][k];
                if(!mp[i][j])
                {
                    if(!b1&&!b2)add(S,p);
                }
                else if(!b1&&!b2)
                {
                    if(mp[i+][j]&&mp[i][j+])add(S+pw[j-]+*pw[j],p);
                }
                else if(!b1&&b2)
                {
                    if(mp[i][j+])add(S,p);
                    if(mp[i+][j])add(S-b2*pw[j]+b2*pw[j-],p);
                }
                else if(b1&&!b2)
                {
                    if(mp[i+][j])add(S,p);
                    if(mp[i][j+])add(S-b1*pw[j-]+b1*pw[j],p);
                }
                else if(b1==b2)
                {
                    if(b1==)
                    {
                        int tmp=;
                        for(int l=j+;l<=m;l++)
                        {
                            if((S>>(l*))%==)tmp++;
                            if((S>>(l*))%==)tmp--;
                            if(!tmp){add(S-pw[l]-pw[j-]-pw[j],p);break;}
                        }
                    }
                    if(b1==)
                    {
                        int tmp=;
                        for(int l=j-;l>=;l--)
                        {
                            if((S>>(l*))%==)tmp--;
                            if((S>>(l*))%==)tmp++;
                            if(!tmp){add(S+pw[l]-*pw[j-]-*pw[j],p);break;}
                        }
                    }
                }
                else if(b1==&&b2==)add(S-*pw[j-]-pw[j],p);
                else if(i==ex&&j==ey)ans+=p;
            }
        }
    }
    printf("%lld",ans);
}