【BZOJ 1611】 [Usaco2008 Feb]Meteor Shower流星雨

1611: [Usaco2008 Feb]Meteor Shower流星雨

Time Limit: 5 Sec Memory Limit: 64 MB

Submit: 1007 Solved: 444

[Submit][Status][Discuss]

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

HINT

Source

Silver

spfa。

预处理出每个位置,然后spfa求出最早到达每个点的时刻即可。

#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <queue>
#define mp make_pair
#define inf 0x3f3f3f3f
#define pa pair<int,int>
using namespace std;
int m,inq[305][305],a[305][305],d[305][305],fx[5][3];
queue<pa> q;
int ok(int x,int y)
{
if (x<0||y<0) return 0;
return 1;
}
int main()
{
scanf("%d",&m);
fx[1][1]=fx[2][1]=0,fx[1][2]=1,fx[2][2]=-1;
fx[3][2]=fx[4][2]=0,fx[3][1]=1,fx[4][1]=-1;
for (int i=0;i<=301;i++)
for (int j=0;j<=301;j++)
a[i][j]=inf,d[i][j]=inf;
for (int i=1;i<=m;i++)
{
int x,y,t;
scanf("%d%d%d",&x,&y,&t);
a[x][y]=min(a[x][y],t);
for (int j=1;j<=4;j++)
{
int nx=x+fx[j][1],ny=y+fx[j][2];
if (ok(nx,ny))
a[nx][ny]=min(a[nx][ny],t);
}
}

d[0][0]=0;
inq[0][0]=1;
q.push(mp(0,0));
int ans=inf;
if (a[0][0]==inf) ans=0;
while (!q.empty())
{
pa p=q.front();
q.pop();
int x=p.first,y=p.second;
inq[x][y]=0;
for (int i=1;i<=4;i++)
{
int nx=x+fx[i][1],ny=y+fx[i][2];
if (ok(nx,ny)&&d[nx][ny]>d[x][y]+1&&a[nx][ny]>d[x][y]+1)
{
d[nx][ny]=d[x][y]+1;
if (a[nx][ny]==inf) ans=min(ans,d[nx][ny]);
if (!inq[nx][ny])
q.push(mp(nx,ny)),inq[nx][ny]=1;
}
}
}
if (ans==inf) ans=-1;
cout<<ans<<endl;
return 0;
}