洛谷4360 锯木厂选址

跟y_immortal巨佬一起写点斜率优化，就当复习了= =

这个推个柿子得到

其中j<k

维护个下凸壳就行了

初值赋错好几次 以为我写跪了。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 20021225
#define ll long long
#define mxn 20010
#define db double
using namespace std;

int n;
ll suf[mxn],s[mxn],w[mxn],f[mxn];
int que[mxn],hd,tl;

db slope(int x,int y)// x < y
{
	return ((db)suf[y]*s[y]-(db)suf[x]*s[x])/((db)s[y]-s[x]);
}

int main()
{
	scanf("%d",&n); ll tot=0,ans=0;
	for(int i=1;i<=n;i++)
		scanf("%lld%lld",&w[n-i+1],&s[n-i+1]);
	for(int i=1;i<=n;i++)	
		s[i] += s[i-1], tot += w[i]*s[i];
	for(int i=n;i;i--)	suf[i] = suf[i+1] + w[i];
	hd =  1; tl = 1; que[hd]=1; ans=(ll)inf*inf;
	for(int i=2;i<=n;i++)
	{
		while(hd<tl&&slope(que[hd],que[hd+1]) > (db)suf[i])	hd++;
		int x = que[hd];// printf("%d\n",x);
		if(x)	ans = min(ans, tot - suf[x]*s[x] -suf[i] *(s[i]-s[x]));
		while(hd<tl&&slope(que[tl-1],que[tl]) < slope(que[tl],i))	tl--;
		que[++tl] = i;
	}
	printf("%lld\n",ans);
	return 0;
}

 

转载于:https://www.cnblogs.com/hanyuweining/p/10321892.html