[ZJOI2007]仓库建设

本文深入探讨了斜率优化动态规划(DP)算法在解决特定类型问题中的应用。通过详细的公式推导和代码实现,展示了如何利用斜率优化来提高DP算法的效率。文章包含关键的数学公式和具体的编程实现细节,适合对算法优化有兴趣的读者。

这题斜率优化DP即可。
一如既往地推公式:
设s1[i]=p[1~i]的和,s2[i]=sigma(p[k]*x[k])(1<=k<=i)

f[i]=min{f[j]+x[i]*(s1[i-1]-s1[j])-s2[i-1]+s2[j]}+c[i]

设k<j<i且j比k更优。

(f[j]+s2[j])-(f[k]+s2[k])/(s1[j]-s1[k])<x[i]

上标:

#include<cstdio>
#define ll long long
#define db double
#define N 1000010
using namespace std;
int n,g[N],len=0,l=0;
ll s1[N],s2[N],f[N],x[N],p[N],c[N];

inline int read()
{
	int x=0; char c=getchar();
	while (c<'0' || c>'9') c=getchar();
	while (c>='0' && c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x;
}

db solve(int x,int y) {return (db)(f[x]+s2[x]-f[y]-s2[y])/(s1[x]-s1[y]);}

int main()
{
	freopen("P2120.in","r",stdin);
	freopen("P2120.out","w",stdout);
//	f[i]=min{f[j]+x[i]*(s1[i-1]-s1[j])-s2[i-1]+s2[j]}+c[i]
//	(f[j]+s2[j])-(f[k]+s2[k])/(s1[j]-s1[k])<x[i]
	n=read();
	for (int i=1;i<=n;i++)
	{
		x[i]=read(),p[i]=read(),c[i]=read();
		s1[i]=s1[i-1]+p[i],s2[i]=s2[i-1]+p[i]*x[i];
	}
	for (int i=1;i<=n;i++)
	{
		while (l<len && solve(g[l+1],g[l])<x[i]) l++;
		f[i]=f[g[l]]+x[i]*(s1[i-1]-s1[g[l]])-s2[i-1]+s2[g[l]]+c[i];
		while (solve(g[len],g[len-1])>=solve(i,g[len]) && len) len--;
		g[++len]=i;
	}
	printf("%lld\n",f[n]);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值