【ATT】Maximal Rectangle

本质上是动态规划。

解法一:

 //O(m*m*n)   
    int maximalRectangle(vector<vector<char> > &matrix) {
        if(matrix.empty()||matrix[0].empty())
            return 0;
        int m = matrix.size();
        int n = matrix[0].size();
        
        vector<vector<int>> cnt(m,vector<int>(n,0)); //cnt[i][j]表示matrix[i][j]往左边延伸‘1’的最大个数。
        
        int i,j;
        
        for(i=0;i<m;i++)
            cnt[i][0] = (matrix[i][0]=='1'?1:0);
        for(i=0;i<m;i++)
        {
            for(j=1;j<n;j++)
            {
                if(matrix[i][j]=='1')
                    cnt[i][j] = cnt[i][j-1]+1;
            }
        }
        
        int max_area = 0;
        
        for(i=0;i<m;i++)  //计算以(i,j)作为矩形右下角的最大矩形面积
        {
            for(j=0;j<n;j++)
            {
                int curRow = i;
                int minCnt = INT_MAX;
                while(curRow>=0)
                {
                    if(cnt[curRow][j]==0)
                        break;
                    minCnt = min(minCnt,cnt[curRow][j]);
                    max_area = max(max_area,minCnt*(i-curRow+1));
                    --curRow;
                }
            }
        }
        
        return max_area;
    }

  

解法二:

