HUST 1644 (思维)

本文介绍了一道来自第四届“恒生杯”程序设计大赛的题目,旨在寻找直方图中最大矩形面积。文章详细解析了解题思路,并提供了完整的C++代码实现。

1644 - 第四届“恒生杯”程序设计大赛决赛 B

时间限制:1秒 内存限制:256兆

70 次提交 18 次通过
题目描述

Given N non-negative integers representing the histogram's bar height where the width of each bar is 1, find the area of largest rectangle in the histogram.

 

Above is a histogram where width of each bar is 1, given height = [2,1,5,6,2,3].

 

 

The largest rectangle is shown in the shaded area, which has area = 10 unit.

 

For example,

Given height = [2,1,5,6,2,3],

 

return 10.

输入

 

The first line contains the number of test cases T. T test cases follow. Each case contains an integer N, followed by integers height[1],...,height[N] on the second line.

 

1 <= T <= 100 

1 <= N <= 100000 

 

1 <= height[i] <= 1000

输出

 

Output T lines, one for each test case, containing the largest rectangle area.


这题是在校内训练的时候给学弟在hust上找的,当时觉得人过的不少,难度应该一般。结果确实没人做出来了,实际上确实是一道不错的题,hust上的题基本上是从leetcode[84]抄过来的,但数据好像不强。

题解和代码如下:

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <stack>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <cmath>
#include <cstdlib>
#include <ctime>
using namespace std;
#define REP(i,k,n) for(int i=k;i<n;i++)
#define REPP(i,k,n) for(int i=k;i<=n;i++)
#define scan(d) scanf("%d",&d)
#define scann(n,m) scanf("%d%d",&n,&m)
#define mst(a,k)  memset(a,k,sizeof(a));
#define LL long long
#define eps 1e-8
using namespace std;
/*
题目要求最大面积,即要找到一对左右边界使得其中面积正好最大。
我们维护一个递增栈:从栈底到栈顶,元素值递增
现在从左到右把柱子扫一遍,如果栈为空或者heigt[i]>=height[top],则入栈
否则:将栈顶元素出栈。
对于这个元素来说,它可形成的面积的左边界就是栈中的它下面的那个元素(不包括它)(因为栈是递增的,所以它下面的那个元素比它小)
右边界即它当前迫使它出栈的这个柱子的前一个。因为如果某两个元素同时在栈中,先出栈的那个元素肯定是比后出栈的元素大的,对不对?所以在
某元素出栈之前,它前面还有元素出栈了(同时和这个元素处在栈中过),那前面出栈的这些元素肯定都是可以与它组成面积的(高度比它高)
这样,对于每一个出栈的元素,我们都可以很轻松的找到它的边界,也就算出了面积。只要在扫的过程中记录下最大值就可以了

但是这里还要注意一点:当把栈顶弹出后,栈为空了,它的左边界即为0,因为栈都被弹空了,说明左边已经没有比它小的数了。

*/
int height[100005];
int main()
{
    int t,n;
    scan(t);
    while(t--)
    {
        scan(n);
        REP(i,0,n) scan(height[i]);
        height[n]=0; //在最后添加一个高度为0的柱子,使得最后栈能“清空”
        int ans=0;
        stack<int>st;
        for(int i = 0; i <= n; i++)
        {
            while(!st.empty() && height[i] < height[st.top()])
            {
                int val = st.top();
                st.pop(); //出栈!!
                //ans = max(ans, height[val]*(i-1-(st.empty()?-1:st.top())));
                if(st.empty())  //如果为空,左边界为0,右边界i-1,长度为i-1+1。
                    ans = max(ans,height[val]*i);
                else    //不为空,左边界新的栈顶+1,右边界i-1,长度i-1-(st.top+1)+1
                    ans = max(ans,height[val]*(i-1-st.top()));
            }
            st.push(i);
        }
        printf("%d\n",ans);
    }
    return 0;
}



内容概要:本文介绍了一种基于蒙特卡洛模拟和拉格朗日优化方法的电动汽车充电站有序充电调度策略,重点针对分时电价机制下的分散式优化问题。通过Matlab代码实现,构建了考虑用户充电需求、电网负荷平衡及电价波动的数学模【电动汽车充电站有序充电调度的分散式优化】基于蒙特卡诺和拉格朗日的电动汽车优化调度(分时电价调度)(Matlab代码实现)型,采用拉格朗日乘子法处理约束条件,结合蒙特卡洛方法模拟大量电动汽车的随机充电行为,实现对充电功率和时间的优化分配,旨在降低用户充电成本、平抑电网峰谷差并提升充电站运营效率。该方法体现了智能优化算法在电力系统调度中的实际应用价值。; 适合人群:具备一定电力系统基础知识和Matlab编程能力的研究生、科研人员及从事新能源汽车、智能电网相关领域的工程技术人员。; 使用场景及目标:①研究电动汽车有序充电调度策略的设计与仿真;②学习蒙特卡洛模拟与拉格朗日优化在能源系统中的联合应用;③掌握基于分时电价的需求响应优化建模方法;④为微电网、充电站运营管理提供技术支持和决策参考。; 阅读建议:建议读者结合Matlab代码深入理解算法实现细节,重点关注目标函数构建、约束条件处理及优化求解过程,可尝试调整参数设置以观察不同场景下的调度效果,进一步拓展至多目标优化或多类型负荷协调调度的研究。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值