BZOJ 3675: [Apio2014]序列分割

转载于 2018-10-19 19:24:00 发布 · 58 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/silenty/p/9818329.html

本文深入探讨了斜率优化在动态规划中的应用，通过具体代码示例解释了如何利用斜率优化来提高算法效率，特别是在处理顺序无关的问题时。文章详细讲解了斜率优化的数据结构实现，以及如何通过比较斜率进行动态规划状态转移。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

难度在于认识到顺序是没有关系的

斜率优化

#include<cstdio>
#include<algorithm>
using namespace std;
int n,m,K,Sum[100005],a[100005];
long long F[2][100005];
struct node{
	long long k,b;
}stack[1000005];
double check(node a,node b){
	return ((double)a.b-b.b)/(b.k-a.k);
}
int main(){
	scanf("%d%d",&m,&K);
	for (int i=1; i<=m; i++) {
		int x;
		scanf("%d",&x);
		if (x) a[++n]=x;
	}
	for (int i=1; i<=n; i++) Sum[i]=Sum[i-1]+a[i];
	for (int k=1; k<=K; k++){
		for (int i=0; i<=n; i++) F[k%2][i]=0;
		int head=1,tail=0;
		for (int i=0; i<=n; i++){
			node line=(node){Sum[i],F[(k-1)%2][i]-1ll*Sum[i]*Sum[i]};
			while (tail-head>0 && check(stack[tail],stack[tail-1])>check(stack[tail],line)) tail--;
			stack[++tail]=line;
			while (tail-head>0 && check(stack[head],stack[head+1])<Sum[i]) head++;
			F[k%2][i]=stack[head].k*Sum[i]+stack[head].b;
		}
	}
	printf("%lld\n",F[K%2][n]);
	return 0;
}

转载于:https://www.cnblogs.com/silenty/p/9818329.html