题解 \(by\;zj\varphi\)

NOIP2013 的原题

最简单的思路就是一个 bfs,可以拿到 \(70pts\)

75pts
#include<bits/stdc++.h>
#define ri signed
#define pd(i) ++i
#define bq(i) --i
#define func(x) std::function<x>
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
#define debug1(x) std::cerr << #x"=" << x << ' '
#define debug2(x) std::cerr << #x"=" << x << std::endl
#define Debug(x) assert(x)
struct nanfeng_stream{
template<typename T>inline nanfeng_stream &operator>>(T &x) {
bool f=false;x=0;char ch=gc();
while(!isdigit(ch)) f|=ch=='-',ch=gc();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=gc();
return x=f?-x:x,*this;
}
}cin;
}
using IO::cin;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
static const int N=31;
static const int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
int n,m,q,ex,ey,sx,sy,tx,ty;
bool vis[N][N][N][N],est[N][N];
struct Que{int x,y,nx,ny,stp;}que[N*N*N*N];
auto bfs=[]() {
int hd=1,tl=0;
que[tl=1].x=ex,que[1].y=ey;
que[1].nx=sx,que[1].ny=sy;
while(hd<=tl) {
int x=que[hd].x,y=que[hd].y;
for (ri i(0);i<4;pd(i)) {
int stx=x+dx[i],sty=y+dy[i],nx=que[hd].nx,ny=que[hd].ny;
if (stx==nx&&sty==ny) nx=x,ny=y;
if (stx<1||stx>n||sty<1||sty>m||!est[stx][sty]||vis[stx][sty][nx][ny]) continue;
++tl;
que[tl].x=stx,que[tl].y=sty,que[tl].nx=nx,que[tl].ny=ny,que[tl].stp=que[hd].stp+1;
vis[stx][sty][nx][ny]=true;
if (nx==tx&&ny==ty) return que[tl].stp;
}
++hd;
}
return -1;
};
inline int main() {
// FI=freopen("y.in","r",stdin);
// FO=freopen("y.out","w",stdout);
cin >> n >> m >> q;
for (ri i(1);i<=n;pd(i))
for (ri j(1);j<=m;pd(j)) cin >> est[i][j];
for (ri z(1);z<=q;pd(z)) {
cin >> ex >> ey >> sx >> sy >> tx >> ty;
if (sx==tx&&sy==ty) printf("0\n");
else {
memset(vis,false,sizeof(vis));
printf("%d\n",bfs());
}
}
return 0;
}
}
int main() {return nanfeng::main();}

常数写小一点,再开个 \(O2\) 洛谷和 loj 就过了。

正解:

对于一个点,它能移动只有当它周围四个方向存在空格子时,所以一个思路就是先将空格子移动到初始格子周围,再向目标点移动,同时维护空格点的位置。

设 \(move_{x,y,t1,t2}\) 表示当前在 \((x,y)\) 的格子,由空格子在它的左/上侧转移到右/下侧需要多少步,这个可以 bfs 预处理出来。

然后在询问时再用一个 bfs 处理出空格子到初始点的四个状态中所需要的步数。

设 \(f_{x,y,t}\) 表示点 \((x,y)\),空格在它的 \(t\) 状态,由初始点转移过来要多少步。

转移就是 \(f_{x,y,t}=\min (f_{xx,yy,tt}+mv_{xx,yy,tt,i}+1)\)

其中 \(t\) 是 \(i\) 对称过去的状态,\((x,y)\) 是 \((xx,yy)\) 周围的点。

最后跑一个最短路,特判一下无解。

