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的更多相关文章

  1. BZOJ 2933([Poi1999]地图-区间Dp)

    2933: [Poi1999]地图 Time Limit: 1 Sec   Memory Limit: 128 MB Submit: 7   Solved: 7 [ Submit][ Status] ...

  2. BZOJ 2929: [Poi1999]洞穴攀行

    2929: [Poi1999]洞穴攀行 Time Limit: 1 Sec  Memory Limit: 128 MBSubmit: 351  Solved: 195[Submit][Status][ ...

  3. bzoj 2929 [Poi1999]洞穴攀行 网络流

    2929: [Poi1999]洞穴攀行 Time Limit: 1 Sec  Memory Limit: 128 MBSubmit: 499  Solved: 295[Submit][Status][ ...

  4. 【刷题】BZOJ 2935 [Poi1999]原始生物

    Description 原始生物的遗传密码是一个自然数的序列K=(a1,...,an).原始生物的特征是指在遗传密码中连续出现的数对(l,r),即存在自然数i使得l=ai且r=ai+1.在原始生物的遗 ...

  5. bzoj 2927: [Poi1999]多边形之战【博弈论】

    先手必胜状态是黑三角在边上 然后其他情况脑补一下,n为偶数先手必胜,可以理解为从某一边取,先手总有办法让后手取得一边有奇数个 #include<iostream> #include< ...

  6. BZOJ 2927: [Poi1999]多边形之战 (博弈)

    题意 有一个凸多边形,顶点编号逆时针从0到n-1.现在这个n边形被剖分成n-2个三角形,给出这n-2个三角形的顶点,保证这是用n-3条不交叉的对角线划分出来的.现在第一个三角形是黑色,其他都是白色.两 ...

  7. bzoj 2935 [Poi1999]原始生物——欧拉回路思路!

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2935 有向图用最小的路径(==总点数最少)覆盖所有边. 完了完了我居然连1999年的题都做不 ...

  8. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  9. 【BZOJ 2802】 2802: [Poi2012]Warehouse Store (贪心)

    2802: [Poi2012]Warehouse Store Description 有一家专卖一种商品的店,考虑连续的n天.第i天上午会进货Ai件商品,中午的时候会有顾客需要购买Bi件商品,可以选择 ...

随机推荐

  1. ast模块

    有这么一个需求,你想从文件中读取字典,方法有很多,这里用的是ast模块 import ast with open("account","r",encoding= ...

  2. JS获取checkbox的个数

    本文算是转载自网络,当时用了他的函数,现在想总结一下,但忘了原文地址了 ================================================================ ...

  3. NoSQL 简介及什么是AICD

    NoSQL 简介 NoSQL(NoSQL = Not Only SQL ),意即"不仅仅是SQL". 在现代的计算系统上每天网络上都会产生庞大的数据量. 这些数据有很大一部分是由关 ...

  4. IOS开发中@2x图片等适应不同分辨率手机

    开发中,例如: nanshanImage.image=[UIImage imageNamed:@'index_pic.png']; 在项目中还保存中index_pic@2x.png的图片,此图为了只适 ...

  5. 在虚机中安装CentOS

    摘要 最近看到.net core 1发布的内容,也想尝试着在lunix上跑一圈.linux这方面的知识一直都没怎么接触过,只在工作中见同事操作过,看到满屏幕的命令行,感觉非常的高大上,趁着现在赶紧学习 ...

  6. Memcached存储命令 - add

    Memcached add 命令用于将 value(数据值) 存储在指定的 key(键) 中. 如果 add 的 key 已经存在,则不会更新数据,之前的值将仍然保持相同,并且您将获得响应 NOT_S ...

  7. pip高级使用技巧以及搭建自己的pypi服务器

    =========================  pip 访问非官方pypi源, 以及代理的设置=========================在Windows下安装某些Python的C ext ...

  8. 微信电脑版微信1.1 for Windows更新 可@人/转发撤回消息/可播小视频

    微信电脑版微信1.1 for Windows发布更新了,版本号为1.1.0.18,群聊可@人/可转发撤回消息/可播小视频,功能越来越接近微信手机版了. 本次更新的一些新特点: 群聊中可以@人. 消息可 ...

  9. hdu.1044.Collect More Jewels(bfs + 状态压缩)

    Collect More Jewels Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Othe ...

  10. 第一节 HTML网页和CSS样式

    1. 第一行 <!DOCTYPE html> 表明网页使用的是HTML5版本 2. 网页的head内容,包含了 title,meta. 3. 网页的标题 title,注意这个显示在浏览器的 ...