BZOJ 1611: [Usaco2008 Feb]Meteor Shower流星雨

 

1611: [Usaco2008 Feb]Meteor Shower流星雨

Description

去年偶们湖南遭受N年不遇到冰冻灾害，现在芙蓉哥哥则听说另一个骇人听闻的消息： 一场流星雨即将袭击整个霸中，由于流星体积过大，它们无法在撞击到地面前燃烧殆尽， 届时将会对它撞到的一切东西造成毁灭性的打击。很自然地，芙蓉哥哥开始担心自己的 安全问题。以霸中至In型男名誉起誓，他一定要在被流星砸到前，到达一个安全的地方 （也就是说，一块不会被任何流星砸到的土地）。如果将霸中放入一个直角坐标系中， 芙蓉哥哥现在的位置是原点，并且，芙蓉哥哥不能踏上一块被流星砸过的土地。根据预 报，一共有M颗流星(1 <= M <= 50,000)会坠落在霸中上，其中第i颗流星会在时刻 T_i (0 <= T_i <= 1,000)砸在坐标为(X_i, Y_i) (0 <= X_i <= 300；0 <= Y_i <= 300) 的格子里。流星的力量会将它所在的格子，以及周围4个相邻的格子都化为焦土，当然 芙蓉哥哥也无法再在这些格子上行走。芙蓉哥哥在时刻0开始行动，它只能在第一象限中， 平行于坐标轴行动，每1个时刻中，她能移动到相邻的（一般是4个）格子中的任意一个， 当然目标格子要没有被烧焦才行。如果一个格子在时刻t被流星撞击或烧焦，那么芙蓉哥哥 只能在t之前的时刻在这个格子里出现。请你计算一下，芙蓉哥哥最少需要多少时间才能到 达一个安全的格子。

Input

* 第1行: 1个正整数：M * 第2..M+1行: 第i+1行为3个用空格隔开的整数：X_i，Y_i，以及T_i

Output

输出1个整数，即芙蓉哥哥逃生所花的最少时间。如果芙蓉哥哥无论如何都无法在流星雨中存活下来，输出-1

Sample Input

4
0 0 2
2 1 2
1 1 2
0 3 5
输入说明:
一共有4颗流星将坠落在霸中，它们落地点的坐标分别是(0, 0)，(2, 1)，(1, 1)
以及(0, 3)，时刻分别为2，2，2，5。

* = 流星落点 # = 行走禁区
t = 0 t = 2 t = 5
5|. . . . . . . 5|. . . . . . . 5|. . . . . . .
4|. . . . . . . 4|. . . . . . . 4|# . . . . . .
3|. . . . . . . 3|. . . . . . . 3|* # . . . . .
2|. . . . . . . 2|. # # . . . . 2|# # # . . . .
1|. . . . . . . 1|# * * # . . . 1|# # # # . . .
0|B . . . . . . 0|* # # . . . . 0|# # # . . . .
-------------- -------------- --------------
0 1 2 3 4 5 6 0 1 2 3 4 5 6 0 1 2 3 4 5 6

Sample Output

5
输出说明:
如果我们观察在t=5时的霸中，可以发现离芙蓉哥哥最近的安全的格子是
(3,0)——不过由于早在第二颗流星落地时，芙蓉哥哥直接跑去(3,0)的路线就被封死了。
离芙蓉哥哥第二近的安全格子为(4,0)，但它的情况也跟(3,0)一样。再接下来的格子就是在
(0,5)-(5,0)这条直线上。在这些格子中，(0,5)，(1,4)以及(2,3)都能在5个单位时间内到达。

5|. . . . . . .
4|. . . . . . .
3|3 4 5 . . . . 某个合法的逃生方案中
2|2 . . . . . . 芙蓉哥哥每个时刻所在地点
1|1 . . . . . .
0|0 . . . . . .
--------------
0 1 2 3 4 5 6

 

 

——华丽的分割线——

这道题目就是一开始求出所有安全的地方，然后模拟就可以了。人先向外扩张一个，然后在用流星毁一下地图，重复到人必死或者安全为止、

有个心态很不好的地方就是、这道题目如果走出300*300这个区域，是算安全的！【卧槽！害我Wa了好多次，这不是坑爹吗！什么心态！

代码：

#include<cstdio>

#include<algorithm>

#include<queue>

using namespace std;



int Dx[]={0,0,1,0,-1};

int Dy[]={0,1,0,-1,0}; 



bool isDestroyed[310][310];

bool mapDestroyed[310][310];

bool inQueue[310][310];



int m;

struct TPointNode{

	int x;

	int y;

	int t;

};

TPointNode point[50010];

int index;

inline bool cmp(TPointNode a,TPointNode b){

	if (a.t<b.t) return true;

	return false;

}



struct PointNode{

	PointNode(){}

	PointNode(int a,int b){

		x=a;y=b;

	}

	int x;

	int y;

};

queue<PointNode> que;

int NowTime;



inline void getDestroyPlace(){

	for (int i=1;i<=m;i++){

		for (int k=0;k<=4;k++){

			int x=Dx[k]+point[i].x,y=Dy[k]+point[i].y;

			if (0<=x && x<=300 && 0<=y && y<=300) isDestroyed[x][y]=true;

		}

	}	

}



inline bool BFS(){

	queue<PointNode> inque;

	while(!que.empty()){

		PointNode NowNode=que.front();

		int nx=NowNode.x,ny=NowNode.y;

		que.pop();

		if (mapDestroyed[nx][ny]==true) continue;

		for (int k=1;k<=4;k++){

			int x=Dx[k]+nx,y=Dy[k]+ny;

			if (0<=x && 0<=y){

				if (isDestroyed[x][y]==false) {

					return true;

				}

				if (x>300 || y>300) return true;

				if (inQueue[x][y]==false){

					inQueue[x][y]=true;

					inque.push(PointNode(x,y));

				}

			}

		}

	}

	que=inque;

	return false;

}



int main(){

	scanf("%d\n",&m);

	for (int i=1;i<=m;i++) scanf("%d%d%d",&point[i].x,&point[i].y,&point[i].t);

	

	getDestroyPlace();

	

	sort(point+1,point+m+1,cmp);

	NowTime=0;

	index=1;

	que.push(PointNode(0,0));

	while(point[index].t==0){

		for (int k=0;k<=4;k++){

			int x=Dx[k]+point[index].x,y=Dy[k]+point[index].y;

			if (0<=x && x<=300 && 0<=y && y<=300) mapDestroyed[x][y]=true;

		}

		index++;

	}

	for (;index<=m;){

		NowTime++;

		if (BFS()) {

			printf("%d\n",NowTime);

			return 0;

		}

		while(point[index].t==NowTime){

			for (int k=0;k<=4;k++){

				int x=Dx[k]+point[index].x,y=Dy[k]+point[index].y;

				if (0<=x && x<=300 && 0<=y && y<=300) mapDestroyed[x][y]=true;

			}

			index++;

		}

	}

	printf("-1\n");

	

	return 0;

}