https://www.zybuluo.com/ysner/note/1293166

题面

戳我

解析

挺不错的一道图论码量题。

可以借此回顾一下\(noip2013\)华容道。

思路和华容道差不多。

照洛谷数据规模看,暴搜可能有\(\frac{1}{3}\)的分数。。。

设\(f[w][i][j]\)表示箱子在点\((i,j)\)时,人能否到达\(w\)方向(即\(0/1/2/3\),表示上下左右)。

那么我们可以先从人出发\(BFS\)一遍,预处理出未推箱子时的\(f[w][i][j]\)。

注意箱子所在点是可以多次更新的。

然后开始推箱子,每次有两种决策:

  • 换个方向
  • 照原方向推一格

决策可不可行肯定要预处理出来,毕竟\(n\)是\(10^6\)级别的。

显然如果能换方向,就说明换方向的前后两点在同一点双联通分量中。

这玩意儿从人的出发点跑一遍\(Tarjan\)就能预处理出来。

推一格这个操作可以直接转移。

注意边界都要赋为\(\#\)。(尤其是每行末尾)

细节挺多,蒟蒻调了两个多小时。。。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
#define ll long long
#define re register
#define il inline
#define pb(a) push_back(a)
#define gc() getchar()
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int N=1644;
int n,m,q,sx,sy,rx,ry,dfn[N*N],low[N*N],top,sta[N*N],tot,Id[N][N];
int mx[4]={1,-1,0,0},my[4]={0,0,1,-1};
bool f[4][N][N],vis[N][N];
char mp[N][N];
struct dat{int x,y,w;};
queue<dat>Q;
vector<int>c[N*N];
il int id(re int x,re int y){return Id[x][y];}
il int gi()
{
re int x=0,t=1;
re char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=gc();
if(ch=='-') t=-1,ch=gc();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=gc();
return x*t;
}
il void Tarjan(re int x,re int y,re int fx,re int fy)
{
re int u=id(x,y);
sta[++top]=u,dfn[u]=low[u]=++tot;
fp(i,0,3)
{
re int xx=x+mx[i],yy=y+my[i],v=id(xx,yy);
if(mp[xx][yy]=='#') continue;
if(xx==fx&&yy==fy) continue;
if(!dfn[v])
{
Tarjan(xx,yy,x,y);
low[u]=min(low[u],low[v]);
if(low[v]>=dfn[u])
{
++tot;re int g;
do{g=sta[top--];c[g].pb(tot);}while(g^v);
c[u].pb(tot);
}
}
else low[u]=min(low[u],dfn[v]);
}
}
il void BFS1()
{
Q.push((dat){rx,ry,0});vis[rx][ry]=1;
while(!Q.empty())
{
re int x=Q.front().x,y=Q.front().y,u=id(x,y);Q.pop();
fp(i,0,3)
{
re int xx=x+mx[i],yy=y+my[i],v=id(xx,yy);
if(mp[xx][yy]=='#'||vis[xx][yy]) continue;
if(xx==sx&&yy==sy) {f[i^1][xx][yy]=1;continue;}
vis[xx][yy]=1;Q.push((dat){xx,yy,0});
}
}
}
il int check(re int x,re int y,re int xx,re int yy)
{
re int u=id(x,y),v=id(xx,yy);
for(re int i=0;i<c[u].size();++i)
for(re int j=0;j<c[v].size();++j)
if(c[u][i]==c[v][j]) return 1;
return 0;
}
il void BFS2()
{
fp(i,0,3) if(f[i][sx][sy]) Q.push((dat){sx,sy,i});
while(!Q.empty())
{
re int x=Q.front().x,y=Q.front().y,w=Q.front().w,u=id(x,y);Q.pop();
fp(i,0,3)
{
if(f[i][x][y]) continue;
re int xx=x+mx[i],yy=y+my[i];
if(mp[xx][yy]=='#') continue;
if(check(x+mx[w],y+my[w],xx,yy)) f[i][x][y]=1,Q.push((dat){x,y,i});
}
{
re int xx=x+mx[w^1],yy=y+my[w^1];///yy=y+mx???????
if(mp[xx][yy]=='#') continue;
if(!f[w][xx][yy]) f[w][xx][yy]=1,Q.push((dat){xx,yy,w});
}
}
}
int main()
{
n=gi();m=gi();q=gi();
memset(mp,'#',sizeof(mp));
fp(i,1,n) scanf("%s",mp[i]+1),mp[i][m+1]='#';
fp(i,1,n)
fp(j,1,m)
{
if(mp[i][j]=='A') rx=i,ry=j,mp[i][j]='.';
else if(mp[i][j]=='B') sx=i,sy=j,mp[i][j]='.';
Id[i][j]=(i-1)*m+j;
}
Tarjan(rx,ry,0,0);
BFS1();
BFS2();
while(q--)
{
re int x=gi(),y=gi(),tag=(sx==x&&sy==y);
fp(i,0,3) tag|=f[i][x][y];
puts(tag?"YES":"NO");
}
return 0;
}

最新文章

  1. 打印Lua的Table对象
  2. delphi7 在虚拟机 vbox里面安装失败
  3. SQL Server选取本周或上一周数据
  4. Docker探索系列1之docker入门安装与操作
  5. HDU-敌兵布阵
  6. 微软TechEd2013大会将在北京、上海召开!
  7. lintcode :Remove Duplicates from Sorted List 删除排序链表中的重复元素
  8. Qt5 文本编辑
  9. linux下的进程、网络、性能监控命令
  10. 使用javassist运行时动态重新加载java类及其他替换选择
  11. eclipse项目导入到android studio中文乱码处理
  12. 保证Android后台不被杀死的几种方法
  13. UVA1374 IDA*
  14. shell for 循环
  15. eureka 服务注册与发现
  16. Python多线程应用示例
  17. 执行次数最多的sql语句
  18. JsonPath教程
  19. Excel批量导入商品,遇到导入失败记录到另一个Excel中供下载查看
  20. 第5天 Java基础语法

热门文章

  1. HTML、CSS常用技巧
  2. 三菱PLC FB库函数调用方法 (Gx Work2版本)
  3. hadoop_exporter
  4. Thawte 企业版代码签名证书
  5. Fiddler抓取https相关设置
  6. Linux笔记:定时任务和文件操作
  7. POJ1422 Air Raid
  8. php的socket通信【转载】
  9. POJ 3268_Silver Cow Party
  10. Bad Hair Day-POJ3250(简单的入栈出栈)