BZOJ 1010: [HNOI2008]玩具装箱toy——斜率优化

BZOJ 1010


第二道斜率优化,对斜率优化有了新的理解。开心(大雾


这题的题面其实就是在说  斜率优化!斜率优化!斜率优化!

好吧就是道裸题



我们考虑更新f[i]时选择的节点


如果k比j更优


则有:

f[j]+(j-i+sum[i]-sum[j]-L)^2

f[j]+(sum[i]-sum[j])^2-2*(i-j+L)*(sum[i]-sum[j])+(i-j+L)^2

令t1=i-j+L,t2=i-k+L;

f[j]+(sum[i]-sum[j])^2-2*t1*(sum[i]-sum[j])+t1^2

f[j]-f[k]+(sum[i]-sum[j])^2-(sum[i]-sum[k])^2+t1^2-t2^2<2*sum[i]*(sum[k]-sum[j]+k-j)


f[j]-f[k]+(sum[i]-sum[j])^2-(sum[i]-sum[k])^2+t1^2-t2^2)/(2*sum[k]-sum[j]+k-j)


。。别看这个公式很烦,其实你如果用Sum[i]去表示sum[i]+i的话很方便,只是这里我写的严谨了点


嗯、、Sum[i]表示的话,公式如下

(f[j]-f[k]+(sum[j]+l)^2+(sum[k]+l)^2)/(2*(sum[k]-sum[j])


好了就是这样


考虑这个东西是个斜率式子,我们要求最小值,所以斜率不可能递减


为什么?因为递减的话你还不如直接从上一个位置更新过来


所以我们可以根据这个性质,维护一个下凸壳


顺便一说,如果是要最大值,则要维护上凸壳


补:

突然发现忘记放代码了。。。。

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define inf 1e9
#define eps 1e-9
#define iter multiset::iterator 
#define ll long long
#define maxn 20010
#define For(i,j,k) for(ll i=j;i<=k;i++)
#define Dow(i,j,k) for(ll i=k;i>=j;i--)
using namespace std;
inline void read(ll &tx){   ll x=0,f=1;char ch=getchar();   while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}  while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}  tx=x*f; }
inline void write(ll x){    if (x<0) putchar('-'),x=-x; if (x>=10) write(x/10);   putchar(x%10+'0');  }
inline void writeln(ll x){write(x);puts("");}
ll n,l,a[100001],q[100001],L,R,f[100001];
inline ll sqr(ll x){return x*x;}
inline double get(ll y,ll x)
{
	return (double)(f[x]-f[y]+sqr(a[x]+l)-sqr(a[y]+l))/(2.0*(a[x]-a[y]));
}
int main()
{
	read(n);read(l);l++;
	For(i,1,n)	read(a[i]);
	For(i,1,n)	a[i]+=a[i-1]+1;
	L=R=1;q[1]=0;
	For(i,1,n)
	{
		while(Lget(q[R],i))	R--;
		q[++R]=i;
	}
	writeln(f[n]);
}	




你可能感兴趣的:(题库,bzoj,dp)