HDU 1506 不会则直方图最大矩形

最新推荐文章于 2018-08-11 20:03:32 发布

原创最新推荐文章于 2018-08-11 20:03:32 发布 · 356 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#dp #HDU

ACM_HDU 专栏收录该内容

64 篇文章

订阅专栏

本文介绍了一种计算由直方图构成的最大矩形面积的高效算法。通过预处理左右边界来减少重复计算，使用__int64类型确保大数值运算准确无误。此算法适用于解决诸如最大矩形面积等计算几何问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;//#include<algorithm>using namespace std;为了调用库函数max

__int64 dp[100005]; //因为高已经很大了乘上底会更大，所以用__int64
__int64 h[100005],l[100005],r[100005];

int main()
{
    __int64 n;
    while(scanf("%I64d",&n))//注意64位的输出格式
    {
        if(n==0)
            return 0;
        memset(dp,0,sizeof(dp));//养成习惯
        memset(h,0,sizeof(h));
        __int64 i;//注意i为了统一也设置成__int64
        for(i=1;i<=n;i++)
            scanf("%I64d",&h[i]);
        for(i=1;i<=n;i++)//初始化两边可以延展的范围
        {
            r[i]=i;
            l[i]=i;
        }
        for(i=1;i<=n;i++)//找左极限必须从1开始，保证前面每一个l[i]需要调用的l[i-1,...,1]都是已经设置好的
        {
            while(h[i]<=h[l[i]-1]&&l[i]>1)//只要不是第一个数，就和前面的数继续比较小于等于就继续向左延伸
                l[i]=l[l[i]-1];
        }
        for(i=n;i>=1;i--)//找右极限必须从n开始,保证r[i]右边的需要调用的r[i+1,...,n]都是已经设置好的
        {
            while(h[i]<=h[r[i]+1]&&r[i]<n)//只要不是最后一个数，就和前面的数继续比较小于等于就继续向右延伸
                r[i]=r[r[i]+1];
        }
        for(i=1;i<=n;i++)
        {
            dp[i]=h[i]*(r[i]-l[i]+1);//底需要加1
        }
        __int64 max=0;
        for(i=1;i<=n;i++)
            if(max<dp[i])
                max=dp[i];//求出组大面积
        printf("%I64d\n",max);
    }
    return 0;
}