Code
#include<bits/stdc++.h>
#define ri signed
#define pd(i) ++i
#define bq(i) --i
#define func(x) std::function<x>
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
#define debug1(x) std::cerr << #x"=" << x << ' '
#define debug2(x) std::cerr << #x"=" << x << std::endl
#define Debug(x) assert(x)
struct nanfeng_stream{
template<typename T>inline nanfeng_stream &operator>>(T &x) {
bool f=false;x=0;char ch=gc();
while(!isdigit(ch)) f|=ch=='-',ch=gc();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=gc();
return x=f?-x:x,*this;
}
}cin;
}
using IO::cin;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
static const int N=31;
static const int dx[]={0,1,0,-1},dy[]={1,0,-1,0};
int mv[N][N][4][4],dis[N][N][4],n,m,q,ex,ey,sx,sy,tx,ty,cnt;
bool vis[N][N],vs[N][N][4],est[N][N];
struct Que{int x,y,d;}qe[N*N*N*N];
auto bfs=[](int x,int y,int f1,int f2) {
int sx=x+dx[f1],sy=y+dy[f1],tx=x+dx[f2],ty=y+dy[f2];
if ((!est[sx][sy])||(!est[tx][ty])) return;
std::queue<Que> que;
Que tmp={sx,sy,0};
que.push(tmp);
memset(vis,false,sizeof(vis));
vis[x][y]=vis[sx][sy]=true;
while(!que.empty()) {
tmp=que.front(),que.pop();
for (ri i(0);i<4;pd(i)) {
int kx=tmp.x+dx[i],ky=tmp.y+dy[i];
if (kx<1||kx>n||ky<1||ky>m) continue;
if (!vis[kx][ky]&&est[kx][ky]) {
vis[kx][ky]=true;
Que nw={kx,ky,tmp.d+1};
if (kx==tx&&ky==ty) {mv[x][y][f1][f2]=nw.d;return;}
que.push(nw);
}
}
}
};
auto init=[]() {
memset(mv,0x3f,sizeof(mv));
for (ri i(1);i<=n;pd(i))
for (ri j(1);j<=m;pd(j))
if (est[i][j])
for (ri k(0);k<4;pd(k))
for (ri l(0);l<4;pd(l))
if (l==k) mv[i][j][k][l]=0;
else bfs(i,j,k,l);
};
auto solve=[](int x,int y) {
cnt=0;
memset(vis,false,sizeof(vis));
vis[x][y]=vis[sx][sy]=true;
Que tmp={x,y,0};
std::queue<Que> que;
que.push(tmp);
while(!que.empty()) {
tmp=que.front(),que.pop();
for (ri i(0);i<4;pd(i)) {
int kx=tmp.x+dx[i],ky=tmp.y+dy[i];
if (kx<1||kx>n||ky<1||ky>m) continue;
if (est[kx][ky]&&!vis[kx][ky]) {
vis[kx][ky]=true;
que.push({kx,ky,tmp.d+1});
} else if (kx==sx&&ky==sy) qe[++cnt]={kx,ky,(i+2)%4},dis[kx][ky][(i+2)%4]=tmp.d;
}
}
};
auto spfa=[]() {
memset(vs,false,sizeof(vs));
std::queue<Que> que;
for (ri i(1);i<=cnt;pd(i)) que.push(qe[i]),vs[qe[i].x][qe[i].y][qe[i].d]=true;
Que tmp;
while(!que.empty()) {
tmp=que.front(),que.pop();
vs[tmp.x][tmp.y][tmp.d]=false;
for (ri i(0);i<4;pd(i)) {
int kx=tmp.x+dx[i],ky=tmp.y+dy[i],td;
if (kx<1||kx>n||ky<1||ky>m) continue;
if (dis[kx][ky][td=(i+2)%4]>dis[tmp.x][tmp.y][tmp.d]+mv[tmp.x][tmp.y][tmp.d][i]+1) {
dis[kx][ky][td]=dis[tmp.x][tmp.y][tmp.d]+mv[tmp.x][tmp.y][tmp.d][i]+1;
if (!vs[kx][ky][td]) {
vs[kx][ky][td]=true;
que.push({kx,ky,td});
}
}
}
}
};
inline int main() {
FI=freopen("y.in","r",stdin);
FO=freopen("y.out","w",stdout);
cin >> n >> m >> q;
for (ri i(1);i<=n;pd(i))
for (ri j(1);j<=m;pd(j)) cin >> est[i][j];
init();
for (ri z(1);z<=q;pd(z)) {
cin >> ex >> ey >> sx >> sy >> tx >> ty;
if (sx==tx&&sy==ty) printf("0\n");
else {
memset(dis,0x3f,sizeof(dis));
solve(ex,ey);
spfa();
int ans=1061109567;
for (ri i(0);i<4;pd(i)) ans=cmin(ans,dis[tx][ty][i]);
printf("%d\n",ans==1061109567?-1:ans);
}
}
return 0;
}
}
int main() {return nanfeng::main();}

