【CEOI2004】锯木厂选址(斜率优化入门题)

#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<iostream>
#include<cstring>
#include<algorithm> 
#define LL long long
using namespace std;
LL n,ans=999999999;
LL w[2001000],dis[2001000],sw[2001000],val[2001000],d[2001000];
LL q[2001000];
double Slpoe(int j,int k)
{
	return (dis[j]*val[j]-dis[k]*val[k])*1.0/(val[j]-val[k]);
} 
int main()
{
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%lld%lld",&w[i],&d[i]);
		val[i]=val[i-1]+w[i];
		dis[i]=dis[i-1]+d[i-1];
		sw[i]=sw[i-1]+d[i-1]*val[i-1];// cout<<sw[i]<<" ";
	}
	val[n+1]=val[n];
	dis[n+1]=dis[n]+d[n];
	sw[n+1]=sw[n]+val[n]*d[n];
//	for(int i=1;i<=n+1;i++) cout<<dis[i]<<' '<<val[i]<<' '<<sw[i]<<endl; 
	int l=1,r=1;
	for(int i=1;i<=n;i++)
	{//cout<<q[l]<<' '<<q[l+1];
		while(l<r&&Slpoe(q[l],q[l+1])<dis[i])l++;//l+1比l更优,向上走 
		LL j=q[l];
		ans=min(ans,sw[n+1]-val[j]*(dis[i]-dis[j])-val[i]*(dis[n+1]-dis[i]));
		while(l<r&&Slpoe(q[r-1],q[r])>Slpoe(q[r],i)) r--;
		q[++r]=i; //cout<<i<<' '<<j<<endl;
	} 
	printf("%lld\n",ans);
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值