【JZOJ5781】秘密通道

Description

有一副n*m的地图，有n*m块地，每块是下列四种中的一种：
墙：用#表示，墙有4个面，分别是前面，后面，左面，右面。
起点：用C表示，为主角的起点，是一片空地。
终点：用F表示，为主角的目的地，是一片空地。
空地：用 . 表示。
其中除了墙不能穿过，其他地方都能走。
 
主角有以下3种操作：
1.移动到相邻的前后左右的地方，花费一个单位时间。
2.向前后左右其中一个方向发射子弹，子弹沿直线穿过，打在最近的一堵墙的一面，然后墙的这面就会形成一个开口通往秘密通道。同一时间最多只能有两个开口，若出现有3个开口，出现时间最早的开口会立即消失。该操作不用时间。
3.可以从一个与开口相邻的空地跳进去，进入秘密通道，从另外一个开口正对的空地跳出来。这个过程花费一个单位时间。

地图四周都是墙，问主角最少用多少时间从C走到F。C和F
只会出现一次。
 

Input

第一行输入两个正整数n，m。
接下来n行，每行m个字符描述地图。
 

Output

输出1个整数，表示最短时间完成路途。如果无解输出nemoguce
 

Sample Input

Input 1
4 4
####
#.F#
#C.#
####
Input 2
6 8
########
#.##..F#
#C.##..#
#..#...#
#.....##
########
Input 3
4 5
#####
#C#.#
###F#
#####
 

Sample Output

Output 1
2
Output 2
4
Output 3
nemoguce

Data Constraint

对于50%的数据，4≤ n,m≤ 15。
对于100%的数据，4≤ n,m≤ 500。
 

Hint

总共用到8次操作，时间之和为4。如下图所示
1.向左射一枪，在(3,1)的右面出现开口。
2.向下射一枪，在(6,2)的上面出现开口。
3.向左从(3,1)进入秘密通道，从(6,2)中出来，到达(5,2)。用1单位时间。
4.向右射一枪，在(5,7)的左面出现开口，(3,1)右面的开口消失。
5.走进(6,2)的开口，出来到(5,6)。用1单位时间。
6.向上射一枪，在(1,6)的下面出现开口。
7.经过秘密通道，走到(2,6)。用1单位时间。
8.走到终点。用1单位时间。

题解：//弱子的题解

考试时的思路是BFS，每走一步向四个方向打一枪，碰到传送门下一步就拓展之前所有的传送门的前面，这样打拿了60分。

正解是将这个地图转换成一个有向图，对于一个空地，如果上下左右为空地，就向它们连一条边权为1的边，然后我们从这个点出发，找到离它最近的一堵墙，向这堵墙前的一个空地连边，边权为这个空地到它的距离，然后向其他三个方向碰到第一堵墙前的空地连边，边权为刚才边权加一（因为使用传送门需要耗费1步），建图时，我们利用x，y坐标对一个点编号，方式随意，只要保证独立就行，建图完成后，在图上从起点向终点跑一边最短路即可。

但是需要注意，例如dis等数组开500*500是小了，我是开了600000才能过，前向星数组还需要乘上一个8，因为对于一个点我们最多有8次拓展。

代码：
 

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
#include<vector>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#define MAXA 600000
#define INF 0x3f3f3f3f
using namespace std;
int cnt,head[MAXA * 8];
struct Rx {
	int next,to,w;
}edge[MAXA * 8];
void Add(int u,int v,int w) {
	edge[++cnt].next = head[u];
	edge[cnt].to = v;
	edge[cnt].w = w;
	head[u] = cnt;
}
struct Node {
	int x,y,dis;
}Temp[5];
bool inq[MAXA];
int Tempidx,dis[MAXA],q[MAXA],h = 1,t = 1;
int n,m,Startx,Starty,Endx,Endy,Ans = 0x3f3f3f3f,Now;
char Map[505][505];
inline int Idx(int x,int y) {
	return (x - 1) * m + y;
}
bool cmp(Node a,Node b) {
	return a.dis < b.dis;
}
void Init(int x,int y) {
	Tempidx = Idx(x,y);
	if(Map[x-1][y] != '#') Add(Tempidx,Idx(x-1,y),1);
	if(Map[x+1][y] != '#') Add(Tempidx,Idx(x+1,y),1);
	if(Map[x][y-1] != '#') Add(Tempidx,Idx(x,y-1),1);
	if(Map[x][y+1] != '#') Add(Tempidx,Idx(x,y+1),1);
	int i = x,j = y;
	while(Map[i][j] != '#')
	    i--;
	Temp[1].x = i + 1;
	Temp[1].y = j;
	Temp[1].dis = x - i;
	i = x;
	while(Map[i][j] != '#')
	    i++;
	Temp[2].x = i - 1;
	Temp[2].y = j;
	Temp[2].dis = i - x;
	i = x;
	while(Map[i][j] != '#')
	    j--;
	Temp[3].x = i;
	Temp[3].y = j + 1;
	Temp[3].dis = y - j;
	j = y;
	while(Map[i][j] != '#')
	    j++;
	Temp[4].x = i;
	Temp[4].y = j - 1;
	Temp[4].dis = j - y;
	sort(Temp + 1,Temp + 5,cmp);
	
	Add(Tempidx,Idx(Temp[1].x,Temp[1].y),Temp[1].dis - 1);
	Add(Tempidx,Idx(Temp[2].x,Temp[2].y),Temp[1].dis);
	Add(Tempidx,Idx(Temp[3].x,Temp[3].y),Temp[1].dis);
	Add(Tempidx,Idx(Temp[4].x,Temp[4].y),Temp[1].dis);
}

void SPFA(int s) {
	int k;
	for(int i=1;i<=n*m;i++)
	    dis[i] = INF;
	dis[s] = 0;
	q[t++] = s;
	inq[s] = 1;
	while(h != t) {
		k = q[h];
		h++;
		inq[k] = 0;
		if(dis[k] == INF) continue;
		for(int i=head[k];i;i=edge[i].next) {
			int j = edge[i].to;
			if(dis[j] > dis[k] + edge[i].w) {
				dis[j] = dis[k] + edge[i].w;
				if(!inq[j]) {
					inq[j] = 1;
					q[t++] = j;
				}
			}
		}
	}
}
int main() {
//	freopen("portal.in","r",stdin);
//	freopen("portal.out","w",stdout);
	ios::sync_with_stdio(false);
	cin.tie(0);
	cin >> n >> m;
	for(int i=1;i<=n;i++)
	    for(int j=1;j<=m;j++) {
	    	cin >> Map[i][j];
	    	if(Map[i][j] == 'C') {
	    		Startx = i;
	    		Starty = j;
			}
			if(Map[i][j] == 'F') {
				Endx = i;
				Endy = j;
			}
		}
	for(int i=1;i<=n;i++)
	    for(int j=1;j<=m;j++)
	        if(Map[i][j] != '#')
	        	Init(i,j);
	SPFA(Idx(Startx,Starty));
	if(dis[Idx(Endx,Endy)] == INF)
	   printf("nemoguce");
	else printf("%d",dis[Idx(Endx,Endy)]);
}