UVA1628 送匹萨 （有关未来费用的区间dp 升级版）

题意：话说正义的使者boshi在修好长城后开了一家匹萨店，邪恶的大膜王jyf为了报复boshi，用他奇妙的膜法分身n个来”骚扰“boshi。 boshi和jyf所在的小镇有一条笔直的街道，0位置是boshi的匹萨店，jyf的分身分别站在x1,x2...xnx1,x2...xn位置上（题目是按照从左到右给出坐标的，boshi走一个单位坐标花费一个单位时间），jyf的订单对boshi提出了这样的要求：如果boshi在时间t送到xixi处，可获得ei−tei−t元钱。由于boshi的匹萨店没有其他店员，所以boshi开着小电驴左右走一定很花时间，到时候jyf就可以倒收boshi的钱辣。 然而机制boshi早就识破了jyf的阴谋，如果发生jyf可能可以倒收他的钱的情况，他跳过这个订单不就行了？那么boshi最多可以坑jyf多少钱呢？

分析：紫书P303有详细讲解，本题与修缮长城一题有所相似。所以解法有相似之处。不同之处就是本体可能会产生负情况，即送餐时间晚了客户会反过来找你要钱所以需要放弃，但修缮长城只有费用，顺手修了肯定是一个不错的选择。依旧将区间两端与位置作为状态不过要添加一维cnt表示还需要送餐的人数。类似地定义：d[i][j][cnt][p]表示已经送完了ij区间（区间内或送餐或放弃）位于p（p==0||p==1）还剩下cnt个客户需要继续送餐。添加的一维成功解决了不知道还有多少的客户需要送餐的问题。这里注意到DP过程中利用的信息都来源于状态，因此定义的状态必须要提供转移足够的信息，这样才能得到所需要的值。

转移方程：d[i][j][cnt][p]=max{d[i][k][cnt-1][1]+pay[knew]//k在余下的右区间  ,  d[k][j][cnt-1][0] +pay[knew]//k在余下的左区间}

参考：https://blog.csdn.net/hahalidaxin/article/details/51017196。
代码：

#include<bits/stdc++.h>
using namespace std;

const int maxn = 100 + 5;

int kase, n;
int p[maxn], v[maxn];
int d[maxn][maxn][maxn][2];
int vis[maxn][maxn][maxn][2];

// already considered s~e, still need to delivery to cnt people.
// pos = 0 means at s, pos = 1 means at e
int dp(int s, int e, int cnt, int pos) {
  if(cnt == 0) return 0;

  int &ans = d[s][e][cnt][pos];
  if(vis[s][e][cnt][pos] == kase) return ans;
  vis[s][e][cnt][pos] = kase;

  ans = 0;
  if(!pos) {
    for(int i = 0; i < s; i++)
      ans = max(ans, v[i] - cnt * abs(p[i] - p[s]) + dp(i, e, cnt - 1, 0));
    for(int i = e + 1; i < n; i++)
      ans = max(ans, v[i] - cnt * abs(p[i] - p[s]) + dp(s, i, cnt - 1, 1));
  }
  else {
    for(int i = 0; i < s; i++)
      ans = max(ans, v[i] - cnt * abs(p[i] - p[e]) + dp(i, e, cnt - 1, 0));
    for(int i = e + 1; i < n; i++)
      ans = max(ans, v[i] - cnt * abs(p[i] - p[e]) + dp(s, i, cnt - 1, 1));
  }
  return ans;
}

int main() {
  int T;
  scanf("%d",&T);
  memset(vis, 0, sizeof(vis));
  for(kase = 1; kase <= T; kase++) {
    scanf("%d", &n);
    for(int i = 0; i < n; i++) scanf("%d", &p[i]);
    for(int i = 0; i < n; i++) scanf("%d", &v[i]);

    int ans = 0;
    for(int k = 1; k <= n; k++)
      for(int i = 0; i < n; i++)
        ans = max(ans, v[i] - k * abs(p[i]) + dp(i, i, k - 1, 0));
    printf("%d\n",ans);
  }
  return 0;
}