POJ - 1324 Holedox Moving (状态压缩+BFS/A*)
有一个n*m(1<=n,m<=20)的网格图,图中有k堵墙和有一条长度为L(L<=8)的蛇,蛇在移动的过程中不能碰到自己的身体。求蛇移动到点(1,1)所需的最小步数。
显然用8个(x,y)来表示蛇的状态是不现实的(用哈希也很难存下,要么爆内存,要么超时),所以首先应当进行状态压缩。可以发现蛇的身体是连续的,因此可以用一个表示方向的向量来储存蛇身体的每个部分在它上个部分的哪个方向(只有头部用(x,y)表示),这样总状态数就变成了n*m*(2^((L-1)*2)),在可接受范围内了。
可优化的部分:
1.可以在图的边界上放上一层墙,避免出界判定。
2.如果在碰撞判定时把“蛇身”和“墙”等同的话,为了防止在撤销蛇身所造成的影响的同时把墙也一并消除,可以用b[x][y]++,--的方法来代替b[x][y]=1。
AC代码:(编码和解码真是个体力活)
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
const int N=+;
const int dx[]= {,,-,};
const int dy[]= {-,,,};
struct P {int x,y;} p[N],r[N];
int n,m,l,k,b[N][N],d[N][N][(<<)+],S,ka;
struct D {int x,y,S;}; int bfs() {
queue<D> q;
memset(d,-,sizeof d);
d[p[].x][p[].y][S]=,q.push({p[].x,p[].y,S});
while(!q.empty()) {
int x=q.front().x,y=q.front().y,S=q.front().S;
q.pop();
if(x==&&y==)return d[x][y][S];
r[]= {x,y};
for(int i=; i<l; ++i) {
int t=S>>((l-i-)*)&;
r[i]= {r[i-].x+dx[t],r[i-].y+dy[t]};
}
for(int i=; i<l; ++i)b[r[i].x][r[i].y]++;
for(int i=; i<; ++i) {
int xx=x+dx[i],yy=y+dy[i],SS=S>>|((i^)<<((l-)*));
if(!b[xx][yy]&&!~d[xx][yy][SS])d[xx][yy][SS]=d[x][y][S]+,q.push({xx,yy,SS});
}
for(int i=; i<l; ++i)b[r[i].x][r[i].y]--;
}
return -;
} int main() {
while(scanf("%d%d%d",&n,&m,&l)&&n) {
printf("Case %d: ",++ka);
for(int i=; i<l; ++i)scanf("%d%d",&p[i].x,&p[i].y);
S=;
for(int i=; i<l; ++i)
for(int j=; j<; ++j)if(p[i-].x+dx[j]==p[i].x&&p[i-].y+dy[j]==p[i].y) {
S=S<<|j;
break;
}
memset(b,,sizeof b);
for(int i=; i<=n+; ++i)b[i][]=b[i][m+]=;
for(int i=; i<=m+; ++i)b[][i]=b[n+][i]=;
scanf("%d",&k);
while(k--) {
int x,y;
scanf("%d%d",&x,&y);
b[x][y]=;
}
printf("%d\n",bfs());
}
return ;
}
bfs
也可以进一步用A*优化,以每个点到点(1,1)的距离为估价函数进行转移,速度大幅提升。
一般A*中优先队列的结点要额外设一个属性g来表示已走过的距离,但也可以省去这个g,用d[x][y][S]来代替g,这样如果要把bfs改成A*的话,只要写出求每个点到点(1,1)距离的bfs代码,然后把原来的bfs代码中的queue换成priority_queue,front换成top就ok了,非常方便。
另外可以优化的一个地方是如果一个点不能到达(1,1),那么直接返回-1就行,不用再浪费时间跑A*了。
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
const int N=+;
const int dx[]= {,,-,};
const int dy[]= {-,,,};
struct P {int x,y;} p[N],r[N];
int n,m,l,k,b[N][N],d[N][N][(<<)+],S,h[N][N],ka;
struct D {
int x,y,S;
bool operator<(const D& b)const {return d[x][y][S]+h[x][y]>d[b.x][b.y][b.S]+h[b.x][b.y];}
}; void bfs() {
queue<P> q;
memset(h,-,sizeof h);
h[][]=,q.push({,});
while(!q.empty()) {
int x=q.front().x,y=q.front().y;
q.pop();
for(int i=; i<; ++i) {
int xx=x+dx[i],yy=y+dy[i];
if(!~h[xx][yy]&&!b[xx][yy])h[xx][yy]=h[x][y]+,q.push({xx,yy});
}
}
} int Astar() {
if(!~h[p[].x][p[].y])return -;
priority_queue<D> q;
memset(d,-,sizeof d);
d[p[].x][p[].y][S]=,q.push({p[].x,p[].y,S});
while(!q.empty()) {
int x=q.top().x,y=q.top().y,S=q.top().S;
q.pop();
if(x==&&y==)return d[x][y][S];
r[]= {x,y};
for(int i=; i<l; ++i) {
int t=S>>((l-i-)*)&;
r[i]= {r[i-].x+dx[t],r[i-].y+dy[t]};
}
for(int i=; i<l; ++i)b[r[i].x][r[i].y]++;
for(int i=; i<; ++i) {
int xx=x+dx[i],yy=y+dy[i],SS=S>>|((i^)<<((l-)*));
if(!b[xx][yy]&&!~d[xx][yy][SS])d[xx][yy][SS]=d[x][y][S]+,q.push({xx,yy,SS});
}
for(int i=; i<l; ++i)b[r[i].x][r[i].y]--;
}
return -;
} int main() {
while(scanf("%d%d%d",&n,&m,&l)&&n) {
printf("Case %d: ",++ka);
for(int i=; i<l; ++i)scanf("%d%d",&p[i].x,&p[i].y);
S=;
for(int i=; i<l; ++i)
for(int j=; j<; ++j)if(p[i-].x+dx[j]==p[i].x&&p[i-].y+dy[j]==p[i].y) {
S=S<<|j;
break;
}
memset(b,,sizeof b);
for(int i=; i<=n+; ++i)b[i][]=b[i][m+]=;
for(int i=; i<=m+; ++i)b[][i]=b[n+][i]=;
scanf("%d",&k);
while(k--) {
int x,y;
scanf("%d%d",&x,&y);
b[x][y]=;
}
bfs();
printf("%d\n",Astar());
}
return ;
}
A*
POJ - 1324 Holedox Moving (状态压缩+BFS/A*)的更多相关文章
- POJ 1324 Holedox Moving (状压BFS)
POJ 1324 Holedox Moving (状压BFS) Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 18091 Acc ...
- poj 1324 Holedox Moving
poj 1324 Holedox Moving 题目地址: http://poj.org/problem?id=1324 题意: 给出一个矩阵中,一条贪吃蛇,占据L长度的格子, 另外有些格子是石头, ...
- POJ 1324 Holedox Moving 搜索
题目地址: http://poj.org/problem?id=1324 优先队列---A*的估价函数不能为蛇头到(1,1)的距离,这样会出错. 看了discuss,有大神说这题A*的估价函数为BFS ...
- POJ 3411 Paid Roads (状态压缩+BFS)
题意:有n座城市和m(1<=n,m<=10)条路.现在要从城市1到城市n.有些路是要收费的,从a城市到b城市,如果之前到过c城市,那么只要付P的钱, 如果没有去过就付R的钱.求的是最少要花 ...
- POJ 1753 Flip Game(状态压缩+BFS)
题目网址:http://poj.org/problem?id=1753 题目: Flip Game Description Flip game is played on a rectangular 4 ...
- poj 1324 状态压缩+bfs
http://poj.org/problem?id=1324 Holedox Moving Time Limit: 5000MS Memory Limit: 65536K Total Submis ...
- 胜利大逃亡(续)(状态压缩bfs)
胜利大逃亡(续) Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total S ...
- POJ 1753 Flip Game (状态压缩 bfs+位运算)
Flip game is played on a rectangular 4x4 field with two-sided pieces placed on each of its 16 square ...
- hdu 3681 Prison Break(状态压缩+bfs)
Problem Description Rompire . Now it’s time to escape, but Micheal# needs an optimal plan and he con ...
随机推荐
- (C#)ListView双击Item事件
/// <summary> /// 双击选择播放列表项进行播放 /// </summary> /// <param name="sender"> ...
- 流量分析系统---kafka集群部署
1.集群部署的基本流程 Storm上游数据源之Kakfa 下载安装包.解压安装包.修改配置文件.分发安装包.启动集群 2.基础环境准备 安装前的准备工作(zk集群已经部署完毕) 关闭防火墙 chk ...
- 服务中的 API 网关(API Gateway)
我们知道在微服务架构风格中,一个大应用被拆分成为了多个小的服务系统提供出来,这些小的系统他们可以自成体系,也就是说这些小系统可以拥有自己的数据库,框架甚至语言等,这些小系统通常以提供 Rest Api ...
- git 使用教程 --基础二
一:分支学习: branch称为分支,默认仅有一个名为master的分支.一般开发新功能流程为:开发新功能时会在分支dev上进行,开发完毕后再合并到master分支. 分支的作用: 创建分支:(创建分 ...
- @MarkFan 口语练习录音 20140401
Hi,everybody 对于未来,我只梦想最好的情况, 并定下最踏实的计划,而绝不花时间在无谓的担心上, 因为我知道,只要把我对自己的承诺付诸实践, 我的未来将不会只是一个梦…… 这是引用考拉小巫的 ...
- 使用Shell脚本查找程序对应的进程ID,并杀死进程
#!/bin/sh NAME='shell.php' echo $NAME ID=`ps -ef | grep "$NAME" | grep -v "$0" | ...
- 深入Struts2的过滤器FilterDispatcher--中文乱码及字符编码过滤器
引用 前几天在论坛上看到一篇帖子,是关于Struts2.0中文乱码的,楼主采用的是spring的字符编码过滤器(CharacterEncodingFilter)统一编码为GBK,前台提交表单数据到Ac ...
- jmeter ant 运行 提示Error occurred during initialization of VM
运行ant提示错误 网上找到的方法 将set HEAP= -Xms512m -Xmx1024m 改成set HEAP= -Xms512m -Xmx512m 保存后运行成功
- Delphi_按字节比较两个文件
1.界面 2.代码 procedure TForm1.btnSelectFile01Click(Sender: TObject); begin if OpenDialog1.Execute then ...
- OpenStack企业私有云新需求(1):Nova 虚机支持 GPU
作者:Sammy Liu 刘世民 本系列会介绍OpenStack 企业私有云的几个需求: GPU 支持 自动扩展(Auto-scaling)支持 混合云(Hybrid cloud)支持 物理机(Bar ...