【NOIP2013】华容道

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#define RG register
#define IL inline
#define pi acos(-1.0)
#define ll long long 
using namespace std;

int gi() {
  char ch=getchar(); int x=0;
  while(ch<'0' || ch>'9') ch=getchar();
  while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}
  return x;
}

const int maxn = 33;

int n,m,que,inf,ans,id_num;
int ex,ey,sx,sy,tx,ty;
int g[maxn][maxn],id[maxn][maxn][6],dis[5050],mov[maxn][maxn][6][6];
int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1};
bool vis[maxn][maxn],in[5050];
int nxt[5050],to[5050],w[5050],h[5050],e_num;
struct NODE {
  int x,y,t;
};

queue<NODE> q1,q2;
queue<int> q;

void add(int x, int y, int v) {
  nxt[++e_num]=h[x];
  to[e_num]=y;
  w[e_num]=v;
  h[x]=e_num;
}

void bfs1(int X, int Y, int k, int l) {
  while(!q1.empty()) q1.pop();
  memset(vis,0,sizeof(vis));
  int kx=X+dx[k],ky=Y+dy[k];
  int lx=X+dx[l],ly=Y+dy[l];
  if(kx==lx && ky==ly) {mov[X][Y][k][l]=1;return;}
  NODE now;
  now.x=kx,now.y=ky,now.t=0;
  q1.push(now);
  while(!q1.empty()) {
    NODE u=q1.front();
    q1.pop();
    for(int p=1; p<=4; p++) {
      NODE v;
      v.x=u.x+dx[p],v.y=u.y+dy[p],v.t=u.t+1;
      if(v.x==X && v.y==Y) continue;
      if(!vis[v.x][v.y] && g[v.x][v.y]) {
	if(v.x==lx && v.y==ly) {mov[X][Y][k][l]=v.t+1;return;}//移动+1
	vis[v.x][v.y]=1;
	q1.push(v);
      }
    }
  }
}

void bfs2() {
  while(!q2.empty()) q2.pop();
  memset(vis,0,sizeof(vis));
  NODE now;
  now.x=ex,now.y=ey,now.t=0;
  q2.push(now);
  vis[now.x][now.y]=1;
  while(!q2.empty()) {
    NODE u=q2.front();
    q2.pop();
    for(int p=1; p<=4; p++) {
      NODE v;
      v.x=u.x+dx[p],v.y=u.y+dy[p],v.t=u.t+1;
      if(!g[v.x][v.y] || vis[v.x][v.y]) continue;
      if(v.x==sx && v.y==sy) {
	if(p==1) dis[id[sx][sy][2]]=min(dis[id[sx][sy][2]],u.t);
	if(p==2) dis[id[sx][sy][1]]=min(dis[id[sx][sy][1]],u.t);
	if(p==3) dis[id[sx][sy][4]]=min(dis[id[sx][sy][4]],u.t);
	if(p==4) dis[id[sx][sy][3]]=min(dis[id[sx][sy][3]],u.t);
	continue;
      }
      vis[v.x][v.y]=1;
      q2.push(v);
    }
  }
}

void spfa() {
  memset(dis,127/3,sizeof(dis));
  bfs2();//将空格移动到目标块的四个方向,并以此为dis的初值进行spfa
  //for(int i=1; i<=4; i++) cout<<dis[id[sx][sy][i]]<<endl;
  q.push(id[sx][sy][1]),in[id[sx][sy][1]]=1;
  q.push(id[sx][sy][2]),in[id[sx][sy][2]]=1;
  q.push(id[sx][sy][3]),in[id[sx][sy][3]]=1;
  q.push(id[sx][sy][4]),in[id[sx][sy][4]]=1;
  while(!q.empty()) {
    int u=q.front();
    in[u]=0,q.pop();
    for(int i=h[u]; i; i=nxt[i]) {
      int v=to[i];
      if(dis[u]+w[i]<dis[v]) {
	dis[v]=dis[u]+w[i];
	if(!in[v]) in[v]=1,q.push(v);
      }
    }
  }
  for(int i=1; i<=4; i++)
    ans=min(ans,dis[id[tx][ty][i]]);
}

int main() {
  freopen("华容道.in","r",stdin);
  freopen("华容道.out","w",stdout);
  n=gi(),m=gi(),que=gi();
  for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++)
      g[i][j]=gi();
  memset(mov,127/3,sizeof(mov));
  inf=mov[0][0][0][0];
  for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++) {
      if(!g[i][j]) continue;
      for(int k=1; k<=4; k++) {
	int kx=i+dx[k],ky=j+dy[k];
	if(g[kx][ky])
	id[i][j][k]=++id_num;
      }
    }//id[i][j][k]表示点(i,j),空格在k方向
  for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++) {
      if(!g[i][j]) continue;
      for(int k=1; k<=4; k++)
	for(int l=1; l<=4; l++) {
	  int kx=i+dx[k],ky=j+dy[k];
	  int lx=i+dx[l],ly=j+dy[l];
	  if(g[kx][ky] && g[lx][ly]) bfs1(i,j,k,l);
	}
    }//mov[i][j][k][l]表示点(i,j),空格在k方向,并往l方向移动一格的步数
  for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++) {
      if(!g[i][j]) continue;
      for(int k=1; k<=4; k++)
	for(int l=1; l<=4; l++) {
	  int kx=i+dx[k],ky=j+dy[k];
	  int lx=i+dx[l],ly=j+dy[l];
	  if(!g[kx][ky] || !g[lx][ly]) continue;
	  if(l==1) add(id[i][j][k],id[lx][ly][2],mov[i][j][k][l]);//上移
	  if(l==2) add(id[i][j][k],id[lx][ly][1],mov[i][j][k][l]);//下
	  if(l==3) add(id[i][j][k],id[lx][ly][4],mov[i][j][k][l]);//右
	  if(l==4) add(id[i][j][k],id[lx][ly][3],mov[i][j][k][l]);//左
	}
    }
  //cout<<mov[1][2][2][2]<<endl;
  while(que--) {
    ex=gi(),ey=gi(),sx=gi(),sy=gi(),tx=gi(),ty=gi();
    if(sx==tx && sy==ty) {puts("0");continue;}
    ans=inf;
    spfa();
    if(ans==inf) puts("-1");
    else printf("%d\n", ans);
  }
  return 0;
}
/*
  puts自带换行
  特判的点没有直接continue
*/


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值