单调栈和暴力拓展解决Largest Rectangle in a Histogram

题目需要求在一些柱体中取得最大的矩形。

先讲一下暴力(雾，不知道叫暴力还是递推好)的方法:通过比他大或等于的元素进行向左向右推广。
先看一下向左推广的代码：

for (int i = 1; i <= n; i++) {	
	while (l[i]>1&&s[i]<=s[l[i]-1]) {
		l[i] = l[l[i] - 1];
	}		
}

用l[i]表示第i个矩形的左边界，即有第l[i]-1个元素小于第i个元素，因为和该元素比较的元素是大于或等于它的，所以比较元素的左边界不一定是此元素的左边界。右边界同理可得，最后以某个矩阵的高作矩形时有长乘以高即可，他们的最大值即为答案。

此方法AC代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<functional>
#include<deque>
#include<cmath>

using namespace std;

typedef long long ll;

const int N = 1e6 + 10;

struct node {
	int val;
	int id;
}date[N];

int s[N];

int l[N], r[N];

deque<node> dq;

int main() {
	int n;
	while (scanf("%d", &n) && n) {
		ll ans = 0;
		for (int i = 1; i <= n; i++) {
			scanf("%d", &s[i]);
			l[i] = r[i] = i;//默认左右边界为它本身
		}
		for (int i = 1; i <= n; i++) {//获得左边界
			while (l[i]>1&&s[i]<=s[l[i]-1]) {
				l[i] = l[l[i] - 1];
			}
		}
		for (int i = n; i >= 1; i--) {//获得又边界
			while (r[i]<n&&s[i]<=s[r[i]+1]) {
				r[i] = r[r[i] + 1];
			}
		}
		for (int i = 1; i <= n; i++) {
			ans = max(ans,(ll)s[i]*(r[i]-l[i]+1));
		}
		printf("%lld\n", ans);
	}
}

然后是比较容易想到的单调栈做法，其实本质是一样的，也是确定左右边界，所以我们每次将某位置的元素入栈的时候，都要比较栈顶(栈顶元素是连续递增的)元素，知道找到小于插入元素的元素的位置，次位置为左/右边界，注意此方法中的左右边界和上面那个方法的左右边界不太一样，这里的左右边界位置的元素是小于插入元素的，上边的左右边界位置的元素是大于或等于插入元素的。

单调栈AC代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<stack>
using namespace std;

const int N = 1e6 + 10;

int a[N], l[N], r[N];
stack<int> s;

int main()
{
	int n;
	while (~scanf("%d", &n) && n) {
		while (!s.empty()) s.pop();//清空栈
		for (int i = 1; i <= n; i++) {
			scanf("%d", &a[i]);
		}
		a[n + 1] = a[0] = -1;
		s.push(0);
		for (int i = 1; i <= n; i++) {
			while (!s.empty() && a[s.top()] >= a[i]) s.pop();//找到左边界
			l[i] = s.top();//这里就是为什么上面要将0入栈的原因，0最小，所以永远在栈底
			s.push(i);
		}
		while (!s.empty()) s.pop();
		s.push(n + 1);
		for (int i = n; i >= 1; i--) {
			while (!s.empty() && a[s.top()] >= a[i]) s.pop();//找右边界
			r[i] = s.top();
			s.push(i);
		}
		long long ans = 0;
		for (int i = 1; i <= n; i++) {
			long long area = 0;
			ans = max(ans, (r[i] - l[i] - 1ll) * a[i]);//注意我们这里得到的左右边界是左右两侧第一个小于此元素的元素的位置
		}
		printf("%lld\n", ans);

	}
	return 0;
}