【整理】斜率or单调队列优化dp

最新推荐文章于 2024-07-12 22:24:20 发布
最新推荐文章于 2024-07-12 22:24:20 发布
阅读量502 收藏
点赞数
分类专栏： 动态规划 分类题集 文章标签： 斜率优化dp 单调队列优化dp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_21841245/article/details/44422795
版权

【1】HDU2993 MAX Average Problem

题意：求一个序列的子区间满足长度大于k且所有数平均值最大

周源论文里的题。。之前有人说周源讲的是错的 其实应该是没什么问题的 可以O(n)求出
不过这题hdu上的数据不知道怎么了 = = 反正我在网上找的ac代码们全都TLE。。。。
总之……意思明白就好 反正也是入门题…… 

#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <cmath> 

using namespace std;

int read()
{
	int sign = 1, n = 0; char c = getchar();
	while(c < '0' || c > '9'){ if(c == '-') sign = -1; c = getchar(); }
	while(c >= '0' && c <= '9') { n = n*10 + c-'0'; c = getchar(); }
	return sign*n;
}

int N, K;
int a[100005];
double sum[100005];

namespace queue{
	int q[100005], head, tail;
	
	inline void init(){ head = 1; tail = 0; }
	inline void push_back(int x){ q[++tail] = x; }
	inline void pop_back() { --tail; }
	inline void pop_front() { ++head; }
	inline int front() { return q[head]; }
	inline int second() { return q[head + 1]; }
	inline int back() { return q[tail]; }
	inline int before() { return q[tail - 1]; }
	inline int size() { return tail - head + 1; } 
}

inline double get_k(int i, int j){ return (sum[i] - sum[j]) / (i - j); }

int main()
{
	while(~scanf("%d%d", &N, &K))
	{
		for(int i = 1; i <= N; ++i)
		{
			a[i] = read();
			sum[i] = sum[i - 1] + a[i];
		}
		
		using namespace queue;
		double ans = 1e-100; init();
		
		for(int i = K; i <= N; ++i)
		{
			int now = i - K;
			while (size() > 1 && get_k(back(), before()) >= get_k(now, back()) ) pop_back();
			push_back(now);
			while (size() > 1 && get_k(front(), i) <= get_k(second(), i)) pop_front();
			ans = fmax(ans, get_k(front(), i));
		}
		
		printf("%.2f\n", ans);
	}
	return 0;
}

【2】HDU3480 Division


题意：把一个集合 分成若干个子集合 每个集合的cost为最大值与最小值差的平方 求最大cost的最小值

和之前摩天轮那道题是一样的。。先排个序 很明显每个集合是要取连续的一段区间
那么dp[i][j]表示前i个人分成了j个区间的最大cost的最小值。。于是dp[i][j] = max{ dp[k-1][j-1] + (sum[i] - sum[k]) ^ 2 }
然后展开以后就可以斜率优化了……还有就是第二维可以滚动。。

说实话这题写的我略痛苦。。因为还是不太熟TAT 

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

int read()
{
	int sign = 1, n = 0; char c = getchar();
	while(c < '0' || c > '9'){ if(c == '-') sign = -1; c = getchar(); }
	while(c >= '0' && c <= '9') { n = n*10 + c-'0'; c = getchar(); }
	return sign*n;
}

const int Nmax = 10005;

int N, M;
int a[Nmax], sqr[Nmax]; 
int dp[Nmax][2], G[Nmax];

namespace queue{
	int q[100005], head, tail;
	
	inline void init(){ head = 1; tail = 0; }
	inline void push_back(int x){ q[++tail] = x; }
	inline void pop_back() { --tail; }
	inline void pop_front() { ++head; }
	inline int front() { return q[head]; }
	inline int second() { return q[head + 1]; }
	inline int back() { return q[tail]; }
	inline int before() { return q[tail - 1]; }
	inline int size() { return tail - head + 1; } 
}

inline double getk(int i, int j)
{
	if(a[i] == a[j]) return 1e200;
	return 1. * (G[i] - G[j]) / (a[i] - a[j]);
}

int main()
{
	for(int T = read(), cas = 0; T--; )
	{
		N = read(), M = read();
		for(int i = 1; i <= N; ++i)
		{
			a[i] = read();
			dp[i][0] = 0x3f3f3f3f;
		} 
		
		sort(a + 1, a + N + 1);
		for(int i = 1; i <= N; ++i) sqr[i] = a[i] * a[i];
		
		using namespace queue;
		for(int j = 1; j <= M; ++j)
		{
			int nw = j & 1, pr = nw ^ 1;
			init(); 
			for(int i = j; i <= N; ++i)
			{
				G[i] = dp[i - 1][pr] + sqr[i];
				while(size() > 1 && getk(before(), back()) >= getk(back(), i)) pop_back();
				push_back(i);
				
				while(size() > 1 && getk(front(), second()) <= 2. * a[i]) pop_front();
				dp[i][nw] = G[front()] + sqr[i] - 2 * a[i] * a[front()];
			}
		}
		
		printf("Case %d: %d\n", ++cas, dp[N][M & 1]);
	}
	return 0;
}


