当时区域赛的一道题。题意大概是这样的,有一个1~N的图,然后你要从1-》N,其中每经过一条边需要消耗你的时间和金钱,每到一个地方可以选择什么都不做,或者买一包盐,卖一包盐,身上不能同时有超过B包盐,然后你还可以从当前所在的世界i转道(i+1)modK的世界,每次转换花费时间1,然后问你在T时间内到达n点的时候你最多能有多少钱。 那个时候不会做是因为不知道如何控制dp的方向,我记得那个时候我已经把状态写好了,转移也是知道的,就是不知道转移的方向,后来发现时间的那一维是降的,于是我们就可以建一个优先队列,把一开始的状态压进去,然后拓展新的状态,这样就可以做出来了,调了很久,因为我心里认为它是无向图。。。。。吸取教训呀要!
#pragma warning(disable:4996)
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<algorithm>
#include<queue>
#include<cmath>
#define maxn 150
using namespace std;
int n, m, B, K, R, T;
int p[6][maxn];
struct Edge
{
int u, v, t, m;
Edge(){}
Edge(int ui, int vi, int ti, int mi) :u(ui), v(vi), t(ti), m(mi){}
};
vector<Edge> G[maxn];
int dp[6][maxn][5][250]; // dp[K][N][B][T] k维空间,第n个点,有b袋盐,剩余时间为T时的最多有多少钱
bool vis[6][maxn][5][250];
struct Node
{
int k, n, b, t;
Node(){}
Node(int ki, int ni, int bi, int ti) :k(ki), n(ni), b(bi), t(ti){}
bool operator < (const Node X) const{
return t < X.t;
}
};
int main()
{
int cas; cin >> cas; int ca = 0;
while (cas--)
{
scanf("%d%d%d%d%d%d", &n, &m, &B, &K, &R, &T);
for (int i = 0; i < K; i++){
for (int j = 1; j <= n; j++){
scanf("%d", &p[i][j]);
}
}
for (int i = 0; i <= n; i++) G[i].clear();
int ui, vi, ti, mi;
for (int i = 0; i < m; i++){
scanf("%d%d%d%d", &ui, &vi, &ti, &mi);
G[ui].push_back(Edge(ui, vi, ti, mi));
}
memset(dp, -1, sizeof(dp));
memset(vis, 0, sizeof(vis));
dp[0][1][0][T] = R;
vis[0][1][0][T] = true;
priority_queue<Node> que;
que.push(Node(0, 1, 0, T));
while (!que.empty()){
Node nxt = que.top(); que.pop();
if (nxt.k == 0 && nxt.n == n) continue;
int val = dp[nxt.k][nxt.n][nxt.b][nxt.t];
for (int i = 0; i < G[nxt.n].size(); i++){
Edge vv = G[nxt.n][i];
if (nxt.k != 0 && (vv.v == 1 || vv.v == n)) continue;
if (nxt.t >= vv.t&&val >= vv.m){
dp[nxt.k][vv.v][nxt.b][nxt.t - vv.t] = max(dp[nxt.k][vv.v][nxt.b][nxt.t - vv.t], val - vv.m);
if (!vis[nxt.k][vv.v][nxt.b][nxt.t - vv.t]){
que.push(Node(nxt.k, vv.v, nxt.b, nxt.t - vv.t));
vis[nxt.k][vv.v][nxt.b][nxt.t - vv.t] = true;
}
if (vv.v == 1 || vv.v == n) continue;
if (nxt.b < B && (val - vv.m) >= p[nxt.k][vv.v]){
dp[nxt.k][vv.v][nxt.b + 1][nxt.t - vv.t] = max(dp[nxt.k][vv.v][nxt.b + 1][nxt.t - vv.t],
val - vv.m - p[nxt.k][vv.v]);
if (!vis[nxt.k][vv.v][nxt.b + 1][nxt.t - vv.t]){
que.push(Node(nxt.k, vv.v, nxt.b + 1, nxt.t - vv.t));
vis[nxt.k][vv.v][nxt.b + 1][nxt.t - vv.t] = true;
}
}
if (nxt.b>0){
dp[nxt.k][vv.v][nxt.b - 1][nxt.t - vv.t] = max(dp[nxt.k][vv.v][nxt.b - 1][nxt.t - vv.t],
val - vv.m + p[nxt.k][vv.v]);
if (!vis[nxt.k][vv.v][nxt.b - 1][nxt.t - vv.t]){
que.push(Node(nxt.k, vv.v, nxt.b - 1, nxt.t - vv.t));
vis[nxt.k][vv.v][nxt.b - 1][nxt.t - vv.t] = true;
}
}
}
}
if (nxt.t < 1) continue;
if (nxt.n == 1 || nxt.n == n) continue;
dp[(nxt.k + 1) % K][nxt.n][nxt.b][nxt.t - 1] = max(dp[(nxt.k + 1) % K][nxt.n][nxt.b][nxt.t - 1],
val);
if (!vis[(nxt.k + 1) % K][nxt.n][nxt.b][nxt.t - 1]){
que.push(Node((nxt.k + 1) % K, nxt.n, nxt.b, nxt.t - 1));
vis[(nxt.k + 1) % K][nxt.n][nxt.b][nxt.t - 1] = true;
}
if (nxt.b < B&&val >= p[(nxt.k + 1) % K][nxt.n]){
dp[(nxt.k + 1) % K][nxt.n][nxt.b + 1][nxt.t - 1] = max(dp[(nxt.k + 1) % K][nxt.n][nxt.b + 1][nxt.t - 1],
val - p[(nxt.k + 1) % K][nxt.n]);
if (!vis[(nxt.k + 1) % K][nxt.n][nxt.b + 1][nxt.t - 1]){
que.push(Node((nxt.k + 1) % K, nxt.n, nxt.b + 1, nxt.t - 1));
vis[(nxt.k + 1) % K][nxt.n][nxt.b + 1][nxt.t - 1] = true;
}
}
if (nxt.b>0){
dp[(nxt.k + 1) % K][nxt.n][nxt.b - 1][nxt.t - 1] = max(dp[(nxt.k + 1) % K][nxt.n][nxt.b - 1][nxt.t - 1],
val + p[(nxt.k + 1) % K][nxt.n]);
if (!vis[(nxt.k + 1) % K][nxt.n][nxt.b - 1][nxt.t - 1]){
que.push(Node((nxt.k + 1) % K, nxt.n, nxt.b - 1, nxt.t - 1));
vis[(nxt.k + 1) % K][nxt.n][nxt.b - 1][nxt.t - 1] = true;
}
}
}
int ans = -1;
for (int i = 0; i <= B; i++){
for (int j = 0; j <= T; j++){
ans = max(ans, dp[0][n][i][j]);
}
}
printf("Case #%d: ", ++ca);
if (ans == -1) puts("Forever Alone");
else printf("%d\n", ans);
}
return 0;
}