直方图最大面积矩形&01矩阵中最大的全0(1)矩阵

求直方图中最大矩形的面积,用这个思路可以进一步解决求01矩阵中最大的全0或全1矩阵

直方图中最大矩形的面积

问题描述

给定一个一维数组,代表直方图中每个区间的高度,求出直方图中最大矩形的面积,图中区间的宽度为1。

图中所示一维数组为[2,1,5,6,2,3],那么所求的矩形应该是5和6区间构成的矩形,面积为5*2=10。

求解算法

这个问题的核心思想是依次计算每个区间能够构成的最宽矩形面积,再从中选出面积最大的
首先考虑一下有序序列,这样的直方图可以用统一的算法来求解每个区间的矩形面积:从后往前遍历,假设最后一个数值为x,它的索引为last,当前遍历到的数值索引为cur,那么当前区间构成的矩形面积就为x*(last-cur+1),再设置一个max变量,记录遍历过程中出现的最大值就行了。
但是问题在于这些数字都是无序的,不好用这种统一的方法来计算矩形面积,所以引入了栈结构,目的就是从序列中挑出有序的部分来计算,确切地说是保持栈内的元素都是升序排列,而且每个元素都入栈一次。

具体做法是:
1、将序列中的元素索引依次入栈,若当前索引对应的数值大于栈顶元素对应的数值时时,将索引顺利入栈;要是小于栈顶元素对应的数值,说明这是逆序,执行第2步;
2、把栈顶元素输出,得到这个索引对应的数值,计算这个数值代表区间的最大矩形面积(此时按照有序时的公式,last就是刚才要入栈但是发现逆序没入成的那个索引值-1,cur则是现在pop出来的索引值),计算完更新max变量,再返回第1步;
3、当所有元素都入栈之后,若栈内还有元素,就依次pop出来计算面积,实质上就是处理剩下的有序序列了。

代码实现

int getMaxArea(vector<int> array){
  int len=array.size();
  if (len==0) return 0;
  stack<int> s;
  int max=array[0];
  s.push(0);
  for (int i=1;i<len;){
    int cur=array[s.top()];
    if (array[i]>=cur){
      s.push(i);
      i++;
    }else{
      int top=array[s.top()];
      int curMax=top*(i-s.top());
      s.pop();
      if (curMax>max) max=curMax;
    }
  }
  if (!s.empty()){
    int last=s.top();
    while (!s.empty()){
      int top=array[s.top()];
      int curMax=top*(last-s.top()+1);
      s.pop();
      if (curMax>max) max=curMax;
    }
  }
  return max;
}

01矩阵中最大的全0或全1矩阵

问题描述

给定一个M*N的二维数组,其元素全部由0和1构成,找出最大的全1矩阵,输出这个矩阵的元素个数。

上图中最大全1矩阵就是红色标出的部分,因此输出元素个数6。

求解算法

有了直方图最大矩形面积的基础,这道题就可以拆解为同样的问题来解决:从第一行开始,计算以每一行为底的矩形构成的直方图中最大矩形面积,再在其中找到最大值,就是要求的矩形
解释一下以每一行为底的矩形构成的直方图,第一行的话,就是这个一维的数组[1,0,1,0,0],那么直方图显而易见就是这些数值,将这个数组送入上面那个函数,求得最大面积是1;遍历到第二行,此时以第二行为底的矩形是2*5的格式,也就是矩阵的上两行,它可以代表一个直方图的分布,那么这个2*5矩阵的最大全1矩形,就相当于所代表的直方图的最大矩形面积。
因此,需要将这个2*5的矩阵转换成一维数组的形式,也就是计算直方图每个区间有多少个1,而且是连续的1(因为不连续的话就构不成一个子矩阵了),所以对第二行的每个a[i][j]做以下处理:
1、若a[i][j]=0,那么它的值保持不变还是0,说明当前没有连续的1;
2、若a[i][j]=1,令a[i][j]=a[i-1][j]+1,说明这个区间又多了一个连续的1。
最后第二行得到的数组是[2,0,2,1,1],将其送入函数,得到最大矩形面积是3;同样的第三行数组是[3,1,3,2,2],得到最大矩形面积是6;第四行是[4,0,0,3,0],最大面积为4。遍历结束,得到最终的最大全1矩阵的面积为6。

代码实现

int getMaxRectangle(vector<vector<int>> array){
  int m=array.size();
  int n=array[0].size();
  int max=0;
  for (int i=0;i<m;i++){
    for (int j=0;j<n;j++){
      if (i!=0){
        array[i][j]=array[i][j]==0?0:array[i-1][j]+1;
      }
    }
    int curMax=getMaxArea(array[i]);
    if (curMax>max) max=curMax;
  }
  return max;
}

参考博客

算法题:直方图和0-1矩阵中最大矩形

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,658评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,482评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,213评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,395评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,487评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,523评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,525评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,300评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,753评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,048评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,223评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,905评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,541评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,168评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,417评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,094评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,088评论 2 352

推荐阅读更多精彩内容