NOIP 模拟 $79\; \rm y$的更多相关文章

  1. NOIP 模拟 $20\; \rm y$

    题解 \(by\;zj\varphi\) 首先发现一共最多只有 \(2^d\) 种道路,那么可以状压,(不要 \(dfs\),会搜索过多无用的状态) 那么设 \(f_{i,j,k}\) 为走 \(i\ ...

  2. NOIP模拟3

    期望得分:30+90+100=220 实际得分:30+0+10=40 T1智障错误:n*m是n行m列,硬是做成了m行n列 T2智障错误:读入三个数写了两个%d T3智障错误:数值相同不代表是同一个数 ...

  3. 7.22 NOIP模拟7

    又是炸掉的一次考试 T1.方程的解 本次考试最容易骗分的一道题,但是由于T2花的时间太多,我竟然连a+b=c都没判..暴力掉了40分. 首先a+b=c,只有一组解. 然后是a=1,b=1,答案是c-1 ...

  4. NOIP模拟 1

    NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. #   用  户  名   ...

  5. 20190725 NOIP模拟8

    今天起来就是虚的一批,然后7.15开始考试,整个前半个小时异常的困,然后一看题,T1一眼就看出了是KMP,但是完了,自己KMP的打法忘的一干二净,然后开始打T2,T2肝了一个tarjan点双就扔上去了 ...

  6. 20190902+0903合集-NOIP模拟

    一直没时间写QwQ 于是补一下. Day 1 晚饭吃的有点恶心…… $1s\,2s\,5s$ 还开 -O2 ?? 有点恐怖. T1 猛的一想: 把外面设成一个点, 向入口连一条权为排队时间的边 从出口 ...

  7. 2021.5.22 noip模拟1

    这场考试考得很烂 连暴力都没打好 只拿了25分,,,,,,,,好好总结 T1序列 A. 序列 题目描述 HZ每周一都要举行升旗仪式,国旗班会站成一整列整齐的向前行进. 郭神作为摄像师想要选取其中一段照 ...

  8. 5.23考试总结(NOIP模拟2)

    5.23考试总结(NOIP模拟2) 洛谷题单 看第一题第一眼,不好打呀;看第一题样例又一眼,诶,我直接一手小阶乘走人 然后就急忙去干T2T3了 后来考完一看,只有\(T1\)骗到了\(15pts\)[ ...

  9. 5.22考试总结(NOIP模拟1)

    5.22考试总结(NOIP模拟1) 改题记录 T1 序列 题解 暴力思路很好想,分数也很好想\(QAQ\) (反正我只拿了5pts) 正解的话: 先用欧拉筛把1-n的素数筛出来 void get_Pr ...

随机推荐

  1. K8s工作流程详解

    在学习k8s工作流程之前,我们得再次认识一下上篇k8s架构与组件详解中提到的kube-controller-manager一个k8s中许多控制器的进程的集合. 比如Deployment 控制器(Dep ...

  2. git各种操作:基本操作 and 多人协作 and 冲突解决

    git基本操作 git 上次文件到远程仓库(参考:https://blog.csdn.net/beiqiaofeng123/article/details/104859326) 如果第一次上传,配置一 ...

  3. 解决umount: /home: device is busy

    取消挂载/home时出现umount: /home: device is busy,        原因是因为有程序在使用/home目录,我们可以使用fuser查看那些程序的进程,        然后 ...

  4. 注释swap分区

    grep  "#" fstab |grep "swap" >/dev/null || sed -i 's/^.*swap/#&/g' fstab

  5. 学习了解PHP中的SeasLog日志扩展

    今天来学习的扩展是和日志相关的一个扩展,对于 PHP 的日志应用来说,除了本身自带的 error_log() . syslog() 之外,在大多数的框架中还会经常见到 monolog 的踪影.当然,我 ...

  6. mysql 选取操作日志(mysql-bin.0000恢复数据)

    my.ini 配置log-bin=mysql-bin 启用日志 用  mysql-bin.0000 mysqlbinlog -d keyboard ../data/mysql-bin.000024 – ...

  7. 网页兼容最新IE声明meta方法

    第三种,总是使用最新版本文档模式. 以下是例子: <meta http-equiv="X-UA-Compatible" content="IE=edge" ...

  8. 3gcms-Flash幻灯片上传后图片模糊解决办法

    很简单,不用纠结,直接修改admin/lib/action/FileAction.class.php 将 $upload->thumbMaxWidth='300'; //以字串格式来传,如果你希 ...

  9. Xshell和Xftp - 下载安装

    简介 Xshell 实际工作运用:连接Linux Xftp 实际工作运用:传输文件到Linux系统 下载安装 三连后评论区留言私发,此贴长期有效!!!

  10. typora+picgo+gitee 图片无法加载问题

    typora+picgo+gitee 图片无法加载问题 待解决