【BZOJ1010】【HNOI2008】玩具装箱toy 动规_斜率优化

第一道斜率优化题，题解仍然不是很透彻。

我在这里可以推荐一篇博客。

http://blog.sina.com.cn/s/blog_5f5353cc0100jx41.html

写得还好，但是可能有一些地方仍说的不是很明白。

我说一下我对斜率优化的理解：

    就是发现某道动态规划的题会TLE，然后它又恰巧符合某些丧心病狂的性质时的一种优化。

    本质是有限度地捏造出一种性质，使得动规方程满足某种邪恶的单调性，从而得到大幅度的优化。

    

    而这个单调性是怎么得到的呢？

    就是我们把原转移状态搞出来，一顿转化，比如把平方拆开等，本来有两种转移j->i和k->i，我们使它原来有一个大小关系DP[j->i]<=DP[k->i]，而最终转化成左边是一个有j和k的式子，右边是一个只有i的式子，这样达到高速判断两个转移哪家强的目的。

    而思想上wyf给我说了一个好理解的版本如下。

    就是它！我们要维护一个凸壳，当i要加入优先队列判断j->i的斜率和k->i的斜率，然后显然k出去了，再判断h->i,嗯，留下来了，就是这样。

    好了，说一些跟本题相关的具体东西。

    就是首先定义f[i]=sum[i]+i,这个sum[i]是前缀和。

然后有当G（j,k）/S(j,k)<=i时j比k优。

#define G(x,y) (dp[x]-dp[y]+(f[x]+c)*(f[x]+c)-(f[y]+c)*(f[y]+c))
#define S(x,y) (2*(f[x]-f[y]))

贴代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 55000
#define G(x,y) (dp[x]-dp[y]+(f[x]+c)*(f[x]+c)-(f[y]+c)*(f[y]+c))
#define S(x,y) (2*(f[x]-f[y]))
using namespace std;
long long c,sum[N],f[N],dp[N];
int n,q[N];
int main()
{
//	freopen("test.in","r",stdin);
	scanf("%d%lld",&n,&c);c++;
	for(int i=1;i<=n;i++)
	{
		scanf("%lld",&sum[i]);
		sum[i]+=sum[i-1];
		f[i]=sum[i]+i;
	}
	int head=1,tail=1;
	q[1]=0;
	for(int i=1;i<=n;i++)
	{
		while(head<tail&&G(q[head+1],q[head])<=S(q[head+1],q[head])*f[i]) head++;
		dp[i]=dp[q[head]]+(f[i]-f[q[head]]-c)*(f[i]-f[q[head]]-c);
		while(head<tail&&G(i,q[tail])*S(q[tail],q[tail-1])<=G(q[tail],q[tail-1])*S(i,q[tail]))tail--;
		q[++tail]=i;
	}
	printf("%lld\n",dp[n]);
	return 0;
}