85. Maximal Rectangle

Given a 2D binary matrix filled with 0's and 1's, find the largest rectangle containing only 1's and return its area.

For example, given the following matrix:

1 0 1 0 01 0 1 1 1 1 1 1 1 1 1 0 0 1 0

Return 6.

1. 基于Largest Rectangle in Histogram

假设把矩阵沿着某一行分开，然后把分开的行作为底面，将自底面往上的矩阵看成一个直方图（histogram）。直方图的中每个项的高度就是从底面行开始往上1的数量。根据 Largest Rectangle in Histogram 就可以求出当前行作为矩阵下边缘的一个最大矩阵。接下来如果对每一行都做一次 Largest Rectangle in Histogram ，从其中选出最大的矩阵，那么它就是整个矩阵中面积最大的子矩阵。

如何计算某一行为底面时直方图的高度呢？如果重新计算，那么每次需要的计算数量就是当前行数乘以列数。然而会发现一些动态规划的踪迹，如果知道上一行直方图的高度，就只需要看新加进来的行（底面）上对应的列元素是不是0，如果是，则高度是0，否则则是上一行直方图的高度加1。利用历史信息，就可以在线行时间内完成对高度的更新。由于 Largest Rectangle in Histogram 的算法复杂度是O(n)。所以完成对一行为底边的矩阵求解复杂度是O(n+n)=O(n)。接下来对每一行都做一次，那么算法总时间复杂度是O(m*n)。

时间复杂度：O(mn)

空间复杂度：O(n)

class Solution {

public:

int maximalRectangle(vector<vector<char>>& matrix) {

if(matrix.size() == 0 || matrix[0].size() == 0) return 0;

int m = matrix.size(), n = matrix[0].size();

int res = 0;

vector<int> heights(n, 0); //记录每一列中从这一列起向上有几个1

for(int i = 0; i < m; i++){

for(int j = 0; j < n; j++){

heights[j] = matrix[i][j]== '0'? 0:heights[j]+1;

}

res = max(res, largestRectangleArea(heights));

}

return res;

}

private:

int largestRectangleArea(vector<int>& heights){

int max_size = 0;

int index = 0;

stack<int>s;

heights.push_back(0);

while(index<heights.size()){

if(s.size()==0 || heights[s.top()] < heights[index]){

s.push(index);

index++;

}else{

int top = s.top();

s.pop();

int size = 0;

if(s.size()==0){

size = heights[top]*index;

}else{

size = heights[top] * (index-s.top()-1);

}

max_size = max(max_size,size);

}

}

return max_size;

}

};