【poj1739】 Tony's Tour
http://poj.org/problem?id=1739 (题目链接)
题意
给出一个n*m的地图,有些是障碍。问从左下角走遍所有非障碍格子一次且仅一次最终到达右下角的路径方案数。
Solution
插头dp。
我们给地图的再加上2行:
.####.
......
那么最后就变成了求一个回路了,思路参见cdq论文。
UPD:并不是每种转移都需要重新求一遍最小表示法,但是那样写虽然常数略大,代码会整洁很多。虽然有强迫症,但是为了常数,还是把代码更新一下吧。
细节
细节见代码。
代码
// poj1739
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define LL long long
#define MOD 4001
#define inf 2147483640
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxd=15,maxs=100010,maxh=4010;
int head[maxh],next[maxs];
int n,m,code[maxd],t[maxd],a[maxd][maxd],size[2],s[2][maxs];
char ch[maxd];
int tot[2][maxs]; void decode(int st) { //解码
for (int i=m;i>=0;i--) {
code[i]=st&7;
st>>=3;
}
}
int encode(int op) { //编码
int cnt=0;int st=0;
if (op) { //有时并不需要重新求最小表示法
memset(t,-1,sizeof(t));t[0]=0;
for (int i=0;i<=m;i++) {
if (t[code[i]]==-1) t[code[i]]=++cnt;
code[i]=t[code[i]];
}
}
for (int i=0;i<=m;i++) st=st<<3|code[i];
return st;
}
void shift() { //换行
for (int i=m;i;i--) code[i]=code[i-1];
code[0]=0;
}
void add(int op,int p,int num) { //更新hash表
int tmp=encode(op),id=tmp%MOD;
for (int i=head[id];i;i=next[i])
if (s[p][i]==tmp) {tot[p][i]+=num;return;}
next[++size[p]]=head[id];s[p][size[p]]=tmp;tot[p][size[p]]=num;
head[id]=size[p];
}
int main() {
while (scanf("%d%d",&n,&m)!=EOF && n && m) {
memset(a,0,sizeof(a));
for (int i=1;i<=n;i++) {
scanf("%s",ch+1);
for (int j=1;j<=m;j++) a[i][j]=ch[j]=='.';
}
a[n+1][1]=a[n+1][m]=1;
for (int i=1;i<=m;i++) a[n+2][i]=1;
n+=2; //加行
int ex=n,ey=m; //终止格显然是无障碍的 int p=0;
tot[0][1]=1;size[0]=1;s[0][1]=0;
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++) {
p^=1; //滚动
size[p]=0;
memset(head,0,sizeof(head)); //清空hash表
for (int k=1;k<=size[p^1];k++) {
decode(s[p^1][k]); //解码
int left=code[j-1],up=code[j]; //左插头,上插头
if (!a[i][j]) { //该格子为障碍格
code[j-1]=code[j]=0;
if (j==m) shift();
add(0,p,tot[p^1][k]);
continue;
}
if (left && up) { //左、上插头都存在,合并连通块
if (left==up) { //即将合并的两部分在同一个联通分量中
if (i==ex && j==ey) { //只能出现在最后一个格子
code[j]=code[j-1]=0;
if (j==m) shift();
add(0,p,tot[p^1][k]);
}
}
else { //不在同一个联通分量,合并
code[j-1]=code[j]=0;
for (int l=0;l<=m;l++) if (code[l]==left) code[l]=up;
if (j==m) shift();
add(1,p,tot[p^1][k]);
}
}
else if (left || up) { //左、上插头存在一个,连通块不变
int tmp;
if (left) tmp=left;
else tmp=up;
if (a[i][j+1]) { //插头插向右边的格子,这一步一定要先写,因为后面那步需要shift,会破坏code
code[j-1]=0;code[j]=tmp;
add(0,p,tot[p^1][k]);
}
if (a[i+1][j]) { //插头插向下边的格子
code[j-1]=tmp;code[j]=0;
if (j==m) shift();
add(0,p,tot[p^1][k]);
}
}
else { //左、上插头都不存在。新建连通块
if (a[i][j+1] && a[i+1][j]) {
code[j-1]=code[j]=10; //随意设一个一定不会与之前连通块重复的编号
add(1,p,tot[p^1][k]);
}
}
}
}
int ans=0;
for (int i=1;i<=size[p];i++) ans+=tot[p][i];
printf("%d\n",ans);
}
return 0;
}
【poj1739】 Tony's Tour的更多相关文章
- 【POJ】【1739】Tony's Tour
插头DP 楼教主男人八题之一! 要求从左下角走到右下角的哈密顿路径数量. 啊嘞,我只会求哈密顿回路啊……这可怎么搞…… 容易想到:要是把起点和重点直接连上就变成一条回路了……那么我们就连一下~ 我们可 ...
- 【SPOJ】1825. Free tour II(点分治)
http://www.spoj.com/problems/FTOUR2/ 先前看了一会题解就自己yy出来了...对拍过后交tle.................. 自己造了下大数据........t ...
- 【Codeforces858F】Wizard's Tour [构造]
Wizard's Tour Time Limit: 50 Sec Memory Limit: 512 MB Description Input Output Sample Input 4 5 1 2 ...
- 【读书笔记】A Swift Tour
素材:A Swift Tour 推荐下载Playground:Download Playground objc 自己较为熟悉,想熟悉下风头正劲的 swift.就先从官方的入门手册开始撸. 每一小节,我 ...
- 【BZOJ3060】[Poi2012]Tour de Byteotia 并查集
[BZOJ3060][Poi2012]Tour de Byteotia Description 给定一个n个点m条边的无向图,问最少删掉多少条边能使得编号小于等于k的点都不在环上. Input ...
- 【CF1053E】Euler tour
[CF1053E]Euler tour 题面 CF 洛谷 大概意思是你有一棵树,然而你并不知道这棵树是啥.给你一个确定了一些位置的欧拉序(就是\(ST\)表求\(LCA\)的那个序列),问你是否存在一 ...
- 【POJ】【1637】Sightseeing tour
网络流/最大流 愚人节快乐XD 这题是给一个混合图(既有有向边又有无向边),让你判断是否有欧拉回路…… 我们知道如果一个[连通]图中每个节点都满足[入度=出度]那么就一定有欧拉回路…… 那么每条边都可 ...
- Python开发【前端】:JavaScript
JavaScript入门 JavaScript一种直译式脚本语言,是一种动态类型.弱类型.基于原型的语言,内置支持类型.它的解释器被称为JavaScript引擎,为浏览器的一部分,广泛用于客户端的脚本 ...
- 【python】函数之内置函数
Python基础 内置函数 今天来介绍一下Python解释器包含的一系列的内置函数,下面表格按字母顺序列出了内置函数: 下面就一一介绍一下内置函数的用法: 1.abs() 返回一个数值的绝对值,可以是 ...
随机推荐
- 【强烈推荐】XCODE的插件之王
有许多关于Xcode的插件,在这里强烈推荐的是Alcatraz插件.因为我们可以通过这个插件来安装其他插件 1.Alcatraz插件. Alcatraz是一个方便我们安装各种那个插件的插件.插件之王? ...
- iOS开发之功能模块--根据需求开发横向的子弹盒View
这个需求是本人工作开发中后期需求要添加的新功能,本人模仿UITableView的代理和数据源方法进行了第一阶段的开发.第二阶段是添加丰富的动画. 这个功能需求描述:能上传添加五个待选头像,五个头像分别 ...
- jQuery插件库代码分享 - 进阶者系列 - 学习者系列文章
这些天将原来在网上找的jQuery插件进行了下整理,特此将代码分享出来给大家. 见下图结构. 对目录结构进行了分类.这里是插件列表. 这里总共收集了20来个插件.还有下面未进行划分的. 下面是DEMO ...
- Ignite 配置更新Oracle JDBC Drive
如果使用Oracle 12C 作为Ignite 的Repository的话,在Repository Createion Wizard的配置过程中,会出现ORA-28040:No matc ...
- 聊下git merge --squash
你经常会面临着将dev分支或者很多零散的分支merge到一个公共release分支里. 但是有一种情况是需要你处理的,就是在你的dev的分支里有很多commit记录.而这些commit是无需在rele ...
- 初刷LeetCode的感受
自从上个月进入实验室的云安全项目组后,因为要接触到实际的代码,在实验室博士的建议下我们项目组的硕士开始刷LeetCode练习编程能力,保持每周抽空刷几道算法题.虽然刷的不多,到现在一共只刷了不到30题 ...
- shell 1到指定数累加
#!/bin/bash read -p "输入尾数:" a expr $(seq -s " + " $a) #seq命令可以指定生成一个数到另一个数之间的所有整 ...
- 005.nginx配置文件
1.替换nginx主配置文件 通过前面的配置,LNMP的环境已经搭建完成,现在我们替换nginx配置文件: [root@huh ~]# cd /usr/local/nginx/conf/[root@h ...
- Windows7开机登录界面背景图修改
首先,准备一张文件大小不超过250K的.jpg格式图片,文件名为:BackgroundDefault,图片分辨率建议与屏幕分辨率一致,因为这张背景无法像壁纸一样拉伸或者填充的. 打开这个路径,[C:\ ...
- Windows共享作为公司文件服务器的案例
1.目录结构 → 主管 部门 → 员工 → Public 2.实现效果 每个部门一个目录 部门主管可以访问自己和部门员工的目录 部门员工只可访问自己的目录 公共目录Public部门所有人都可访问 3. ...