poj Meteor Shower - 搜索
Time Limit: 1000MS | Memory Limit: 65536K | |
Total Submissions: 16313 | Accepted: 4291 |
Description
Bessie hears that an extraordinary meteor shower is coming; reports say that these meteors will crash into earth and destroy anything they hit. Anxious for her safety, she vows to find her way to a safe location (one that is never destroyed by a meteor) . She is currently grazing at the origin in the coordinate plane and wants to move to a new, safer location while avoiding being destroyed by meteors along her way.
The reports say that M meteors (1 ≤ M ≤ 50,000) will strike, with meteor i will striking point (Xi, Yi) (0 ≤ Xi ≤ 300; 0 ≤ Yi ≤ 300) at time Ti (0 ≤ Ti ≤ 1,000). Each meteor destroys the point that it strikes and also the four rectilinearly adjacent lattice points.
Bessie leaves the origin at time 0 and can travel in the first quadrant and parallel to the axes at the rate of one distance unit per second to any of the (often 4) adjacent rectilinear points that are not yet destroyed by a meteor. She cannot be located on a point at any time greater than or equal to the time it is destroyed).
Determine the minimum time it takes Bessie to get to a safe place.
Input
* Line 1: A single integer: M
* Lines 2..M+1: Line i+1 contains three space-separated integers: Xi, Yi, and Ti
Output
* Line 1: The minimum time it takes Bessie to get to a safe place or -1 if it is impossible.
Sample Input
4 0 0 2 2 1 2 1 1 2 0 3 5
Sample Output
5
Source
1 /* 2 * poj.org 3 * Problem#3669 4 * Accepted 5 * Time:266ms 6 * Memory:620k 7 */ 8 #include<iostream> 9 #include<cstdio> 10 #include<fstream> 11 #include<cstring> 12 #include<queue> 13 #include<algorithm> 14 #define INF 0xfffffff 15 #define smin(a, b) a = min(a, b); 16 using namespace std; 17 typedef bool boolean; 18 template<typename T> 19 inline void readInteger(T& u){ 20 char x; 21 while(!isdigit((x = getchar()))); 22 for(u = x - '0'; isdigit((x = getchar())); u = (u << 1) + (u << 3) + x - '0'); 23 ungetc(x, stdin); 24 } 25 template<typename T>class Matrix{ 26 public: 27 T *p; 28 int lines; 29 int rows; 30 Matrix():p(NULL){ } 31 Matrix(int lines, int rows):lines(lines), rows(rows){ 32 p = new T[(lines * rows)]; 33 } 34 T* operator [](int pos){ 35 return (p + pos * lines); 36 } 37 }; 38 int n; 39 Matrix<boolean> visited; 40 Matrix<int> down; 41 inline void init(){ 42 visited = Matrix<boolean>(303, 303); 43 down = Matrix<int>(303, 303); 44 fill(down.p, down.p + 303 * 303, INF); 45 memset(visited.p, false, sizeof(boolean) * 303 * 303); 46 readInteger(n); 47 for(int i = 1, a, b, c; i <= n; i++){ 48 readInteger(a); 49 readInteger(b); 50 readInteger(c); 51 smin(down[a][b], c); 52 if(a - 1 >= 0) smin(down[a - 1][b], c); 53 smin(down[a + 1][b], c); 54 if(b - 1 >= 0) smin(down[a][b - 1], c); 55 smin(down[a][b + 1], c); 56 } 57 } 58 typedef class Point{ 59 public: 60 int x; 61 int y; 62 int step; 63 Point(const int x = 0, const int y = 0, const int step = 0):x(x), y(y), step(step){} 64 }Point; 65 int result = -1; 66 queue<Point> que; 67 const int move[2][4] = {{1, -1, 0, 0}, {0, 0, 1, -1}}; 68 void solve(){ 69 if(down[0][0] == INF){ 70 result = 0; 71 return; 72 } 73 que.push(Point()); 74 visited[0][0] = true; 75 while(!que.empty()){ 76 Point e = que.front(); 77 que.pop(); 78 for(int i = 0; i < 4; i++){ 79 Point eu = Point(e.x + move[0][i], e.y + move[1][i], e.step + 1); 80 if(eu.x >= 0 && eu.y >= 0 && eu.step < down[eu.x][eu.y] && !visited[eu.x][eu.y]){ 81 if(down[eu.x][eu.y] == INF){ 82 result = eu.step; 83 return; 84 } 85 visited[eu.x][eu.y] = true; 86 que.push(eu); 87 } 88 } 89 } 90 } 91 int main(){ 92 init(); 93 solve(); 94 printf("%d", result); 95 return 0; 96 }