[SDOI2008] Sue 的小球（三维版区间dp）：第一弹

[SDOI2008] Sue 的小球

题目描述

Sue 和 Sandy 最近迷上了一个电脑游戏，这个游戏的故事发在美丽神秘并且充满刺激的大海上，Sue 有一支轻便小巧的小船。然而，Sue 的目标并不是当一个海盗，而是要收集空中漂浮的彩蛋，Sue 有一个秘密武器，只要她将小船划到一个彩蛋的正下方，然后使用秘密武器便可以在瞬间收集到这个彩蛋。然而，彩蛋有一个魅力值，这个魅力值会随着彩蛋在空中降落的时间而降低，Sue 要想得到更多的分数，必须尽量在魅力值高的时候收集这个彩蛋，而如果一个彩蛋掉入海中，它的魅力值将会变成一个负数，但这并不影响 Sue 的兴趣，因为每一个彩蛋都是不同的，Sue 希望收集到所有的彩蛋。

然而 Sandy 就没有 Sue 那么浪漫了，Sandy 希望得到尽可能多的分数，为了解决这个问题，他先将这个游戏抽象成了如下模型：

将大海近似的看做 $x$ 轴，以 Sue 所在的初始位置作为坐标原点建立一个竖直的平面直角坐标系。

一开始空中有 $N$ 个彩蛋，对于第 $i$ 个彩蛋，他的初始位置用整数坐标 $(x_i,y_i)$ 表示，游戏开始后，它匀速沿 $y$ 轴负方向下落,速度为 $v_i$ 单位距离/单位时间。Sue 的初始位置为 $(x_0,0)$，Sue 可以沿 $x$ 轴的正方向或负方向移动，Sue 的移动速度是 $1$ 单位距离/单位时间，使用秘密武器得到一个彩蛋是瞬间的，得分为当前彩蛋的 $y$ 坐标的千分之一。

现在，Sue 和 Sandy 请你来帮忙，为了满足 Sue 和 Sandy 各自的目标，你决定在收集到所有彩蛋的基础上，得到的分数最高。

输入格式

第一行为两个整数 $N$, $x_0$ 用一个空格分隔，表示彩蛋个数与 Sue 的初始位置。

第二行为 $N$ 个整数 $x_i$，每两个数用一个空格分隔，第 $i$ 个数表示第 $i$ 个彩蛋的初始横坐标。

第三行为 $N$ 个整数 $y_i$，每两个数用一个空格分隔，第 $i$ 个数表示第 $i$ 个彩蛋的初始纵坐标。

第四行为 $N$ 个整数 $v_i$，每两个数用一个空格分隔，第 $i$ 个数表示第 $i$ 个彩蛋匀速沿 $y$ 轴负方向下落的的速度。

输出格式

一个实数，保留三位小数，为收集所有彩蛋的基础上，可以得到最高的分数。

样例 #1

样例输入 #1

3 0
-4 -2 2
22 30 26
1 9 8

样例输出 #1

0.000

提示

对于 $30\%$ 的数据， $N\le 20$。

对于 $60\%$ 的数据， $N\le 100$。

对于 $100\%$ 的数据，$-10^4\le x_i,y_i,v_i\le 10^4$，$N\le 1000$

思路

因为 Sue 肯定会向左或向右去接（你不可能路过不收），所以 Sue 是在区间左端或者右端的。

那么就是一个区间 dp 了，此时我们开三维，f[i][j][0/1] 表示区间 [i,j] 的 Sue此时在左端、右端（0是左端）的最小损耗值，又因为题目求最大值，因此我们得把总数算出来（也就是加上初始 $y$）。为什么是求损耗值呢？因为速度跟损耗值很有关系。

