【简●解】[HNOI2005]星际贸易

前提是你 提交于 2019-11-29 00:08:00

【大意】

太多了,懒得打,贴\(LG\)的图了。。。

【分析】

开始拿到这道题有点慌:怎么限制条件这么多,再读读题。

注意一个东西,就是贸易额与费用是独立分开的,并且题目保证只有一种方案获得最大贸易额。

所以我们\(dp\)也可以分开\(dp\)

对于贸易额,每走到一个星球,我可以选择卖或不卖,经典的背包问题,直接\(f[i][j]=min\{f[i-1][j-A_i]+B_i\}\)

然后从后往前扫记录决策点的转移,这些星球是必到的,即必定会在该星球上维修。

现在来看费用。

\(g[i][j]\)表示到第\(i\)个星球,剩余反物质燃料为\(j\)的最小费用,暴力枚举转移点,则有\(g[i][j]=\{g[k][l]+(j-l+2)\times P_i+T_i\}(0\le k<i,L_i-L_k\le L_0,l\le j)\)

\(O(n^4)\)转移。。。

然后可推出一个更优的式子:\(g[i][j]=min\{g[k][j+2]+F_i,g[i][j-1]+P_i\}\)

\(O(n^3)\)转移。。。

然后,,,我就嫖题解了。。。。

噢,原来可以单调队列优化啊,对于每个\(j\)开一个单调队列,若当前的\(i\)为必到的星球,就清空队列中所有的元素,然后把\(i\)放进去。

\(O(n^2)\)转移!!!

【Code】

#include <queue> #include <vector> #include <cstdio> #include <cstring> #include <algorithm> #define Re register using namespace std; const int INF = 0x3f3f3f3f; const int MAXR = 4000 + 5; const int N = 2000 + 5; const int M = 2000 + 5; inline int read(){     int f = 1, x = 0; char ch;     do { ch = getchar(); if (ch == '-') f = -1; } while (ch < '0' || ch > '9');     do {x = (x << 3) + (x << 1) + ch - '0'; ch = getchar(); } while (ch >= '0' && ch <= '9');      return f * x; } int n, m, r, L0, maxm, maxr, A[N], B[N], L[N], P[N], F[N], f[N][M], g[N][M], vis[N];  deque <int> Q[MAXR]; int main(){     n = read(), m = read(), r = read(), L0 = read();      if (r > (n << 1)) r = (n << 1);     for (int i = 1;i <= n; ++i) {         A[i] = read(), B[i] = read(), L[i] = read(), P[i] = read(), F[i] = read();         if (L[i] - L[i - 1] > L0) {             puts("Poor Coke!");             return 0;         }     }     memset(f, 128, sizeof f);     f[0][0] = 0;     for (int i = 1;i <= n; ++i) {         for (int j = 0;j <= m; ++j) {             if (f[i - 1][j] >= 0) f[i][j] = f[i - 1][j];             if (j >= A[i]) {                 if (f[i][j] < f[i - 1][j - A[i]] + B[i]) {                     f[i][j] = f[i - 1][j - A[i]] + B[i];                 }             }         }     }          for (int i = 1;i <= m; ++i) if (f[n][i] > f[n][maxm]) maxm = i;           for (int i = n, val = maxm;i >= 1; --i) {         if (val - A[i] >= 0 && f[i][val] == f[i - 1][val - A[i]] + B[i]){             vis[i] = 1;             val -= A[i];         }     }          memset(g, 0x3f, sizeof g);     g[0][r] = 0;     Q[r].push_back(0);     for (int i = 1;i <= n; ++i) {         for (int j = 0;j <= r; ++j) {             if (P[i] > 0 && j > 0) {                 g[i][j] = min(g[i][j], g[i][j - 1] + P[i]);             }                          while (!Q[j].empty() && L[i + 1] - L[Q[j].front()] > L0) Q[j].pop_front();                          if (!Q[j + 2].empty()) {                 g[i][j] = min(g[i][j], g[Q[j + 2].front()][j + 2] + F[i]);             }                          if (vis[i]) Q[j].clear();                          while (!Q[j].empty() && g[Q[j].back()][j] >= g[i][j]) Q[j].pop_back();                          Q[j].push_back(i);         }     }     for (int i = 1;i <= r; ++i) if (g[n][i] < g[n][maxr]) maxr = i;          if (g[n][maxr] == INF) printf("Poor Coke!\n");     else printf("%d %d\n", f[n][maxm], f[n][maxm] - g[n][maxr]);              return 0; }
标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!