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 longusing 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;}