[HNOI2008][玩具装箱toy][斜率优化]

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;

int const N= 50010;
int n, L;
typedef long long LL;
LL dp[N], sum[N], A[N], B[N]; 
int que[N], data[N];

/*
dp[i]= min{ dp[j]+ ( s[i]- s[j]+ i- j- 1- L )^2 }
令 A[i]= s[i]- i; B[j]= s[j]+ j+ 1+ L
原方程化为
dp[i]= min{ dp[j]+ ( A[i]- B[j] )* ( A[i]- B[j] ) } 
*/

LL G( int j, int k ){
	return dp[j]- dp[k]+ B[j]* B[j]- B[k]* B[k];
}

LL S( int j, int k ){
	return 2* ( B[j]- B[k] );
}

int main(){
	scanf("%d%d",&n,&L);
	for( int i= 1; i<= n; ++i ){
		scanf("%d", data+ i );
		sum[i]= sum[i-1]+ data[i];
	}
	for( int i= 1; i<= n; ++i ){
		A[i]=  sum[i]+ i;
		B[i]=  sum[i]+ i+ 1+ L;
	}B[0]= 1+ L;
	
	int head= 0, tail= 0; que[0]= 0;
	for( int i= 1; i<= n; ++i ){
		while( head< tail && G( que[head+ 1], que[head] )<= A[i]* S( que[head+ 1], que[head] ) )
		head++;
		
		dp[i]= dp[ que[head] ]+ ( A[i]- B[ que[head] ] )* ( A[i]- B[ que[head] ] );

		while( head< tail && G( que[tail], que[tail- 1] )* S( i, que[tail] )>=
		                     S( que[tail], que[tail- 1] )* G( i, que[tail] ) ) tail--;
		                     
		que[++tail]= i;
	}
	
	cout << dp[n] << endl;
	
	return 0;
}
 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值