Poj 2559 最大矩形面积

题目：

https://cn.vjudge.net/problem/POJ-2559

大意：给定一些长为1，高为x的连续矩形，求最大连续矩形面积

算法：单调栈；
题解：
维护一个高度单调递增的栈，栈中存结构体x，y，x表示高度，y表示以x为高度，向左扩展的最大宽度；
对于某一高度为h的矩形
1.当栈顶元素s.x小于h时，将｛h，1｝直接入栈；
2.当栈顶元素s.x大于h时，将栈顶元素弹出，并更新答案，重复步骤，直至栈顶元素小于h；
此时维护值k表示已弹出元素y的和，更新答案时，将s.x*(k+s.y)；
k维护的是从s到右的大于s.x的连续宽度，所以以s.x为高的矩形向右扩展，最多能扩展k；
y为定义，记录答案；
3.将｛h，k+1｝入栈；
k是从h向左大于h的连续宽度，是以h为高，向左扩展的最大宽度；再加上h的宽度1；
4.全部入栈完毕后，执行出栈操作；
仍然如2维护k值，更新答案；

样例：7  2 1 4 5 1 3 3     ans=8;
入栈  2  S <2,1> 
        1   pop<2,1> k=0,ans=(2,0)=2;
            S <1,2>
        4 S<1,2> <4,1>
        5 S<1,2> <4,1> <5,1>
        1 pop <5,1> k=0; ans=(2,5)=5;k=1; 
          pop <4,1> k=1; ans=(5,8)=8;k=2;
          pop <1,2> k=2; ans=(4,8)=8;k=4;
          push<1,5>
        3 S <1,5> <3,1>
        3 pop<3,1> k=0;ans=(3,8)=8;k=1;
          push<3,2> 
push S<1,5> <3,2>
        pop<3,2> k=0,ans=(6,8);k=2;
        pop<1,5> k=2,ans=(7,8);k=7; 

代码：

#include<iostream>
#include<cstdio>
#include<stack>
#include<algorithm>
using namespace std;
const int N=100000;
int n,h[N];
long long k,ans,now;
struct node{
    long long x,y;
};
stack<node> S;
void init(){
    now=0;
    k=0;
    while(!S.empty()) S.pop();
    ans=0;
}
int main(){
//  freopen("data.txt","r",stdin);
//freopen("my.out","w",stdout);
    while(scanf("%d",&n)!=EOF){

        if(n==0) return 0;
        init();
        for(int i=1;i<=n;i++) scanf("%d",&h[i]);
        S.push((node){h[1],1});
        for(int i=2;i<=n;i++){
            k=0;//cout<<i<<endl;
            while(1){
                node u=S.top();
                if(u.x>=h[i]){
                    now=u.x*(k+u.y);
                    ans=max(ans,now);
                    k+=u.y;
                    S.pop();
                }
                if(S.empty()){
                    S.push((node){h[i],k+1});
                    break;
                }
                if(u.x<h[i]){
                    S.push((node){h[i],k+1});
                    break;
                }
            }
        }
        k=0;
        while(!S.empty()){
            node u=S.top();
            now=u.x*(k+u.y);
            ans=max(ans,now);
            k+=u.y;
            S.pop();
        }
        printf("%lld\n",ans);
    }
}