那么我们的状态转移是什么呢？
f l , r , 0 = min ⁡ { f l + 1 , r , 0 + ( x l + 1 − x l ) × ( s 1 , l + s r + 1 , n ) , f l , r , 0 } f l , r , 0 = min ⁡ { f l + 1 , r , 1 + ( x r − x l ) × ( s 1 , l + s r + 1 , n ) , f l , r , 0 } f l , r , 1 = min ⁡ { f l , r − 1 , 0 + ( x r − x l ) × ( s 1 , l − 1 + s r , n ) , f l , r , 0 } f l , r , 1 = min ⁡ { f l , r − 1 , 1 + ( x r − x r − 1 ) × ( s 1 , l − 1 + s r , n ) , f l , r , 0 } f_{l,r,0}=\min\{f_{l+1,r,0}+(x_{l+1}-x_l)\times (s_{1,l}+s_{r+1,n}),f_{l,r,0}\}\\ f_{l,r,0}=\min\{f_{l+1,r,1}+(x_{r}-x_l)\times (s_{1,l}+s_{r+1,n}),f_{l,r,0}\}\\ f_{l,r,1}=\min\{f_{l,r-1,0}+(x_{r}-x_l)\times (s_{1,l-1}+s_{r,n}),f_{l,r,0}\}\\ f_{l,r,1}=\min\{f_{l,r-1,1}+(x_{r}-x_{r-1})\times (s_{1,l-1}+s_{r,n}),f_{l,r,0}\} fl,r,0​=min{fl+1,r,0​+(xl+1​−xl​)×(s1,l​+sr+1,n​),fl,r,0​}fl,r,0​=min{fl+1,r,1​+(xr​−xl​)×(s1,l​+sr+1,n​),fl,r,0​}fl,r,1​=min{fl,r−1,0​+(xr​−xl​)×(s1,l−1​+sr,n​),fl,r,0​}fl,r,1​=min{fl,r−1,1​+(xr​−xr−1​)×(s1,l−1​+sr,n​),fl,r,0​}
注意初始化是 f[x0][x0][0]=f[x0][x0][1]=0，然后len=2开始计算，因为 len=1 我们已经算过了。

代码

//我们可以让 f[0/1][i][j]表示收集区间 [i,j] 的点，0表示在左边收集，1表示在右边收集

//而且Sue只能在边界的地方移动

//要求分数最高，此时就得让损耗最少，此时我们的f就是损耗

#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

#define int long long

const int N = 1e3+10;

int f[N][N][2];
int n,x0;
int ans;
struct E{
    int x,y,v;
    bool operator<(const E& t)const{
        return x<t.x;
    }
}e[N];

signed main(){
    cin>>n>>x0;
    e[n+1].x=x0;
    for(int i=1;i<=n;i++)cin>>e[i].x;
    for(int i=1;i<=n;i++)cin>>e[i].y,ans+=e[i].y;
    for(int i=1;i<=n;i++)cin>>e[i].v;
    n++;
    
    sort(e+1,e+1+n);
    
    memset(f,0x3f,sizeof f);
    
    for(int i=1;i<=n;i++){
        if(e[i].x==x0){
            f[i][i][0]=f[i][i][1]=0;
            break;
        }
    }
    
    for(int i=1;i<=n;i++){
        e[i].v+=e[i-1].v;
    }
    
    for(int len=2;len<=n;len++){
        for(int l=1;l+len-1<=n;l++){
            int r=l+len-1;
            f[l][r][0]=min(f[l+1][r][0]+(e[l+1].x-e[l].x)*(e[l].v+e[n].v-e[r].v),
                            f[l+1][r][1]+(e[r].x-e[l].x)*(e[l].v+e[n].v-e[r].v));
            f[l][r][1]=min(f[l][r-1][1]+(e[r].x-e[r-1].x)*(e[l-1].v+e[n].v-e[r-1].v),
                            f[l][r-1][0]+(e[r].x-e[l].x)*(e[l-1].v+e[n].v-e[r-1].v));
        }
    }
    
    
    printf("%.3lf",(ans-min(f[1][n][0],f[1][n][1]))/1000.0);
    
    return 0;
    
}