BZOJ 2574: [Poi1999]Store-Keeper
Description
推箱子. \(n,m\leqslant 100\)
Sol
Tarjan+边双连通分量+BFS.
直接搜索的复杂度是 \(n^6\) 记录人的位置,箱子的位置和转移.
箱子的位置相当于一个障碍.
先灌水,把移动到箱子周围,每次的状态都是人在箱子旁边,这样复杂度就降为了 \(4n^4\) .
每次BFS分两步,改变方向和移动箱子.
如果箱子所在位置不是割点,那么往所有方向都联通.
否则想改变方向的话,两点必须是双连通才行.
这样复杂度就降低了...判断的时候可以预处理出来.
这样判断可以变成 \(O(1)\) .但是我没有这么写,直接暴力判断两点是否有相同的双连通分量.
Code
#include<cstdio>
#include<cstring>
#include<utility>
#include<queue>
#include<vector>
#include<iostream>
using namespace std; #define debug(a) cout<<#a<<"="<<a<<" "
#define mpr make_pair
#define x first
#define y second
typedef pair< int,int > pr;
const int N = 105; int n,m,dx,dy,sx,sy,tx,ty,ans;
int id[N][N],a[N][N],isg[N*N],to[N*N][4],v[N][N],vis[N*N*4];
int dfsn[N*N],low[N*N],cnt,bcnt;
pr stk[N*N];int top;
vector< int > bl[N*N];
char s[N];
int mvx[]={ -1,1,0,0 };
int mvy[]={ 0,0,-1,1 };
struct S{ int x,y,r,d; };
queue< S > q; void Print(int a[N][N]){
cout<<"------------------------------"<<endl;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) printf("%d%c",a[i][j]," \n"[j==m]);
}
void insert(int x,int b){
for(int i=0,lim=bl[x].size();i<lim;i++) if(bl[x][i] == b) return;
bl[x].push_back(b);
}
void Tarjan(int u,int fa){
low[u]=dfsn[u]=++cnt;
for(int i=0,xx,yy;i<4;i++) if(to[u][i] && to[u][i]!=fa){
if(!dfsn[to[u][i]]){
stk[++top]=mpr(u,to[u][i]);
Tarjan(to[u][i],u);
low[u]=min(low[u],low[to[u][i]]);
if(dfsn[u]<=low[to[u][i]]){
++bcnt,isg[u]=1;
do{
xx=stk[top].x,yy=stk[top--].y;
insert(xx,bcnt),insert(yy,bcnt);
}while(!((xx==u && yy==to[u][i]) || (xx==to[u][i] && yy==u)));
}
}else low[u]=min(low[u],dfsn[to[u][i]]);
}
}
int Move(int x,int y,int r1,int r2){
if(!a[x+mvx[r2]][y+mvy[r2]]) return 0;
int now=id[x+mvx[r1]][y+mvy[r1]],gt=id[x+mvx[r2]][y+mvy[r2]];
if(!isg[id[x][y]]) return 1;
else{
for(int i=0,j,limi=bl[now].size(),limj=bl[gt].size();i<limi;i++)
for(j=0;j<limj;j++) if(bl[now][i] == bl[gt][j]) return 1;
}return 0;
}
void Fill(int fx,int fy){
if(fx == sx && fy == sy) return;
v[fx][fy]=1;
for(int i=0;i<4;i++) if(a[fx+mvx[i]][fy+mvy[i]] && !v[fx+mvx[i]][fy+mvy[i]] && (fx+mvx[i]!=sx || fy+mvy[i]!=sy))
Fill(fx+mvx[i],fy+mvy[i]);
}
int BFS(){
Fill(dx,dy);
// Print(v);
for(int i=0;i<4;i++) if(a[sx+mvx[i]][sy+mvy[i]]){
// debug(sx+mvx[i]),debug(sy+mvy[i])<<endl;
if(v[sx+mvx[i]][sy+mvy[i]]) vis[id[sx][sy]<<2|i]=1,q.push((S){ sx,sy,i,0 });
}
for(S now;!q.empty();){
now=q.front(),q.pop();
// cout<<now.x<<" "<<now.y<<" "<<now.r<<" "<<now.d<<endl;
if(now.x == tx && now.y == ty) return ans=now.d,1;
for(int i=0;i<4;i++){
if(now.r!=i && Move(now.x,now.y,now.r,i) && !vis[id[now.x][now.y]<<2|i]){
vis[id[now.x][now.y]<<2|i]=1;
if(!vis[id[now.x+mvx[i^1]][now.y+mvy[i^1]]<<2|i] && a[now.x+mvx[i^1]][now.y+mvy[i^1]]){
vis[id[now.x+mvx[i^1]][now.y+mvy[i^1]]<<2|i]=1;
q.push((S){ now.x+mvx[i^1],now.y+mvy[i^1],i,now.d+1 });
}
}else if(now.r == i){
if(!vis[id[now.x+mvx[i^1]][now.y+mvy[i^1]]<<2|i] && a[now.x+mvx[i^1]][now.y+mvy[i^1]]){
vis[id[now.x+mvx[i^1]][now.y+mvy[i^1]]<<2|i]=1;
q.push((S){ now.x+mvx[i^1],now.y+mvy[i^1],i,now.d+1 });
}
}
}
}return 0;
}
int main(){
// freopen("in.in","r",stdin);
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%s",s+1);
for(int j=1;j<=m;j++){
id[i][j]=i*m+j;
if(s[j] == 'S') a[i][j]=0;
else a[i][j]=1;
if(s[j] == 'M') dx=i,dy=j;
if(s[j] == 'P') sx=i,sy=j;
if(s[j] == 'K') tx=i,ty=j;
}
}
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++){
if(a[i-1][j]) to[id[i][j]][0]=id[i-1][j];
if(a[i][j-1]) to[id[i][j]][1]=id[i][j-1];
if(a[i+1][j]) to[id[i][j]][2]=id[i+1][j];
if(a[i][j+1]) to[id[i][j]][3]=id[i][j+1];
}
Tarjan(id[tx][ty],0); // Print(a); if(BFS()) printf("%d\n",ans);
else puts("NO");
return 0;
}
BZOJ 2574: [Poi1999]Store-Keeper的更多相关文章
- BZOJ 2933([Poi1999]地图-区间Dp)
2933: [Poi1999]地图 Time Limit: 1 Sec Memory Limit: 128 MB Submit: 7 Solved: 7 [ Submit][ Status] ...
- BZOJ 2929: [Poi1999]洞穴攀行
2929: [Poi1999]洞穴攀行 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 351 Solved: 195[Submit][Status][ ...
- bzoj 2929 [Poi1999]洞穴攀行 网络流
2929: [Poi1999]洞穴攀行 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 499 Solved: 295[Submit][Status][ ...
- 【刷题】BZOJ 2935 [Poi1999]原始生物
Description 原始生物的遗传密码是一个自然数的序列K=(a1,...,an).原始生物的特征是指在遗传密码中连续出现的数对(l,r),即存在自然数i使得l=ai且r=ai+1.在原始生物的遗 ...
- bzoj 2927: [Poi1999]多边形之战【博弈论】
先手必胜状态是黑三角在边上 然后其他情况脑补一下,n为偶数先手必胜,可以理解为从某一边取,先手总有办法让后手取得一边有奇数个 #include<iostream> #include< ...
- BZOJ 2927: [Poi1999]多边形之战 (博弈)
题意 有一个凸多边形,顶点编号逆时针从0到n-1.现在这个n边形被剖分成n-2个三角形,给出这n-2个三角形的顶点,保证这是用n-3条不交叉的对角线划分出来的.现在第一个三角形是黑色,其他都是白色.两 ...
- bzoj 2935 [Poi1999]原始生物——欧拉回路思路!
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2935 有向图用最小的路径(==总点数最少)覆盖所有边. 完了完了我居然连1999年的题都做不 ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- 【BZOJ 2802】 2802: [Poi2012]Warehouse Store (贪心)
2802: [Poi2012]Warehouse Store Description 有一家专卖一种商品的店,考虑连续的n天.第i天上午会进货Ai件商品,中午的时候会有顾客需要购买Bi件商品,可以选择 ...
随机推荐
- JavaScript的闭包原理
什么是js(JavaScript)的闭包原理,有什么作用? 一.定义 官方解释:闭包是一个拥有许多变量和绑定了这些变量的环境的表达式(通常是一个函数),因而这些变量也是该表达式的一部分. 个人的理解是 ...
- json跨域
很有意思的两种连接 ,效果相同. 不同之处: aehyok({"result":"我是远程js带来的数据"}); <script type="t ...
- Docker学习笔记 — 配置国内免费registry mirror
Docker学习笔记 — 配置国内免费registry mirror Docker学习笔记 — 配置国内免费registry mirror
- ApiHelper
public static class ApiHelper { public static ReslultInfo<K> Try<K>(Action<ReslultInf ...
- OC-成员变量的作用域
#import <Foundation/Foundation.h> @interface Person : NSObject { int _no; @public // 在任何地方都能直接 ...
- 总结Spring、Hibernate、Struts2官网下载jar文件
一直以来只知道搭SSH需要jar文件,作为学习的目的,最好的做法是自己亲自动手去官网下.不过官网都是英文,没耐心一般很难找到下载入口,更何 况版本的变化也导致不同版本jar文件有些不一样,让新手很容易 ...
- JQuery-EasyUI与EXTjs有什么区别?
一.ExtJS1.ExtJS可以用来开发RIA也即富客户端的AJAX应用,是一个用javascript写的,主要用于创建前端用户界面,是一个与后台技术无关的前端ajax框架.因此,可以把ExtJS用在 ...
- Properties类的使用方法
它提供了几个主要的方法: 1. getProperty ( String key),用指定的键在此属性列表中搜索属性.也就是通过参数 key ,得到 key 所对应的 value. 2. load ( ...
- spring配置文件详解--真的蛮详细
spring配置文件详解--真的蛮详细 转自: http://book.51cto.com/art/201004/193743.htm 此处详细的为我们讲解了spring2.5的实现原理,感觉非常 ...
- 摄像头/光驱故障:由于其配置信息(注册表中的)不完整或已损坏,Windows 无法启动这个硬件设备。 (代码 19)
原因:出现这条提示一般都是因为注册表错误引起的,重装驱动或应用程序无法解决. 设备管理器中相关设备-故障如图: 以下方法可以解决问题: 1.在任务栏上点“开始”菜单-“运行”,输入”regedit“ ...