【3】HDU3530 Subsequence


题意：找出一个数列中的最长的一个子区间 满足最大值和最小值之差在[M, K]之间

开两个单调队列一个维护递增一个维护递减就可以了。。然后每次判断队首合不合法。。。并更新答案。。。
这题我最开始队列初始化为a[1]就一直会挂 = =我也不知道为什么 委屈 

#include <cstdio>
#include <iostream>
#include <cstdlib>

using namespace std;

int read()
{
	int sign = 1, n = 0; char c = getchar();
	while(c < '0' || c > '9'){ if(c == '-') sign = -1; c = getchar(); }
	while(c >= '0' && c <= '9') { n = n*10 + c-'0'; c = getchar(); }
	return sign*n;
}

int N, M, K;
int a[100005];

struct queue{
	int Q[100005], head, tail;
	
	inline void init(){ head = 1; tail = 0; }
	inline void push_back(int x){ Q[++tail] = x; }
	inline void pop_back() { --tail; }
	inline void pop_front() { ++head; }
	inline int front() { return Q[head]; }
	inline int back() { return Q[tail]; }
	inline int size() { return tail - head + 1; } 
}q[2];

int main()
{
	while(~scanf("%d%d%d", &N, &M, &K))
	{
		for(int i = 1; i <= N; ++i) a[i] = read();
		
		int left = 0, ans = 0; q[0].init(); q[1].init();
		for(int i = 1; i <= N; ++i)
		{
			while(q[0].size() && a[q[0].back()] <= a[i]) q[0].pop_back(); 
			q[0].push_back(i);
			while(q[1].size() && a[q[1].back()] >= a[i]) q[1].pop_back(); 
			q[1].push_back(i);
			
			while(a[q[0].front()] - a[q[1].front()] > K)
			{
				if(q[0].front() < q[1].front()){ left = q[0].front(); q[0].pop_front(); } 
				else { left = q[1].front(); q[1].pop_front(); }
			}
			
			if(a[q[0].front()] - a[q[1].front()] >= M) ans = max(ans, i - left);
		}
		printf("%d\n", ans);
	} 
	
	return 0;
}

【4】HYSBZ1911 特别行动队


题意：……中文题就不写了……

然后这题的公式还是蛮好推的。。但是因为特别长所以要注意细节。。。
A是小于0的 我们要求的是最大值 那么就维护一个上凸包（其实也可以转化成下凸包。。都一样）
f[i] = max( f[j] + a*(s[i]-s[j])^2 + b*(s[i]-s[j]) + c) 
     = ( f[j] + a*s[j]^2  - b*s[j] ) + ( a*s[i]^2 +b*s[i]+c )  - ( 2*a*s[i]*s[j] )
分三块。。很明显是不是。。 

#include <cstdio>
#include <iostream>

using namespace std;

typedef long long LL;

const int Nmax = 1000005; 
int N;
LL a[Nmax], sum[Nmax], sqr[Nmax], A, B, C;
LL dp[Nmax], F[Nmax];

namespace queue{
	int q[Nmax], head, tail;
	
	inline void init(){ head = 1; tail = 0; }
	inline void push_back(int x){ q[++tail] = x; }
	inline void pop_back() { --tail; }
	inline void pop_front() { ++head; }
	inline int front() { return q[head]; }
	inline int second() { return q[head + 1]; }
	inline int back() { return q[tail]; }
	inline int before() { return q[tail - 1]; }
	inline int size() { return tail - head + 1; } 
}

double get_k(int i, int j)
{
	if(sum[i] == sum[j]) return 1e200;
	return 1. * (F[i] - F[j]) / (sum[i] - sum[j]);
}

int main()
{
	ios::sync_with_stdio(false);
	cin >> N >> A >> B >> C;
	for(int i = 1; i <= N; ++i)
	{
		cin >> a[i];
		sum[i] = sum[i - 1] + a[i];
		sqr[i] = sum[i] * sum[i];
		dp[i] = A * sqr[i] + B * sum[i] + C;
	} 
	
	using namespace queue; init(); 
	for(int i = 1; i <= N; ++i)
	{
		int K = 2 * A * sum[i];
		while(size() > 1 && get_k(front(), second()) >= K) pop_front();
		dp[i] = max(dp[i], F[front()] + dp[i] - 2 * A * sum[front()] * sum[i]);
		F[i] = dp[i] + A * sqr[i] - B * sum[i];
		while(size() > 1 && get_k(before(), back()) <= get_k(back(), i)) pop_back();
		push_back(i);
	}
	
	cout << dp[N] << endl;
	
	return 0;
}

