bzoj 3125: CITY

最新推荐文章于 2018-04-28 07:30:00 发布

fyc_kabuto

最新推荐文章于 2018-04-28 07:30:00 发布

阅读量204

点赞数 1

分类专栏：插头dp

本文链接：https://blog.csdn.net/qq_36808030/article/details/79075775

版权

插头dp 专栏收录该内容

6 篇文章 0 订阅

订阅专栏

题意：

给一个地图，有些地方只能按某个方向通过，问哈密顿回路个数。

题解：

直接插头dp，钦定的地方特判一下就好了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const int mod=100037;
int n,m,map[15][15],X=-1,Y;
int hash[mod],list[2][mod],now=0,num[2];
LL ans=0,state[2][mod];
char str[15];
int get(int s,int p) {return (s>>((p-1)*2))&3;}
void change(int &s,int p,int c)
{
    s^=get(s,p)<<((p-1)*2);
    s^=c<<((p-1)*2);
}
void add(int st,LL sum)
{
    int s=st%mod;
    while(hash[s]!=-1&&list[now][hash[s]]!=st) (s+=1)%=mod;
    if(hash[s]==-1) hash[s]=++num[now],list[now][num[now]]=st,state[now][num[now]]=sum;
    else state[now][hash[s]]+=sum;
}
void solve()
{
    state[0][1]=num[0]=1;list[0][1]=ans=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            now^=1;num[now]=0;memset(hash,-1,sizeof(hash));
            for(int k=1;k<=num[!now];k++)
            {
                int st=list[!now][k],p=get(st,j),q=get(st,j+1);
                LL sum=state[!now][k];
                if(!map[i][j]) {if(!p&&!q) add(st,sum);continue;}
                if(map[i][j]==2)
                {
                    if(p&&!q&&map[i][j+1]) change(st,j,0),change(st,j+1,p),add(st,sum);
                    continue;
                }
                if(map[i][j]==3)
                {
                    if(!p&&q&&map[i+1][j]) change(st,j,q),change(st,j+1,0),add(st,sum);
                    continue;
                }
                if(!p&&!q)
                {
                    if(map[i+1][j]&&map[i][j+1])
                        change(st,j,1),change(st,j+1,2),add(st,sum);
                }
                else if(!p&&q)
                {
                    if(map[i][j+1]) add(st,sum);
                    if(map[i+1][j]) change(st,j,q),change(st,j+1,0),add(st,sum);
                }
                else if(p&&!q)
                {
                    if(map[i+1][j]) add(st,sum);
                    if(map[i][j+1]) change(st,j,0),change(st,j+1,p),add(st,sum);
                }
                else if(p==1&&q==2) {if(i==X&&j==Y) ans+=sum;}
                else if(p==2&&q==1) {change(st,j,0);change(st,j+1,0);add(st,sum);}
                else if(p==1&&q==1)
                {
                    int top=1;
                    for(int pos=j+2;pos<=m+1;pos++)
                    {
                        int tmp=get(st,pos);
                        if(tmp==1) top++;if(tmp==2) top--;
                        if(top==0) {change(st,j,0);change(st,j+1,0);change(st,pos,1);add(st,sum);break;}
                    }
                }
                else if(p==2&&q==2)
                {
                    int top=1;
                    for(int pos=j-1;pos;pos--)
                    {
                        int tmp=get(st,pos);
                        if(tmp==2) top++;if(tmp==1) top--;
                        if(top==0) {change(st,j,0);change(st,j+1,0);change(st,pos,2);add(st,sum);break;}
                    }
                }
            }
        }
        for(int j=1;j<=num[now];j++) list[now][j]<<=2;
    }
}
int qs(char s)
{
    if(s=='.') return 1;
    if(s=='-') return 2;
    if(s=='|') return 3;
    return 0;
}
int main()
{
    scanf("%d %d",&n,&m);
    memset(map,0,sizeof(map));
    for(int i=1;i<=n;i++)
    {
        scanf("%s",str+1);
        for(int j=1;j<=m;j++)
        {
            if(str[j]=='.') X=i,Y=j;
            map[i][j]=qs(str[j]);
        }
    }
    if(X==-1) {printf("0");return 0;}
    solve();printf("%lld",ans);
}