ZOJ 2278 Fight for Food(LIS)

题意:给出主人公的初始位置,老鼠的出现位置和时间,主人公每一个单位时间只能向邻近的格子走一步,如果走到某个格子上的时候老鼠正好出现在那个位置,那么就可以抓住它,求主人公最多能抓到多少老鼠.

看起来三维的结构,其实只需要时间这一维就够了,因为两个格子的最短距离必然是曼哈顿距离,距离的计算可以是常数的,那么就相当于转换成了类似LIS的问题了.

预处理的时候可以把每个点到其他所有点的所需时间都找出来存到dis[i1][j1][i2][j2]里.

然后按时间对老鼠排个序,因为主人公的出现时间相当于0,所以我这里做的是不升序,再进行普通的O(n^2)的LIS算法dp[i] = max(dp[j]) + 1, 0<=j<i.

因为这里的P可能有30000, 普通的LIS会超时,所以必须要优化.

这里提几个优化:

1.如果主人公到不了某个格子(被石头挡住了),那么出现在这个位置上的老鼠就可以删除掉了.

2.预处理的时候计算出每个格子能到的最远格子的所需时间max_dis[i][j],再在LIS算法中维护一个maxv[i]代表dp[0]到dp[i]的最大值,在LIS算法的内层循环里:

由于老鼠出现时间是降序的,当两个格子i, j的老鼠出现时间差大于max_dis[i][j]的时候,就意味着这个点可以走到任何的k<=j的格子位置上,这里就可以break出来了.

只需要求一下dp[i] = max(dp[i], dp[j] + 1)就肯定是最优了.

#include <cstdio>#include <queue>#include <algorithm>#include <memory.h>#include <cmath>using namespace std;const int INF = 0x6f6f6f6f;const int MAX = 30052;struct Node{int i, j, step;Node(int i_ = 0, int j_ = 0, int step_ = 0):i(i_), j(j_), step(step_){}};struct State{int i, j, time, dp;State(int i_ = 0, int j_ = 0, int time_ = 0, int dp_ = 0):i(i_), j(j_), time(time_), dp(dp_){}bool operator<(const State & rhs)const{return time > rhs.time;}}s[MAX];char mat[11][11];bool vis[11][11];int dis[11][11][11][11];int max_dis[11][11], maxv[MAX];int dir[4][2] = {{0, 1}, {0, -1}, {1, 0}, {-1, 0}};int N, M, P, si, sj;void bfs(int si, int sj, int dis[11][11]){memset(vis, false, sizeof(vis));queue<Node> q;q.push(Node(si, sj, 0));vis[si][sj] = true;while(q.size()){Node f = q.front();q.pop();dis[f.i][f.j] = f.step;max_dis[si][sj] = f.step;for(int i = 0; i < 4; ++i){int ni = f.i + dir[i][0], nj = f.j + dir[i][1];if(ni >= 0 && ni < N && nj >= 0 && nj < M && mat[ni][nj] != '#' && !vis[ni][nj]){vis[ni][nj] = true;q.push(Node(ni, nj, f.step + 1));}}}}void preprocess(){memset(dis, 0x6f, sizeof(dis));for(int i = 0; i < N; ++i){for(int j = 0; j <= M; ++j){if(mat[i][j] != '#'){if(mat[i][j] == 'L'){si = i, sj = j;}bfs(i, j, dis[i][j]);}}}}int main(int argc, char const *argv[]){while(scanf("%d%d", &N, &M) == 2){memset(max_dis, 0, sizeof(max_dis));memset(maxv, 0, sizeof(maxv));for(int i = 0; i < N; ++i){scanf("%s", mat[i]);}preprocess();int TP;P = 0;scanf("%d", &TP);for(int i = 0; i < TP; ++i){int x, y, T;scanf("%d%d%d", &x, &y, &T);s[P] = State(x - 1, y - 1, T);if(dis[si][sj][x - 1][y - 1] != INF){//can't reach that grid++P;}}s[P++] = State(si, sj, 0);sort(s, s + P);for(int i = 0; i < P; ++i){s[i].dp = 1;int j = i - 1;for(; j >= 0 && s[j].time - s[i].time < max_dis[s[i].i][s[i].j]; --j){if(dis[s[i].i][s[i].j][s[j].i][s[j].j] <= s[j].time - s[i].time && s[j].dp + 1 > s[i].dp){s[i].dp = s[j].dp + 1;}}if(j >= 0 && s[i].dp < maxv[j] + 1){s[i].dp = maxv[j] + 1;}if(i == 0)maxv[i] = s[i].dp;else maxv[i] = max(maxv[i - 1], s[i].dp);}printf("%d\n", s[P - 1].dp - 1);}return 0;}