【5】HDU 3507 Print Article


题意：把N个字母分成若干行 每个字母有个ci 然后每行的cost为sigma(ci) + M 求总cost的最小值

和上面那题差不多 而且公式还简单些 = = 

#include <cstdio>
#include <iostream>

using namespace std;

typedef long long LL;

const int Nmax = 500005; 
int N, M;
LL a[Nmax], sum[Nmax], sqr[Nmax];
LL dp[Nmax], F[Nmax];

namespace queue{
	int q[Nmax], head, tail;
	
	inline void init(){ head = 1; tail = 0; q[head] = 0;}
	inline void push_back(int x){ q[++tail] = x; }
	inline void pop_back() { --tail; }
	inline void pop_front() { ++head; }
	inline int front() { return q[head]; }
	inline int second() { return q[head + 1]; }
	inline int back() { return q[tail]; }
	inline int before() { return q[tail - 1]; }
	inline int size() { return tail - head + 1; } 
}

double get_k(int i, int j)
{
	if(sum[i] == sum[j]) return 1e200;
	return 1. * (F[i] - F[j]) / (sum[i] - sum[j]);
}

int main()
{
	ios::sync_with_stdio(false);
	while(cin >> N >> M)
	{
		for(int i = 1; i <= N; ++i)
		{
			cin >> a[i];
			sum[i] = sum[i - 1] + a[i];
			sqr[i] = sum[i] * sum[i];
			dp[i] = sqr[i] + M;
		} 
		
		using namespace queue; init(); 
		for(int i = 1; i <= N; ++i)
		{
			int K = 2 * sum[i];
			while(size() > 1 && get_k(front(), second()) <= K) pop_front();
			dp[i] = min(dp[i], F[front()] + dp[i] - 2 * sum[front()] * sum[i]);
			F[i] = dp[i] + sqr[i];
			while(size() > 1 && get_k(before(), back()) >= get_k(back(), i)) pop_back();
			push_back(i);
		}
		
		cout << dp[N] << endl;
	}

	return 0;
}

【6】HDU 3669 Cross the Wall


题解：http://blog.csdn.net/qq_21841245/article/details/44461527


