POJ 2559 最大矩形面积 Histogram

#include <bits/stdc++.h>
#define Maxn 100010
using namespace std;
//sta[]是单调栈，栈中存储的元素高度递增，p为栈的栈顶位置
int h[Maxn],w[Maxn],n,sta[Maxn],p;//h:矩形的hight,w:矩形的width
long long ans=0;
int main()
{
    freopen("data1.in","r",stdin);
    while(scanf("%d",&n)!=EOF){
        memset(w,0,sizeof(w));
        for(int i=1; i<=n; ++i){scanf("%d",&h[i]);} //读入n个矩形的高度，每个矩形的默认宽度：1
        sta[0]=h[n+1]=p=0;  //边界的处理：sta[0]=0为了方便第1个矩形入栈，h[n+1]=0为了将栈中所有矩形出栈
        for(int i=1; i<=n+1; ++i){
            if(sta[p]<h[i]){ //当前矩形比栈顶矩形高度高的话入栈，不要写=，h[0]=0,而h[n+1]也是0,是出栈用的
                sta[++p]=h[i];
                w[p]=1;
            }else{
                int width=0;
                while(sta[p]>h[i]){//当前矩形比栈顶矩形高度低则出栈，不要写=，h[n+1]是0是出栈用的
                    width+=w[p];
                    ans=max(ans,(long long)width*sta[p]);
                    p--;
                }
                sta[++p]=h[i],w[p]=width+1; //更新当前矩形宽带(与前面矩形构成的最大宽带，木桶的短板)，同时将其入栈
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}