【2020.8.23NOIP模拟赛】传送爸爸【SPFA】

题目背景
$wdyhy$非常喜欢他的爸爸们，所以只要让爸爸们永远呆在$wdyhy$的迷宫里，$wdyhy$就可以永远爱着他们啦！

---

题目描述

$wdyhy$有一个 $R$ 行 $C$ 列的迷宫，每一个小格有一个字符。

KaTeX parse error: Expected 'EOF', got '#' at position 1: #̲(number sign) 表⽰一个墙块，

$.(dot)$ 表⽰一块空地，

$S(uppercaseletters)$ 表⽰你现在的位置，

$C(uppercaseletterc)$ 表⽰爸爸现在的位置。

你只能通过空地，并且，只有当两块空地有相临边时，你才可以从其中一个走向另一个。特别的，描述在地图里的矩形区域完全被墙块包围。

为了能够更快的到达爸爸的位置，你从迦勒底获得了一个传送枪，它的操作方式如下：任意时候它都可以向上下左右四个方向发射传送门。当一个传送门以某一特定方向发送时，它会沿着此方向飞出直到碰到第一个墙块。当发生上述情况时，一个传送门就会贴在面对着你那面的墙块上。

任意时刻最多只能存在两个传送门。如果迷宫里已经有了两个传送门，那么你可以使用传送枪的另一个功能马上移除一个（由你选定）。向一个已经存在的传送门发射时，后发射的传送门会替代原有的,所以每个墙块的每一面最多只会存在一个传送门。注意，可以存在两个传送门在同一墙块上，只要他们所朝方向不同。

当迷宫中存在两个传送门时，你可以通过它们传送自己。当你站在一个传送门旁时，你可以走进传送门然后从另一个传送门走出来到它相邻的空地，这样耗费的时间与你走相邻两个空地的时间是相同的。

你可以假设发射传送门不会耗费任何时间，走相邻两个空地和通过传送门的时间都是 $1$ 。 那么,你最少需要多长时间才能到爸爸那里。

---

输入格式

输入文件第一行包括两个整数：地图的行数 R 和列数 C。接下来的 R 行描述这个地图。每一行包括 C 个字符：KaTeX parse error: Expected 'EOF', got '#' at position 1: #̲,.,S或$C$（含义已在上文描述）。

输出格式

输出一个整数——从你现在位置出发到达爸爸的位置所需要的最少时间。

---

输入输出样例

输入 #1
4 4
.#.C
.#.#
…
S…

输出 #1
4

---

说明/提示

---

解题思路
这道题可以用 SPFA 解决，就是建图十分复杂。。。（还是我太菜了）

我们可以知道，一个地方可以直接走到另外一个地方，也可以走传送门。
直接走：处理起来比较简单，枚举四个方向就好了。。。
走传送门：我们枚举四个方向，分别走到第一个碰到墙壁的地方，记录墙壁前的那个点x，将现在所在位置与x相连。。。边权值为四个方向中最短的哪一个。。、

为了不TLE，我们可以用记忆化搜索，如果这个点曾来过，且方向一致，可以直接$return$之前求出来的值。。。

PS：很耗油也经常用的矩阵标号，第i行j列为：$（i-1)\ast$总列数$+j$

---

代码

#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const int dx[4]={1,0,-1,0};
const int dy[4]={0,1,0,-1};
int r,c,x,y,k,s,h,e,t,head[10000010],bsy[1010][1010],v[10000010],q[10000000],dis[10000000],rem[4][1000010],ch[10000010],end[4][1000010];
char aa;
struct cc{
	int x,next,w;
}a[10000000];
void add(int x,int y,int w){
	a[++k].x=y;
	a[k].next=head[x];
	a[k].w=w;
	head[x]=k;
}
bool check(int x,int y){
	if(x>0&&y>0&&x<=r&&y<=c)
		return 1;
	else return 0;
}
int dfs(int x,int y,int k){//现在在第i行j列，方向为k
	if(bsy[x][y]||!check(x,y))
	{
		ch[k]=(x-dx[k]-1)*c+(y-dy[k]);//记录当方向为k，一直走所到的墙壁前的位置
		end[k][(x-1)*c+y]=ch[k];//end表示从(x-1)*c+y出发，方向为k，一直走所到的墙壁前的位置
		return rem[k][(x-1)*c+y]=0;//end表示从(x-1)*c+y出发，方向为k，能走几步
	}
	if(rem[k][(x-1)*c+y])//记忆化
	{
		ch[k]=end[k][(x-1)*c+y];
		return rem[k][(x-1)*c+y];
	}
	int re=dfs(x+dx[k],y+dy[k],k)+1;
	end[k][(x-1)*c+y]=end[k][(x+dx[k]-1)*c+(y+dy[k])];
	return rem[k][(x-1)*c+y]=re;
}
void SPFA(){
	memset(dis,0X7f,sizeof(dis));
	v[s]=1;
	dis[s]=0;
	q[1]=s;
	h=0,t=1;
	while(h<t)
	{
		h++;
		int x=q[h];
		for(int i=head[x];i;i=a[i].next)
		{
			if(dis[x]+a[i].w<dis[a[i].x])
			{
				dis[a[i].x]=dis[x]+a[i].w;
				if(!v[a[i].x])
				{
					v[a[i].x]=1;
					q[++t]=a[i].x;
				}
			}
		}
		v[x]=0;
	}
}
int main(){
	scanf("%d%d",&r,&c);
	/*for(int i=1;i<=c;i++)
		bsy[0][i]=bsy[r+1][i]=1;
	for(int i=1;i<=r;i++)
		bsy[i][0]=bsy[i][c+1]=1;*/
	for(int i=1;i<=r;i++)
	{
		for(int j=1;j<=c;j++)
		{
			cin>>aa;
			if(aa=='#') bsy[i][j]=1;
			else if(aa=='S') s=(i-1)*c+j;
			else if(aa=='C') e=(i-1)*c+j;
		}
	}
	for(int i=1;i<=r;i++)
	{
		for(int j=1;j<=c;j++)
		{
			if(!bsy[i][j])//要不是墙
			{
				int far=2147483600,str=(i-1)*c+j;
				for(int k=0;k<4;k++)
				{
					far=min(far,dfs(i+dx[k],j+dy[k],k)+1);//求出四个方向离墙的距离的最小值
					if(!bsy[i+dx[k]][j+dy[k]]&&check(i+dx[k],j+dy[k]))
						add(str,(i+dx[k]-1)*c+(j+dy[k]),1);//正常走，建边
				}
				for(int k=0;k<4;k++)//传送,建边
				{
					if(str!=ch[k])
						add(str,ch[k],far);
				}
			}
		}
	}
	SPFA();
	printf("%d",dis[e]);
}