题目:
给一个直方图,求直方图中的最大矩形的面积。例如,下面这个图片中直方图的高度从左到右分别是2, 1, 4, 5, 1, 3, 3, 他们的宽都是1,其中最大的矩形是阴影部分。
Input
输入包含多组数据。每组数据用一个整数n来表示直方图中小矩形的个数,你可以假定1 <= n <= 100000. 然后接下来n个整数h1, …, hn, 满足 0 <= hi <= 1000000000. 这些数字表示直方图中从左到右每个小矩形的高度,每个小矩形的宽度为1。 测试数据以0结尾。
Output
对于每组测试数据输出一行一个整数表示答案。
Sample Input
7 2 1 4 5 1 3 3
4 1000 1000 1000 1000
0
Sample Output
8
4000
解题思路:
对于每一个小矩形,在这个高度要想面积尽量大肯定要向左和向右尽量延伸,每当向左和向右遇到比它小的矩形的时候停下来。所以我们维护一个单调递减栈,也就是入栈元素大于栈顶元素才入栈,这样把小矩形从左到右入栈的时候,每当遇见入栈矩形比栈顶小的时候出栈,直到栈顶小于入栈矩形,那么每一次出栈的矩形到此时正在准备入栈的矩形之间的距离就是这个出栈的矩形向右延申的最大长度。向左延申的最大长度就是把小矩形从右到左入栈。
所以用两次维护单调栈,记录每个矩形向右和向左延申的最大距离再找出最大面积。
代码:
#include<iostream>
using namespace std;
long long getmax(long long a,long long b)
{
if(a>=b)
return a;
else
return b;
}
long long h[100005];
int stackl[100005];
int stackr[100005];
int main()
{
int n;
while(cin>>n)
{
if(n==0)
break;
else
{
for(int i=1;i<=n;i++)
cin>>h[i];
}
h[n+1]=0;
h[0]=0;
int top=0;
int r[n+5];
for(int i=1;i<=n+1;i++)
{
while(1)
{
if(top==0||h[i]>=h[stackl[top]])
{
top++;
stackl[top]=i;
break;
}
else
{
r[stackl[top]]=i;
//maxl=getmax(maxl,(i-stackl[top])*h[stackl[top]]);
top--;
}
}
}
top=0;
int *l=new int[n+5];
for(int i=n;i>=0;i--)
{
while(1)
{
if(top==0||h[i]>=h[stackr[top]])
{
top++;
stackr[top]=i;
break;
}
else
{
l[stackr[top]]=i;
//maxr=getmax(maxr,(stackr[top]-i)*h[stackr[top]]);
top--;
}
}
}
long long max=0;
for(int i=1;i<=n;i++)
{
max=getmax(max,(r[i]-l[i]-1)*h[i]);
}
cout<<max<<endl;
}
return 0;
}