bzoj2037 [Sdoi2008]Sue的小球（区间dp）

首先把所有点（包括起点）按x从小到大排序。然后我们肯定是从起点开始贪心的走到左边离我最近的或者右边离我最近的收掉。因此我们可以区间dp，
f[0][i][j]表示取完了i~j的彩蛋，目前在i的最大得分。
f[1][i][j]表示取完了i~j的彩蛋，目前在j的最大得分。
因为我们目前的决策会对以后的所有彩蛋得分也有影响，因此可以把这个影响也计算到当前状态中。即
$f[0][i][j]=y[i]+max\{f[0][i+1][j]+(x_{i+1}-x_i)*w[i+1][j],f[1][i+1][j]+(x_j-x_i)*w[i+1][j]\}$
其中w[i][j]表示除去i~j以外的点的v的和。
时间复杂度$O(n^2)$

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 1010
#define ll long long
ll const inf=1LL<<60;
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,sum[N];ll f[2][N][N];
struct node{
    int x,y,v;
    friend bool operator<(node a,node b){return a.x<b.x;}
}a[N];
int main(){
//  freopen("a.in","r",stdin);
    n=read();int x0=read();
    for(int i=1;i<=n;++i) a[i].x=read();
    for(int i=1;i<=n;++i) a[i].y=read();
    for(int i=1;i<=n;++i) a[i].v=read();
    ++n;a[n].x=x0;a[n].y=0;a[n].v=0;sort(a+1,a+n+1);
    for(int i=1;i<=n;++i) sum[i]=sum[i-1]+a[i].v;
    for(int i=1;i<=n;++i)
        for(int j=i;j<=n;++j)
            if(i==j&&a[i].x==x0&&!a[i].v) f[0][i][j]=f[1][i][j]=0;
            else f[0][i][j]=f[1][i][j]=-inf;
    for(int len=2;len<=n;++len)
        for(int i=1;i<=n;++i){
            int j=i+len-1;if(j>n) break;
            f[0][i][j]=a[i].y+max(f[0][i+1][j]-(ll)(a[i+1].x-a[i].x)*(sum[i]+sum[n]-sum[j]),f[1][i+1][j]-(ll)(a[j].x-a[i].x)*(sum[i]+sum[n]-sum[j]));
            f[1][i][j]=a[j].y+max(f[0][i][j-1]-(ll)(a[j].x-a[i].x)*(sum[i-1]+sum[n]-sum[j-1]),f[1][i][j-1]-(ll)(a[j].x-a[j-1].x)*(sum[i-1]+sum[n]-sum[j-1]));
        }
    printf("%.3lf",max(f[0][1][n],f[1][1][n])*0.001);
    return 0;
}