NOIP2013 华容道 (棋盘建图+spfa最短路)
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
#define N 35
#define maxn 5000
#define mod 1000000007
#define ll long long
using namespace std; int n,m,q,cte;
int ok[N][N],id[N][N],head[maxn],dis[maxn],inq[maxn],st[];
int xx[]={-,,,};
int yy[]={,,,-};
struct EDGE{
int to,nxt,val;
}edge[maxn*];
struct node{
int x,y;
};
node ins(int s1,int s2) {node k;k.x=s1;k.y=s2;return k;}
bool check(int x,int y)
{
if(x<||y<||x>n||y>m||!ok[x][y]) return false;
else return true;
}
void edge_add(int u,int v,int w)
{
cte++;
edge[cte].to = v;
edge[cte].nxt=head[u];
edge[cte].val= w;
head[u]=cte;
}
void bfs(int ex,int ey,int sx,int sy,int d)
{
queue<node>q;
memset(dis,-,sizeof(dis));
dis[id[ex][ey]]=,dis[id[sx][sy]]=;
q.push(ins(ex,ey));
int x,y,fx,fy;
while(!q.empty())
{
node k=q.front();q.pop();
x=k.x,y=k.y;
for(int i=;i<;i++)
{
fx=x+xx[i],fy=y+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]!=-) continue;
dis[id[fx][fy]]=dis[id[x][y]]+;
q.push(ins(fx,fy));
}
}
for(int i=;i<;i++)
{
fx=sx+xx[i],fy=sy+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]==-||i==d) continue;
edge_add(id[sx][sy]+d,id[sx][sy]+i,dis[id[fx][fy]]);
}
edge_add(id[sx][sy]+d,id[ex][ey]+(d+)%,);
}
void bfsS(int ex,int ey,int sx,int sy)
{
queue<node>q;
memset(dis,-,sizeof(dis));
q.push(ins(ex,ey)),dis[id[ex][ey]]=;
int x,y,fx,fy;
while(!q.empty())
{
node k=q.front();q.pop();
x=k.x,y=k.y;
for(int i=;i<;i++)
{
fx=x+xx[i],fy=y+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]!=-||(fx==sx&&fy==sy)) continue;
dis[id[fx][fy]]=dis[id[x][y]]+;
q.push(ins(fx,fy));
}
}
}
void build_edge()
{
int ct=;
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(ok[i][j]) id[i][j]=ct,ct+=;
}
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(!check(i,j)) continue;
if(check(i-,j)) bfs(i-,j,i,j,);
if(check(i,j+)) bfs(i,j+,i,j,);
if(check(i+,j)) bfs(i+,j,i,j,);
if(check(i,j-)) bfs(i,j-,i,j,);
}
}
void spfa(int S)
{
queue<int>q;
memset(dis,0x3f,sizeof(dis));
q.push(S),dis[S]=,inq[S]=;
while(!q.empty())
{
int u=q.front();q.pop();
for(int j=head[u];j!=-;j=edge[j].nxt)
{
int v=edge[j].to;
if(dis[v]>dis[u]+edge[j].val)
{
dis[v]=dis[u]+edge[j].val;
if(!inq[v]) q.push(v),inq[v]=;
}
}
inq[u]=;
}
} int solve(int ex,int ey,int sx,int sy,int tx,int ty)
{
bfsS(ex,ey,sx,sy);
if(sx==tx&&sy==ty) return ;
for(int i=;i<;i++) st[i]=dis[id[sx+xx[i]][sy+yy[i]]];
int fx,fy,ans=inf;
for(int i=;i<;i++)
{
fx=sx+xx[i],fy=sy+yy[i];
if(!check(fx,fy)||st[i]==-) continue;
spfa(id[sx][sy]+i);
for(int j=;j<;j++)
{
if(!check(tx+xx[j],ty+yy[j])||dis[id[tx][ty]+j]==-) continue;
ans=min(ans,st[i]+dis[id[tx][ty]+j]);
}
}
if(ans>=inf) return -;
else return ans;
} int main()
{
freopen("aa.in","r",stdin);
scanf("%d%d%d",&n,&m,&q);
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
scanf("%d",&ok[i][j]);
memset(head,-,sizeof(head));
build_edge();
int ex,ey,sx,sy,tx,ty;
while(q--)
{
scanf("%d%d%d%d%d%d",&ex,&ey,&sx,&sy,&tx,&ty);
printf("%d\n",solve(ex,ey,sx,sy,tx,ty));
}
return ;
}
---恢复内容开始---
题目大意:类似于华容道游戏,但所有棋子都是11大小的,棋盘上仅有一个空格,其它位置填满了棋子,有些棋子是固定的不能被移动,每个格子内只能有一个棋子,棋子只能移动到相邻的空格上,每移动任意棋子一次视为一次操作,求让指定棋子移动到指定位置的最小操作数
据说爆搜能拿70pts....
观察爆搜为什么会T,主要是因为搜出了许多无用的状态
考虑优化这个过程,空格在指定棋子的上下左右分别记为一种状态,相当于一个位置记录了4个状态
为了方便转移,直接把每个格子的四个状态记成连续的,这样方便转移
每种状态,它能转移到其它状态的情况分为两种:
<1>指定格子的某个方向有空格,转移到同一个格子的另一个方向有空格,指定格子上的棋子的位置不动
这种情况,我们对于每个格子的四个方向分别bfs即可,因为除了空格子,其它位置都被填满了,所以空格子在被填满的棋盘里移动,和棋子在空棋盘里的移动方式是一样的(貌似是一句废话但我仍然思考了很久...)
<2>指定棋子朝着这个状态指向的空格移动
这种情况,显然操作次数为1,转移的状态呢,则是从 当前格子的指向 转移到 目标格子的反方向 (就是那个(d+2)%4)
虽然有500次询问,但图都是一样的,所以建一次图就行
那么,如果一开始空格子不在指定棋子边上呢?
我们把空格子移动到指定棋子的四个方向就行了,注意不要让空格子移动到指定棋子的格子上
再分别跑一次SPFA最短路即可
总结:之前还有一道棋盘建图题(荷叶塘),但那道题建图更简单更好想,思路都是先bfs建边,然后跑最短路
注意特判移动次数为0的情况!
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
#define N 35
#define maxn 5000
#define mod 1000000007
#define ll long long
using namespace std; int n,m,q,cte;
int ok[N][N],id[N][N],head[maxn],dis[maxn],inq[maxn],st[];
int xx[]={-,,,};
int yy[]={,,,-};
struct EDGE{
int to,nxt,val;
}edge[maxn*];
struct node{
int x,y;
};
node ins(int s1,int s2) {node k;k.x=s1;k.y=s2;return k;}
bool check(int x,int y)
{
if(x<||y<||x>n||y>m||!ok[x][y]) return false;
else return true;
}
void edge_add(int u,int v,int w)
{
cte++;
edge[cte].to = v;
edge[cte].nxt=head[u];
edge[cte].val= w;
head[u]=cte;
}
void bfs(int ex,int ey,int sx,int sy,int d)
{
queue<node>q;
memset(dis,-,sizeof(dis));
dis[id[ex][ey]]=,dis[id[sx][sy]]=;
q.push(ins(ex,ey));
int x,y,fx,fy;
while(!q.empty())
{
node k=q.front();q.pop();
x=k.x,y=k.y;
for(int i=;i<;i++)
{
fx=x+xx[i],fy=y+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]!=-) continue;
dis[id[fx][fy]]=dis[id[x][y]]+;
q.push(ins(fx,fy));
}
}
for(int i=;i<;i++)
{
fx=sx+xx[i],fy=sy+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]==-||i==d) continue;
edge_add(id[sx][sy]+d,id[sx][sy]+i,dis[id[fx][fy]]);
}
edge_add(id[sx][sy]+d,id[ex][ey]+(d+)%,);
}
void bfsS(int ex,int ey,int sx,int sy)
{
queue<node>q;
memset(dis,-,sizeof(dis));
q.push(ins(ex,ey)),dis[id[ex][ey]]=;
int x,y,fx,fy;
while(!q.empty())
{
node k=q.front();q.pop();
x=k.x,y=k.y;
for(int i=;i<;i++)
{
fx=x+xx[i],fy=y+yy[i];
if(!check(fx,fy)||dis[id[fx][fy]]!=-||(fx==sx&&fy==sy)) continue;
dis[id[fx][fy]]=dis[id[x][y]]+;
q.push(ins(fx,fy));
}
}
}
void build_edge()
{
int ct=;
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(ok[i][j]) id[i][j]=ct,ct+=;
}
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(!check(i,j)) continue;
if(check(i-,j)) bfs(i-,j,i,j,);
if(check(i,j+)) bfs(i,j+,i,j,);
if(check(i+,j)) bfs(i+,j,i,j,);
if(check(i,j-)) bfs(i,j-,i,j,);
}
}
void spfa(int S)
{
queue<int>q;
memset(dis,0x3f,sizeof(dis));
q.push(S),dis[S]=,inq[S]=;
while(!q.empty())
{
int u=q.front();q.pop();
for(int j=head[u];j!=-;j=edge[j].nxt)
{
int v=edge[j].to;
if(dis[v]>dis[u]+edge[j].val)
{
dis[v]=dis[u]+edge[j].val;
if(!inq[v]) q.push(v),inq[v]=;
}
}
inq[u]=;
}
} int solve(int ex,int ey,int sx,int sy,int tx,int ty)
{
bfsS(ex,ey,sx,sy);
if(sx==tx&&sy==ty) return ;
for(int i=;i<;i++) st[i]=dis[id[sx+xx[i]][sy+yy[i]]];
int fx,fy,ans=inf;
for(int i=;i<;i++)
{
fx=sx+xx[i],fy=sy+yy[i];
if(!check(fx,fy)||st[i]==-) continue;
spfa(id[sx][sy]+i);
for(int j=;j<;j++)
{
if(!check(tx+xx[j],ty+yy[j])||dis[id[tx][ty]+j]==-) continue;
ans=min(ans,st[i]+dis[id[tx][ty]+j]);
}
}
if(ans>=inf) return -;
else return ans;
} int main()
{
freopen("aa.in","r",stdin);
scanf("%d%d%d",&n,&m,&q);
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
scanf("%d",&ok[i][j]);
memset(head,-,sizeof(head));
build_edge();
int ex,ey,sx,sy,tx,ty;
while(q--)
{
scanf("%d%d%d%d%d%d",&ex,&ey,&sx,&sy,&tx,&ty);
printf("%d\n",solve(ex,ey,sx,sy,tx,ty));
}
return ;
}
---恢复内容结束---
NOIP2013 华容道 (棋盘建图+spfa最短路)的更多相关文章
- POJ 1556 - The Doors - [平面几何+建图spfa最短路]
题目链接:http://poj.org/problem?id=1556 Time Limit: 1000MS Memory Limit: 10000K Description You are to f ...
- poj 2135 Farm Tour 最小费用最大流建图跑最短路
题目链接 题意:无向图有N(N <= 1000)个节点,M(M <= 10000)条边:从节点1走到节点N再从N走回来,图中不能走同一条边,且图中可能出现重边,问最短距离之和为多少? 思路 ...
- Invitation Cards(邻接表+逆向建图+SPFA)
Time Limit: 8000MS Memory Limit: 262144K Total Submissions: 17538 Accepted: 5721 Description In ...
- hdu 4885 (n^2*log(n)推断三点共线建图)+最短路
题意:车从起点出发,每次仅仅能行驶L长度,必需加油到满,每次仅仅能去加油站或目的地方向,路过加油站就必需进去加油,问最小要路过几次加油站. 開始时候直接建图,在范围内就有边1.跑最短了,再读题后发现, ...
- hdu 4885 (n^2*log(n)判断三点共线建图)+最短路
题意:车从起点出发,每次只能行驶L长度,必需加油到满,每次只能去加油站或目的地方向,路过加油站就必需进去加油,问最小要路过几次加油站. 开始时候直接建图,在范围内就有边1.跑最短了,再读题后发现,若几 ...
- CF786B Legacy 线段树优化建图 + spfa
CodeForces 786B Rick和他的同事们做出了一种新的带放射性的婴儿食品(???根据图片和原文的确如此...),与此同时很多坏人正追赶着他们.因此Rick想在坏人们捉到他之前把他的遗产留给 ...
- B - Legacy CodeForces - 787D 线段树优化建图+dij最短路 基本套路
B - Legacy CodeForces - 787D 这个题目开始看过去还是很简单的,就是一个最短路,但是这个最短路的建图没有那么简单,因为直接的普通建图边太多了,肯定会超时的,所以要用线段树来优 ...
- HDU1281(二分图最大匹配,棋盘建图,找关键点)
棋盘游戏 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submis ...
- 【[NOI2011]智能车比赛】(建图+spfa+坑爹精度)
过了这题我就想说一声艹,跟这个题死磕了将近6个小时,终于是把这个题死磕出来了.首先看到这个题的第一反应,和当初做过的一个房间最短路比较相似,然后考虑像那个题那样建边,然后跑最短路.(具体建边方法请参考 ...
随机推荐
- vue 配置页面动态的 title
- POJ 1811 Prime Test( Pollard-rho整数分解经典题 )
链接:传送门 题意:输入 n ,判断 n 是否为素数,如果是合数输出 n 的最素因子 思路:Pollard-rho经典题 /************************************** ...
- 实现路由器自动登录校园网(edu)
准备工作: (1)一个可以刷openwrt固件的路由器,如大多人使用的crazybox版本的路由. (2)一个可用的edu账号. (3)一个浏览器(firfox,chrome) 下面开始: 一:刷op ...
- Jquery 根据HTML内容选择元素
选择所有包含 "is" 的 元素: $("p:contains(is)")
- N - 畅通工程再续
N - 畅通工程再续 思路:zz #include<cmath> #include<cstdio> #include<cstring> #include<io ...
- bzoj 1266 1266: [AHOI2006]上学路线route
1266: [AHOI2006]上学路线route Time Limit: 3 Sec Memory Limit: 162 MBSubmit: 2356 Solved: 841[Submit][S ...
- POJ 2133
类似于DP一样做,但这题有个大坑,自己看DIS吧.... #include <iostream> #include <cstdio> #include <cstring& ...
- HDU 3886
一开始又往打表想了....不过,打表确实不好处理,转DFS. DFS有几个要注意的问题,1.对于枚举以零开始的数.我纠结了很久,最终学习别人的方法,设一个BOOL,并且假设最高一位有零,很方便.2.当 ...
- C#高级编程八十三天----程序集的含义
程序集的含义 一.程序集是包括一个或多个类型定义文件和资源文件的集合.它同意我们分析可重用类型的逻辑表示和物理表示. 相当于你定义了一个项目XXProject,项目存在非常多文件(类,窗口,接口,资源 ...
- js 数组克隆
在JS中,对于对象类型的变量的引用是通过类似C之类的"指针"的方式来操作的,及假设多个变量引用同一个对象,则随意一个变量对对象的操作都会影响到其他的变量. 这个可能不是非常明确,以 ...