题目链接

https://www.lydsy.com/JudgeOnline/problem.php?id=1814

题解

插头dp，枚举边界线的插头，分类讨论转移点上插头和左插头的情况，来转移这个点的插头。

代码

#include <cstdio>
#include <cstring>
#include <algorithm>const int maxn=20;
const int hash_mod=1007;
const int maxd=maxn*maxn;
const int maxk=50000;int read()
{int x=0,f=1;char ch=getchar();while((ch<'0')||(ch>'9')){if(ch=='-'){f=-f;}ch=getchar();}while((ch>='0')&&(ch<='9')){x=x*10+ch-'0';ch=getchar();}return x*f;
}struct func
{int x;long long y;
};struct hash_table
{int now[hash_mod+10],pre[maxk+2],cnt;func val[maxk+2];int clear(){cnt=0;memset(now,0,sizeof now);return 0;}long long add(int s,long long v){int trans=s%hash_mod,j=now[trans];while(j){if(val[j].x==s){val[j].y+=v;return val[j].y;}j=pre[j];}val[++cnt].x=s;val[cnt].y=v;pre[cnt]=now[trans];now[trans]=cnt;return v;}func get(int x){return val[x];}
};hash_table ht[2];
int n,m,rv[maxn+4],v[maxn+2][maxd+2],now,last;
char s[maxn+4];int decode(int x)
{for(int i=0; i<=m; ++i){rv[i]=x&3;x>>=2;}return 0;
}int encode()
{int res=0;for(int i=m; ~i; --i){res=(res<<2)+rv[i];}return res;
}int findl(int pos)
{int now=-1;while(1){if(rv[pos]){now+=(rv[pos]==1)?1:-1;}if(!now){return pos;}--pos;}
}int findr(int pos)
{int now=1;while(1){if(rv[pos]){now+=(rv[pos]==1)?1:-1;}if(!now){return pos;}++pos;}
}int expand(int x,int y,func from)
{decode(from.x);int west=rv[m],north=rv[y];if((!y)&&west){return 0;}if(v[x][y]){if(west||north){return 0;}ht[now].add(encode(),from.y);}else if(west&&north){rv[m]=rv[y]=0;if(west<north){int flag=0;for(int i=0; i<=m; ++i){if(rv[i]){flag=1;break;}}if(!((!flag)&&(x==n))){return 0;}}if(west==north){if(west==1){rv[findr(y+1)]=1;}else{rv[findl(y-1)]=2;}}ht[now].add(encode(),from.y);}else if((!west)&&(!north)){rv[y]=1;rv[m]=2;ht[now].add(encode(),from.y);}else if((!west)&&north){ht[now].add(encode(),from.y);rv[m]=north;rv[y]=0;ht[now].add(encode(),from.y);}else if(west&&(!north)){ht[now].add(encode(),from.y);rv[y]=west;rv[m]=0;ht[now].add(encode(),from.y);}return 0;
}int work()
{now=0;last=1;ht[0].clear();ht[0].add(0,1);for(int i=1; i<=n; ++i){for(int j=0; j<m; ++j){std::swap(now,last);ht[now].clear();for(int k=1; k<=ht[last].cnt; ++k){expand(i,j,ht[last].get(k));}}}return 0;
}int main()
{n=read();m=read();for(int i=1; i<=n; ++i){scanf("%s",s);for(int j=0; j<m; ++j){if(s[j]=='*'){v[i][j]=1;}}}work();printf("%lld\n",ht[now].add(0,0));return 0;
}

BZOJ 1814 Ural 1519 Formula 1相关推荐

bzoj 1814: Ural 1519 Formula 1【插头dp】
设f[i][j][s]为轮廓线推到格子(i,j),状态为s的方案数括号表示一段线的左端和右端,表示成左括号和右括号,状压的时候用1和2表示,0表示已经闭合下面的蓝线是黄色格子的轮廓线,dp转移要把 ...
bzoj1814 Ural 1519 Formula 1(插头dp模板题)
1814: Ural 1519 Formula 1 Time Limit: 1 Sec Memory Limit: 64 MB Submit: 924 Solved: 351 [Submit][S ...
【BZOJ1814】Ural 1519 Formula 1 插头DP
[BZOJ1814]Ural 1519 Formula 1 题意:一个 m * n 的棋盘,有的格子存在障碍,求经过所有非障碍格子的哈密顿回路个数.(n,m<=12) 题解:插头DP板子题,刷板 ...
Ural 1519. Formula 1 优美的插头DP
今天早上学了插头DP的思想和最基础的应用,中午就开始敲了,岐哥说第一次写不要看别人代码,利用自己的理解一点点得写出来,这样才锻炼代码能力!于是下午慢慢地构思轮廓,一点点地敲出主体代码,其实是很磨蹭的, ...
bzoj 4451 : [Cerc2015]Frightful Formula FFT
4451: [Cerc2015]Frightful Formula Time Limit: 10 Sec Memory Limit: 64 MB Submit: 177 Solved: 57 [S ...
[COGS 1512][URAL 1519]一级方程式赛车
按格插头DP #include <iostream> #include <cstdio> #include <cstring> #include <algor ...
BZOJ 4451: [Cerc2015]Frightful Formula
Description 给你一个n*n矩阵的第一行和第一列,其余的数通过如下公式推出: \(f_{i,j}=a\times f_{i,j-1}+b\times f_{i-1,j}+c\) 求\(f_ ...
bzoj 4451: [Cerc2015]Frightful Formula 数学+排列组合
题意给你一个n*n矩阵的第一行和第一列,其余的数通过如下公式推出: F[i,j]=a*f[i,j-1]+b*f[i-1,j]+c 求f[n][n]%(10^6+3) 2<=n<=2000 ...
[FFT || 递推] BZOJ 4451 [Cerc2015]Frightful Formula
f1,i f_{1,i} 的贡献是 f1,i∗an−i∗bn−1∗Cn−i2n−i−2 f_{1,i}*a^{n-i}*b^{n-1}*C_{2n-i-2}^{n-i} fi,1 f_{i,1} 同理 ...

BZOJ 1814 Ural 1519 Formula 1

题目链接

题解

代码

BZOJ 1814 Ural 1519 Formula 1相关推荐

最新文章

热门文章