bzoj 2037: [Sdoi2008]Sue的小球

题意：

彩蛋有一个魅力值，这个魅力值会随着彩蛋在空中降落的时间而降低，Sue要想得到更多的分数，必须尽量在魅力值高的时候收集这个彩蛋，而如果一个彩蛋掉入海中，它的魅力值将会变成一个负数，Sue的初始位置为(x0, 0)，Sue可以沿x轴的正方向或负方向移动，Sue的移动速度是1单位距离/单位时间，使用秘密武器得到一个彩蛋是瞬间的，得分为当前彩蛋的y坐标的千分之一。收集到所有彩蛋的基础上，得到的分数最高。

题解：

直接dp就可以了。
$f[i][j][k]$表示取完了第i个球到第j个球，且站在左边/又边的最小代价。
最后用总得分减去就可以了。
我个sb这种题wa三次。
code：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
struct node{
    int x,y,v;
}a[1010];
int f[1010][1010][3];
int n,X;
int lsum[1010],rsum[1010];
bool cmp(node a,node b){return a.x<b.x;}
void cmin(int &x,int y){x=min(x,y);}
int main()
{
    scanf("%d %d",&n,&X);
    for(int i=1;i<=n;i++) scanf("%d",&a[i].x);
    for(int i=1;i<=n;i++) scanf("%d",&a[i].y);
    for(int i=1;i<=n;i++) scanf("%d",&a[i].v);
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++) lsum[i]=lsum[i-1]+a[i].v;
    for(int i=n;i>=1;i--) rsum[i]=rsum[i+1]+a[i].v;
    int t=n+1;
    for(int i=1;i<=n;i++) if(a[i].x>=X){t=i;break;}
    memset(f,127,sizeof(f));
    if(t!=n+1&&X==a[t].x) f[t][t][0]=f[t][t][1]=0;
    else
    {
        if(t!=1) f[t-1][t-1][0]=(X-a[t-1].x)*lsum[n];
        if(t!=n+1) f[t][t][1]=(a[t].x-X)*lsum[n];
    }
    for(int i=t;i>=1;i--)
        for(int j=i;j<=n;j++)
        for(int k=0;k<=1;k++)
            {
                if(f[i][j][k]>1000000000) continue;
                if(k==0)
                {
                    if(i!=1) cmin(f[i-1][j][0],f[i][j][k]+(a[i].x-a[i-1].x)*(lsum[i-1]+rsum[j+1]));
                    if(j!=n) cmin(f[i][j+1][1],f[i][j][k]+(a[j+1].x-a[i].x)*(lsum[i-1]+rsum[j+1]));
                }
                if(k==1)
                {
                    if(i!=1) cmin(f[i-1][j][0],f[i][j][k]+(a[j].x-a[i-1].x)*(lsum[i-1]+rsum[j+1]));
                    if(j!=n) cmin(f[i][j+1][1],f[i][j][k]+(a[j+1].x-a[j].x)*(lsum[i-1]+rsum[j+1]));
                }
            }
    double ans=0;
    for(int i=1;i<=n;i++) ans+=(double)a[i].y;
    double ans1=(double)min(f[1][n][0],f[1][n][1]);
    printf("%.3lf",(ans-ans1)/1000.0);
}