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



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值