bzoj1358: [Baltic2009]Beetle

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=1358

题解

　　AC的人这么少？
　　远程膜$Claris$。
　　$Claris$说你可以枚举最后走了几个点，然后就成了裸的费用提前计算DP。每次的损失就可以直接得到了。
　　$f[i][j][0/1]$当前已经喝了i到j的所有水，现在在左/右端点的最小损失。
　　因为枚举了走几个点，所以每次转移的损失可以直接算出来。
　　不用担心水会减到零的问题，因为如果减成了负数，这个答案肯定会被一种没减成负数的答案更新掉。
　　如上。

代码

//动态规划
#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 350
#define inf ((long long)(1)<<40)
#define ll long long
using namespace std;
ll f[maxn][maxn][2], cnt, start, x[maxn], tmp[maxn], INIT, N, M;
void init()
{
    int i, t=1;
    scanf("%lld%lld",&N,&M);
    INIT=0;
    for(i=1;i<=N;i++)scanf("%lld",x+i),tmp[i]=x[i],INIT=x[i]==0?M:INIT;
    if(!INIT)N++,tmp[N]=x[N]=0;
    sort(x+1,x+N+1);
    for(i=2;i<=N;i++)if(x[i]!=x[i-1])x[++t]=x[i];
    for(i=1;i<=N;i++)if(x[i]==0)start=i;
}
inline void upd(ll &x, ll &y, ll t){y=min(y,x+t);}
ll dp()
{
    ll i, j, l, ans=inf;
    for(l=1;l<=cnt;l++)for(i=1;i<=N-l+1;i++)f[i][i+l-1][0]=f[i][i+l-1][1]=inf;
    f[start][start][0]=f[start][start][1]=0;
    for(l=1;l<cnt;l++)
        for(i=1;i<=N-l+1;i++)
        {
            j=i+l-1;
            upd(f[i][j][0],f[i][j+1][1],(x[j+1]-x[i])*(cnt-l));
            upd(f[i][j][0],f[i-1][j][0],(x[i]-x[i-1])*(cnt-l));
            upd(f[i][j][1],f[i][j+1][1],(x[j+1]-x[j])*(cnt-l));
            upd(f[i][j][1],f[i-1][j][0],(x[j]-x[i-1])*(cnt-l));
        }
    for(i=1;i<=N-cnt+1;i++)ans=min(ans,min(f[i][i+cnt-1][0],f[i][i+cnt-1][1]));
    return ans;
}
int main()
{
    ll ans=-inf;
    init();
    for(cnt=1;cnt<=N;cnt++)ans=max(ans,M*cnt-dp());
    if(!INIT)ans-=M;
    printf("%lld",ans);
    return 0;
}