紫书上动态规划例题,也是2003的Finals题。
紫书上说的十分浅显易懂。
对于这道题,我们首先要明白每一秒我们可以干什么。我们可以去前一个车站,可以去下一个车站,也可以在原地等一秒。
而影响我们选择哪种方案的要素为当前时间和所处的车站。
所以这两个就成为我们要不断调节的元素。我们采用dp[i][j]代表当前是第i分钟,所处第j个车站,我们希望停留的时间尽量少,所以我们dp里面每个位置记录的就是当前情况下所花费的最小等待时间。
我们首先令dp[T][n]为0,因为你必须在第T时刻恰好到达n车站。然后就是不断向前递推了。
每一秒在某个车站都满足,dp[i][j] = dp[i + 1][j] + 1.这样就相当于等了一秒。
但是如果当前有地铁经过的话,我们可以选择上与不上,此时递推式满足dp[i][j] = min(dp[i][j], dp[i + t[j]][j + 1])去下一个车站或者是反向的dp[i][j] = min(dp[i][j], dp[i - t[j - 1]][j - 1])去前一个车站。
至于如何判断当前时刻当前车站是否有地铁是很简单了。
如果最后发现第一个车站是没有更新过的,那就说明没有办法在指定时间内从车站1到车站n,直接输出impossible、
代码如下:
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<utility>
#include<stack>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
#include<set>
#include<map>
using namespace std;
const int maxn = 300;
const int INF = 0x3f3f3f3f;
int n, T, m1, m2;
int d1[maxn], d2[maxn], t[maxn];
int dp[maxn][maxn];
bool has_train[maxn][maxn][2];
int main() {
#ifndef ONLINE_JUDGE
freopen("in.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif
int Case = 1;
while(scanf("%d", &n) != EOF) {
if(n == 0)
break;
memset(has_train, 0, sizeof(has_train));
memset(dp, 0, sizeof(dp));
scanf("%d", &T);
for(int i = 1; i < n; i++) {
scanf("%d", &t[i]);
}
scanf("%d", &m1);
for(int i = 0; i < m1; i++) {
scanf("%d", &d1[i]);
}
scanf("%d", &m2);
for(int i = 0; i < m2; i++) {
scanf("%d", &d2[i]);
}
for(int i = 0; i < m1; i++) {
int tmp = d1[i];
for(int j = 1; j <= n; j++) {
has_train[tmp][j][0] = 1;
tmp += t[j];
// if(tmp > T)
// break;
}
}
for(int i = 0; i < m2; i++) {
int tmp = d2[i];
for(int j = n; j >= 1; j--) {
has_train[tmp][j][1] = 1;
tmp += t[j - 1];
// if(tmp > T)
// break;
}
}
for(int i = 1; i <= n - 1; i++)
dp[T][i] = INF;
dp[T][n] = 0;
for(int i = T - 1; i >= 0; i--) {
for(int j = 1; j <= n; j++) {
dp[i][j] = dp[i + 1][j] + 1;
if(j < n && has_train[i][j][0] && i + t[j] <= T)
dp[i][j] = min(dp[i][j], dp[i + t[j]][j + 1]);
if(j > 1 && has_train[i][j][1] && i + t[j - 1] <= T)
dp[i][j] = min(dp[i][j], dp[i + t[j - 1]][j - 1]);
}
}
printf("Case Number %d: ", Case++);
if(dp[0][1] >= INF)
printf("impossible\n");
else
printf("%d\n", dp[0][1]);
}
return 0;
}