MoeO3
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
7.29-7.30集训——dp优化
DL24zhanghao的博客
07-29 189
dp优化 决策单调性 四边形不等式 一种特殊的1D/1D决策单调性 斜率优化 wqs二分
斜率优化DP
小蒟蒻yyb的博客
10-21 1128
斜率优化DP题外话考试的时候被这个玩意弄得瑟瑟发抖 大概是Day4(http://www.cnblogs.com/cjyyb/p/7705308.html) 小蒟蒻表示根本不会做….. 然后自己默默地搞了一下斜率优化算是开始吗??其实我讲的会非常非常非常简单,,,而且绝对没有一张图(因为我绘图水平太菜) 貌似没太多友善的题目可以用来搞….算了 虚一点,缥缈一点的来说吧….其实我就是写给自己
Hdu 3401 题解 单调队列优化DP
Hawo11的博客
05-06 489
Hdu3401
POJ 1180 斜率优化DP（单调队列）
weixin_34055910的博客
04-26 93
Batch Scheduling Time Limit:1000MS Memory Limit:10000K Total Submissions:4347 Accepted:1992 Description There is a sequence of N jobs to be processed on one machine...
【单调队列优化】CF319C——Kalila and Dimna in the Logging Industry
十里寒塘路♣烟花一半醒
12-27 505
前言 本来挺简单的一道题... 我发现自己不太擅长简化、抽象模型，老是不知道怎么把学过的技巧用上去_(:зゝ∠)_ 题目 Kalila and Dimna are two jackals living in a huge jungle. One day they decided to join a logging factory in order to make money. The m...
POJ1180 [IOI2002]任务安排（斜率优化dp入门好题）
The_OIer的博客
11-12 2830
题目描述 There is a sequence of N jobs to be processed on one machine. The jobs are numbered from 1 to N, so that the sequence is 1,2,…, N. The sequence of jobs must be partitioned into one or more batche...
poj 1180 Batch Scheduling(斜率优化dp或决策单调)
madaidao的专栏
11-05 1037
题目链接
DP优化技巧-斜率优化（基础版）
最新发布
Never_gonna的博客
07-12 347
DP优化技巧-斜率优化（基础版） 本文尚未触及到用平衡树、CDQ分治与斜率优化结合，这将会在后文讲解到
CodeForces 311B/CSU 1963 Cats Transport/Feed The Rabbit（斜率优化dp）
alpc_qleonardo
07-23 796
突然冒出来一道斜率优化的题，我都快忘了……赶紧翻原来写的博客……         题意也不是很好理解……有n座山，m只猫，每只猫都在固定一座山上，每座山线性排列，告诉你山之间的距离，猫i在某个时刻ti出现，然后总共有p个饲养员，他们从第一座山出发，单位时间走单位距离，而且都是一只往前走不停，每次路过可以把出现的猫领走，然后饲养员出发时间可以不同。问如何安排才能够使得所有猫的总等待时间最少，输出这个时间。不管你有没有题意了，没懂自己再去看题……
『DP·斜率优化·二分·CDQ分治』再谈神仙题「任务安排」①-④
Ronaldo7_ZYB的博客
05-10 339
前言 这个人很菜，在寒假的时候只写了任务安排1-2的题解；因此任务安排1-2是由当年的zyb写的。现在对原来的博客进行了补充，写了任务安排3-4.其中任务安排1是一个基础的DP，任务安排2是一个典型的斜率优化；3和4则是对斜率优化的一个延伸，需要利用特殊的方法来进行维护。 任务安排 Description There is a sequence of N jobs to be processed ...
#HDU2829#Lawrence(DP+ 斜率优化)
ElaineZXY的博客
05-23 998
Lawrence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 4109    Accepted Submission(s): 1873 Problem Description T. E. Lawrence was a
各种DP的总结
SemiWaker的博客
03-06 3002
各种DP的总结By SemiWaker
【数位dp】[Scoi2014] bzoj3598 方伯伯的商场之旅
每天心塞一点点
03-13 2512
题目点这里 和方伯伯的斗争终于结束了。。。我也是快要死了。。。。（请自动忽视那两道非人哉的题 = =） 之前听学长说这题很水的数位dp ：） 对 水到我都做不起了。  = =题解看了三遍啊！！！还是没研究明白这踏马是什么鬼！！ ……总有种我数位dp白学了的感觉（其实确实也是白学了）  思路……当然不是我想的………… 思路请见省选rank3大神：http
【二维树状数组优化dp】[Scoi2014] bzoj3594方伯伯的玉米田
每天心塞一点点
03-13 1995
终于方伯伯系列就只差逛商场了（两天的压轴题是人能做的？） 这道题记得当时我去省选的时候还把n^2的算法写出来了的……现在想想真不容易 = = 状态很好定义 dp[i][j] 表示已经处理完i根玉米 拔高了j次 剩下的玉米的最大值 我去省选的时候想的转移是这样的（其实我不确定） dp[i][j] = max(dp[i-1][j], dp[k][j-1]) (a[k] - a[i] =
【动态规划】【双调TSP与MTSP问题】 hdu2224 & hdu4281
每天心塞一点点
04-08 1898
一个月前写的了。。。。 首先双调TSP：hdu2224 给出平面上的N个点（1~N） 从 1 号点出发走到最右端的点 N 再返回 1  要求中途不能重复走某个点而且必须把所有点全走一遍 求最小的总路径 这个东西不是NPC 算导上有介绍O（N ^ 2）的算法。 dp[ i ][ j ]来表示从 1 到 i 以及从 1 到 j 的最小总路径 （i 转移： ①当 i ②
【数位dp？】CodeForces 288E Polo the Penguin and Lucky Numbers
每天心塞一点点
03-27 1699
题目点这里 那么好吧终于把这道题折腾完了 以及我并不觉得这是dp = = 题意：所有数位仅含4和7的数为lucky number 设[L, R]中第i大的lucky number为ai 求 a1*a2 + a2*a3 + a3*a4 +... an-1*an L和R的范围是10的10000次方 好吧 这数据范围就是数位dp的标准范围了 = =（不！其实一般的数位dp
【动态规划】[NOIP2014] codevs3729飞扬的小鸟
每天心塞一点点
09-12 1150
去年这题炸了，然后一年过去了，我又写炸了 人老了，写不动代码了…… 虽然复杂度是O(NM)没错，但是常数越写越大简直快要飞起来了 = = 其实感觉可以开个辅助数组会降不少常数…… 题目点这里 dp[i][j]表示到达(i, j)这个位置需要的最小点击次数 三个转移： 向上点击不到顶：step[i][j] = min(step[i][j], step[i - 1][j
动态规划优化：四边形不等式、斜率优化与单调队列
"圆弧与贝塞尔曲线互相转换" 的文章似乎与给定的标签和内容不符，因为提供的标签和内容主要涉及动态规划(DP)的优化策略，包括四边形不等式、斜率优化、单调队列优化以及状态压缩。以下是根据这些标签和部分内容生成...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值