//Two solutions; 
//O(m*n)
    int maxRecByRow(vector<int>& height)
    {
        int max_area = 0;
        stack<int> st;
        int i=0;
        
        while(i<height.size())
        {
            if(st.empty()||height[st.top()]<=height[i])
                st.push(i++);
            else
            {
                int tp = st.top();
                st.pop();
                int cur_area = height[tp]*(st.empty()?i:i-st.top()-1);
                max_area = max(max_area,cur_area);
            }
        }
        return max_area;
        
    }
    int maximalRectangle(vector<vector<char> > &matrix) {
        // Start typing your C/C++ solution below
        // DO NOT write int main() function
        if(matrix.empty()||matrix[0].empty())
            return 0;
        int m = matrix.size();
        int n = matrix[0].size();
        vector<vector<int>> heights(m,vector<int>(n+1,0)); //这里多一列,填充0,作为sentinal
//heights[i][j]表示第i行作为baseline,第j个矩形方阵的高度。 int i,j; for(i=0;i<m;i++) { for(j=0;j<n;j++) { if(matrix[i][j]=='1') heights[i][j] = (i==0?1:heights[i-1][j]+1); } } int max_area = 0; for(i=0;i<m;i++) max_area = max(max_area,maxRecByRow(heights[i])); return max_area; }

  

 

分析来自于:http://www.cnblogs.com/lichen782/p/leetcode_maximal_rectangle.html

LeetCode 笔记系列 18 Maximal Rectangle [学以致用]

题目: Given a 2D binary matrix filled with 0's and 1's, find the largest rectangle containing all ones and return its area.

leetcode的题目真是越来越经典了。比如这个题目,就是一道男默女泪的题。

一般人拿到这个题目,除非做过类似的,很难一眼找出一个方法来,更别说找一个比较优化的方法了。

首先一个难点就是,你怎么判断某个区域就是一个矩形呢?

其次,以何种方式来遍历这个2D的matrix呢?

一般来说,对这种“棋盘式”的题目,像什么Queen啦,象棋啦,数独啦,如果没有比较明显的遍历方式,可以采用一行一行地遍历。(好像废话哦。。。)

然后,当遍历到(i, j)的时候,该做什么样的事情呢?想想,嗯,那我可不可以简单看看,以(i,j)为矩形左上角,能不能形成一个矩形,能不能形成多个矩形?那形成的矩形中,我们能不能找一个最大的呢?(有同学问,为毛你要以这个点为左上角,不为左下角,或者其他脚哩?因为我们打算从左到右,从上到下一行一行遍历嘛,这样就不会漏掉,说不定还能做一些优化)

首先,如果(i, j)是0,那肯定没法是矩形了。

如果是1,那么我们怎么找以它为左上角的矩形呢?呼唤画面感!

。。。你TM在逗我?==b

图中圈圈表示左上角的1,那么矩形的可能性是。。。太多啦,怎么数嘛!

我们可以试探地从左上角的1所在的列开始,往下数数,然后呢,比如在第一行,例如是蓝色的那个矩形,我们看看在列上,它延伸了多远,这个面积是可以算出来的。

然后继续,第二行,例如是那个红色的矩形,再看它延伸到多远,哦,我们知道,比第一行近一些,我们也可以用当前离第一行的行数,乘以延伸的距离,得到当前行表示的矩形面积。

但是到了第一个虚线的地方,它远远超过了上面的其他所有行延伸的距离了,注意它的上方都是空心的哦,所以,我们遇到这种情况,计算当前行和左上角1围成的面积的时候,只能取所有前面最小的延伸距离乘以当前离第一行的行数。其实,这对所有情况都是这样的,是吧?于是,我们不是就有方法遍历这些所有的矩形了嘛。

代码如下:

复制代码
 1     /**
 2      * 以给出的坐标作为左上角,计算其中的最大矩形面积
 3      * @param matrix
 4      * @param row 给出坐标的行
 5      * @param col 给出坐标的列
 6      * @return 返回最大矩形的面积
 7      */
 8     private int maxRectangle(char[][] matrix, int row, int col) {
 9         int minWidth = Integer.MAX_VALUE;
10         int maxArea = 0;
11         for (int i = row; i < matrix.length && matrix[i][col] == '1'; i++) {
12             int width = 0;
13             while (col + width < matrix[row].length
14                     && matrix[i][col + width] == '1') {
15                 width++;
16             }
17             if (width < minWidth) {// 如果当前宽度小于了以前的最小宽度,更新它,为下面的矩形计算做准备
18                 minWidth = width;
19             }
20             int area = minWidth * (i - row + 1);
21             if (area > maxArea)
22                 maxArea = area;
23         }
24         return maxArea;
25     }
复制代码

这样,我们再对每个点都调用一下上面的这个方法,不是就能求出最大面积了么。

解法一:

 View Code

这个需要O(n3),所以没有通过大集合的测试。

leetcode的讨论组给出了一个比较难理解的方法,这里就不采用了。

说说第三个方法。前一个笔记,我们讨论了柱状图的最大矩形面积,那可以O(n)的,学以致用呀!btw,leetcode的这两题也是挨一块儿的,用心良苦。。。。

如果我们把每一行看成x坐标,那高度就是从那一行开始往上数的1的个数。带入我们的maxAreaInHist方法,在O(n2)时间内就可以求出每一行形成的“柱状图”的最大矩形面积了。它们之中最大的,就是我们要的答案。

代码如下:

复制代码
 1 public int maximalRectangle2(char[][] matrix) {
 2         int m = matrix.length;
 3         int n = m == 0 ? 0 : matrix[0].length;
 4         int[][] height = new int[m][n + 1];
 5         //actually we know that height can just be a int[n+1], 
 6         //however, in that case, we have to write the 2 parts together in row traverse,
 7         //which, leetcode just doesn't make you pass big set
 8         //所以啊,leetcode是喜欢分开写循环的,即使时间复杂度一样,即使可以节约空间
 9         int maxArea = 0;
10         for(int i = 0; i < m; i++){
11             for(int j = 0; j < n; j++) {
12                 if(matrix[i][j] == '0'){
13                     height[i][j] = 0;
14                 }else {
15                     height[i][j] = i == 0 ? 1 : height[i - 1][j] + 1;
16                 }
17             }
18         }
19         for(int i = 0; i < m; i++){
20             int area = maxAreaInHist(height[i]);
21             if(area > maxArea){
22                 maxArea = area;
23             }
24         }
25         return maxArea;
26      }
27      
28      private int maxAreaInHist(int[] height){
29          Stack<Integer> stack = new Stack<Integer>();
30          int i = 0;
31          int maxArea = 0;
32          while(i < height.length){
33              if(stack.isEmpty() || height[stack.peek()] <= height[i]){
34                  stack.push(i++);
35              }else {
36                  int t = stack.pop();
37                  maxArea = Math.max(maxArea, height[t] * (stack.isEmpty() ? i : i - stack.peek() - 1));
38              }
39          }
40          return maxArea;
41      }
复制代码

这里有一个和leetcode相关的细节。就是本来在计算height数组的时候,我们没有必要分配成代码中的那个样子,一维就可以了,然后在遍历每一行的时候计算当前行的height数组,然后再计算maxArea。这种情况下还是过不了大集合,所以不得不为每一行都保存一个height,先期计算该二维数组。

总结:

1. 学到的新知识要用;

2. 画面感和逻辑分析都很重要,不可偏非。

posted @ 2013-10-25 20:51  summer_zhou  阅读(183)  评论(0编辑  收藏  举报