3437: 小P的牧场

斜率优化什么的真心无脑。

令sum[i]=sigma(b[k])(k<=i) mul[i]=sigima(k*b[k])(k<=i)。

则f[i]=min(f[j]+i*(sum[i]-sum[j])-mul[i]+mul[j])(j<i)

单调队列维护就好了。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=1000000+5;
typedef long long ll;
ll f[N],sum[N],mul[N],a[N],b;
int q[N];
double slop(int j,int k){
	return (double)(f[j]+mul[j]-f[k]-mul[k])/(double)(sum[j]-sum[k]);
}
int main(){
	int n;scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
	for(int i=1;i<=n;i++){
		scanf("%lld",&b);
		sum[i]=sum[i-1]+b;
		mul[i]=mul[i-1]+b*i;
	}
	int h=1,t=1;q[1]=0;
	for(int i=1;i<=n;i++){
		while(h<t&&slop(q[h+1],q[h])<i)h++;
		int j=q[h];
		f[i]=f[j]+i*(sum[i]-sum[j])-mul[i]+mul[j]+a[i];
		while(h<t&&slop(q[t-1],q[t])>slop(q[t],i))t--;
		q[++t]=i;
	}
	printf("%lld",f[n]);
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值