数据结构基础【单调栈】

单调栈

定义：

• 单调递增栈：出栈时为单调递增序列。
• 单调递增栈：出栈时为单调递增序列。

单调栈的经典运用：
POJ2559Largest Rectangle in a Histogram
题目大意：
如下图所示，在一条水平线上方有若千个矩形，求包含于这些矩形的并集内部的最大矩形的面积(在下图中，答案就是阴影部分的面积)，矩形个数≤105。

思路：
我们建立一个栈，用来保存若干个矩阵，这些矩阵的高度是单调递增的，然后我们从左到右一次扫描每个矩阵。
如果当前矩形比栈顶矩形高，直接进栈，否则不断取出栈顶，直至栈顶空或者栈顶矩形高度比当前矩形小.在出栈程中，我们累计被弹出的矩形的宽度之和，每弹出一个矩形，就用它的高度乘上累计的宽度去更新答案.整到出栈过程结束后，我们把一个高度当成当前矩形高度，宽度为累计值的新矩形入栈.
整个扫描结束后，我们把栈中剩余的矩形依次弹出，按照与上面相同的方法更新答案。为了简化程序实现，也可以增加一 个高度为0的矩形a[n+ 1]，以避免在扫描结束后栈中有剩余矩形。

#include<iostream>
#include<cstring>
#include<stack>
#include<cstdio>
using namespace std;
typedef long long ll;
const int N = 1e5 + 10;
ll  ans, a[N], s[N], w[N];
int n, p;
int main(){
	while(~scanf("%d", &n) && n){
		for(int i = 1; i <= n; i++){
			scanf("%lld", &a[i]);
		}
		ans = 0;
		a[n + 1] = 0; //避免在扫描完以后有剩余的矩形 
		p = 0; //栈顶指针 
		for(int i = 1; i <= n + 1; i++){
			if(a[i] > s[p]){
				s[++p] = a[i]; //找到比栈顶元素高的就进栈 
				w[p] = 1; 
			}
			else{//否则不断取出栈顶，直至栈顶空或者栈顶矩形高度比当前矩形小
				int width = 0;
				while(s[p] > a[i]){
					width += w[p];  //向左寻找可以组成的矩形的宽度 
					ans = max(ans, (ll)width * s[p]);
					p--;
				}
				s[++p] = a[i];
				w[p] = width + 1;
			}
		}
		printf("%lld\n", ans);
	}
	return 0;
}