BZOJ2037 Sue的小球 【区间dp】

WOJ2039

---

SOL

区间dp得好好理解一下了。。

• 我们注意到$y$坐标的限制是不存在的。$minimize:Delta=\sum t[i]\ast v[i]$

• 考虑子问题，按$x$坐标排序后，一段点被走完后只会停在最左端或者最右端，再“扩展段长”。

• 不放设：$f[0/1][i][j]，\text{表示i到j走完，停在左/右端点,总共的delta值}$。

• 区间dp转移即可。注意到，这里记录了左右端点，转移的时候不需要"折返跑"（即$f[0][i][j]$由$f[0/1][i][j+1]$转移的情况）。

---

CODE

#include<bits/stdc++.h>
using namespace std;
#define sf scanf
#define pf printf
#define ll long long
#define cs const
#define db double
#define ri register int
#define gc getchar()
#define in red()
inline int red(){
	int num=0,f=1;char c=gc;
	for(;!isdigit(c);c=gc)if(c=='-')f=-1;
	for(;isdigit(c);c=gc)num=(num<<1)+(num<<3)+(c^48);
	return num*f;
}
cs int N=1100;
db sum[N],x[N],v[N],tot,f[2][N][N],x0;
int ord[N];
inline bool cmp(int a,int b){return x[a]<x[b];}
int n;
signed main(){
	n=in;x0=in;
	for(ri i=1;i<=n;++i)x[i]=in,ord[i]=i;
	for(ri i=1;i<=n;++i)tot+=(in);
	for(ri i=1;i<=n;++i)v[i]=in;
	sort(ord+1,ord+n+1,cmp);
	for(ri i=1;i<=n;++i)sum[i]=sum[i-1]+v[ord[i]];
	sort(x+1,x+n+1);
	for(ri i=1;i<=n;++i)for(ri j=1;j<=n;++j)for(ri k=0;k<2;++k)f[k][i][j]=1e20;
	int id=lower_bound(x+1,x+n+1,x0)-x;
	f[0][id][id]=f[1][id][id]=sum[n]*(x[id]-x0);
	--id;
	f[0][id][id]=f[1][id][id]=sum[n]*(x0-x[id]);
	for(ri k=2;k<=n;++k){
		for(ri i=1;i<=n;++i){
			int j=i+k-1;
			if(j>n)break;
			f[0][i][j]=min(f[0][i][j],f[0][i+1][j]+(sum[n]-sum[j]+sum[i])*(x[i+1]-x[i]));
			f[0][i][j]=min(f[0][i][j],f[1][i+1][j]+(sum[n]-sum[j]+sum[i])*(x[j]-x[i]));
			f[1][i][j]=min(f[1][i][j],f[0][i][j-1]+(sum[n]-sum[j-1]+sum[i-1])*(x[j]-x[i]));
			f[1][i][j]=min(f[1][i][j],f[1][i][j-1]+(sum[n]-sum[j-1]+sum[i-1])*(x[j]-x[j-1]));
		}
	}
	printf("%.3lf",(tot-min(f[0][1][n],f[1][1][n]))/1000.0);

	return 0;
}