大致意思就是现在你要不断的奔跑到不同的地点去接球，每一秒可以移动一个单位长度，而你接到一个球的动作是瞬间的，收益是y[i]-t*v[i] 然后呢，要求分数最高。
起初看这个题目QWQ完全没有任何思路,大概只能想到......
先按照x排序（记得把起始位置也加进去）
然后令f[l][r]表示收集完l~r的球，最后在l的最大收益
g[l][r]收集完l~r的球，最后在r的最大收益
然后...然后....然后....
我就去看题解了。
好了 进入正题。
首先我们定义
f[l][r]表示收集完l~r的球，最后在l的最小损失
g[l][r]收集完l~r的球，最后在r的最小损失
最后用总收益减去损失
在按照x排完序之后
进行区间dp，由小区间转到大区间
f[l][r]可以从f[l+1][r]和g[l+1][r]转移而来
g[l][r]可以从f[l][r-1]和g[l][r-1]转移而来
我们可以这么理解
每当我们去接下一个球的时候，其他球在向下掉，相当于我们损失了这些的收益
那么时间就是x之差的绝对值，然后提前用前缀和预处理v
就可以直接算出损失了多少收益了
f[l][r]=min(f[l][r],f[l+1][r]+(sum[n]-sum[r]+sum[l])abs(a[l+1].x-a[l].x));
f[l][r]=min(f[l][r],g[l+1][r]+(sum[n]-sum[r]+sum[l])abs(a[r].x-a[l].x));
g[l][r]=min(g[l][r],f[l][r-1]+(sum[n]-sum[r-1]+sum[l-1])*abs(a[l].x-a[r].x));
g[l][r]=min(g[l][r],g[l][r-1]+(sum[n]-sum[r-1]+sum[l-1])*abs(a[r-1].x-a[r].x));
转移式子就不过多解释了
然后最后用ans-min(f[1][n],g[1][n])再 /1000就行
最后注意初始化的时候 嗯
QWQ我的写法和很多题解都不一样 不过也过了QWQ不太知道是为什么
for (int i=1;i<=n;i++) f[i][i]=abs(a[i].x-start)sum[n],g[i][i]=abs(a[i].x-start)sum[n];
上代码 嗯
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 1010;
struct Node{
int v,x,y;
};
Node a[maxn];
int f[maxn][maxn]; //i~j接完 最后在i
int g[maxn][maxn]; // i~j接完，最后在j
int sum[maxn];
int start,n;
double ans;
bool cmp(Node a,Node b)
{
return a.x<b.x;
}
int main()
{
n=read();
start=read();
memset(f,127/3,sizeof(f));
memset(g,127/3,sizeof(g));
for (int i=1;i<=n;i++) a[i].x=read();
for (int i=1;i<=n;i++) a[i].y=read(),ans+=a[i].y;
for (int i=1;i<=n;i++) a[i].v=read();
n++;
a[n].x=start;
sort(a+1,a+1+n,cmp);
for (int i=1;i<=n;i++) sum[i]=sum[i-1]+a[i].v;
for (int i=1;i<=n;i++) f[i][i]=abs(a[i].x-start)*sum[n],g[i][i]=abs(a[i].x-start)*sum[n];
for (int i=2;i<=n;i++)
for (int l=1;l<=n-i+1;l++)
{
int r = l+i-1;
f[l][r]=min(f[l][r],f[l+1][r]+(sum[n]-sum[r]+sum[l])*abs(a[l+1].x-a[l].x));
f[l][r]=min(f[l][r],g[l+1][r]+(sum[n]-sum[r]+sum[l])*abs(a[r].x-a[l].x));
g[l][r]=min(g[l][r],f[l][r-1]+(sum[n]-sum[r-1]+sum[l-1])*abs(a[l].x-a[r].x));
g[l][r]=min(g[l][r],g[l][r-1]+(sum[n]-sum[r-1]+sum[l-1])*abs(a[r-1].x-a[r].x));
}
ans=ans-min((double)f[1][n],(double)g[1][n]);
printf("%.3lf",ans/1000);
return 0;
}