题目:

Given a 2D matrix matrix, find the sum of the elements inside the rectangle defined by its upper left corner (row1, col1) and lower right corner (row2, col2).


The above rectangle (with the red border) is defined by (row1, col1) = (2, 1) and (row2, col2) =(4, 3), which contains sum = 8.

Example:

Given matrix = [
[3, 0, 1, 4, 2],
[5, 6, 3, 2, 1],
[1, 2, 0, 1, 5],
[4, 1, 0, 1, 7],
[1, 0, 3, 0, 5]
] sumRegion(2, 1, 4, 3) -> 8
sumRegion(1, 1, 2, 2) -> 11
sumRegion(1, 2, 2, 4) -> 12

Note:

  1. You may assume that the matrix does not change.
  2. There are many calls to sumRegion function.
  3. You may assume that row1 ≤ row2 and col1 ≤ col2.

分析:

给定一个二维矩阵,求其子矩阵内所有元素的和。

如果每次调用sumRegion,遍历范围内所有元素的和,这种办法可行但是并不可取。我们先来看下面这个图。

上图是求解一个矩形面积的图画 展示,可以清楚的理解各个矩形块之间的关系,我们把矩形内的元素和比作矩形的面积。定义dp为二维数组,dp[i][j]是以(i, j)元素为右下角,可以求得一个子矩阵内所有元素的和。那么我们把dp[i][j]的值想象成上图中矩形的面积,大矩形的面积=去除右边一列的面积+去除下边一行的面积-重复的地方(左上角的面积)+当前小矩形块的面积。

此时当前的小矩形块就是所给定的矩阵(i, j)位置的元素。

去除右边一列的面积也就是以(i, j-1)元素为右下角的矩阵的元素和,实际上就是dp[i][j-1]。

去除下边一行的面积就是以(i-1, j)元素为右下角的矩阵的元素和,实际上就是dp[i-1][j]。

左上角的面积则是以(i-1, j-1)元素为右下角的矩阵的元素和,实际上就是dp[i-1][j-1]。

借此我们可以在O(mn)的时间内求的dp矩阵,接下来要求任意给定子矩阵的所有元素和,看下图。

同样还是将矩阵内的元素和比作矩形的面积,此时右下角为我们所要求的矩形面积=大矩形面积-去除右边区域的面积+去除下边区域的面积+重复的地方(左上角的面积)。

此时我们就可以利用上边所求的dp数组来快速求的结果。由于给了范围,也就是x1,y1,x2,y2,通过上面的图我们不难看出各个面积区域之间的关系。

实际上最后的结果也就是右下角的矩形面积。那么

大矩形面积也就是以x2, y2为右下角的矩阵的元素和,实际上就是dp[x2][y2]。

去除右边区域的面积就是以x2, y1-1为右下角的矩阵的元素和,实际上就是dp[x2][y1-1]。

去除下边区域的面积就是以x1-1, y2为右下角的矩阵的元素和,实际上就是dp[x1-1][y2]。

左上角的面积就是以x1-1, y1-1为右下角的矩阵的元素和,实际上就是dp[x1-1][y1-1]。

最后的结果就是dp[x2][y2] - dp[x2][y1-1] - dp[x1-1][y2] + dp[x1-1][y1-1]。

小技巧就是开辟的dp数组长宽可以比原数组多1,这样就不用了处理边界条件了。

程序:

C++

class NumMatrix {
public:
NumMatrix(vector<vector<int>>& matrix) {
if(matrix.empty())
return;
int m = matrix.size();
int n = matrix[].size(); dp = vector<vector<int>>(m+, vector<int>(n+, ));
for(int i = ; i <= m; ++i){
for(int j = ; j <= n; ++j){
dp[i][j] = dp[i-][j] + dp[i][j-] - dp[i-][j-] + matrix[i-][j-];
}
}
} int sumRegion(int row1, int col1, int row2, int col2) {
return dp[row2+][col2+] - dp[row2+][col1] - dp[row1][col2+] + dp[row1][col1];
}
private:
vector<vector<int>> dp;
};

Java

class NumMatrix {

    public NumMatrix(int[][] matrix) {
if(matrix.length == 0 || matrix == null)
return;
int m = matrix.length;
int n = matrix[0].length;
dp = new int[m+1][n+1];
for(int i = 1; i <= m; ++i){
for(int j = 1; j <= n; ++j){
dp[i][j] = dp[i-1][j] + dp[i][j-1] - dp[i-1][j-1] + matrix[i-1][j-1];
}
} } public int sumRegion(int row1, int col1, int row2, int col2) {
return dp[row2+1][col2+1] - dp[row2+1][col1] - dp[row1][col2+1] + dp[row1][col1];
}
private int[][] dp;
}

最新文章

  1. 【洛谷&#183;P2320】鬼谷子的钱袋
  2. Ext.NET-基础篇
  3. selenium 配合sikuli script操作高德地图
  4. Eclipse设置代码模版
  5. 学习使用monkey 测试
  6. 4.接口隔离原则(Interface Segregation Principle)
  7. WCF编程系列(一)初识WCF
  8. 段(SEGMENT)
  9. redis memcache
  10. UE4里的渲染线程
  11. add,update,list.jsp源码
  12. C#:几种数据库的大数据批量插入
  13. 二叉树(LeetCode) C++相关知识代码 系列1
  14. Linux环境下Hadoop集群搭建
  15. cookie和session的讲解
  16. ----regular expression in js----
  17. Django 笔记(四)模板标签 ~ 自定义过滤器
  18. ubuntu 下搭建redis和php的redis的拓展
  19. vue-cli 搭建的项目处理不同环境下请求不同域名的问题
  20. shell监控之列出1小时内cpu占用最多的10个进程

热门文章

  1. .NET Core 3.1之深入源码理解HealthCheck(二)
  2. 《C++Primer》第五版习题解答--第四章【学习笔记】
  3. DevOps is Hard、DevSecOps is Even Harder . --- Enterprise Holdings
  4. The command &#39;/bin/sh -c unzip -o php-7.2.23-src.zip&#39; returned a non-zero code: 1
  5. cogs 2109. [NOIP 2015] 运输计划 提高组Day2T3 树链剖分求LCA 二分答案 差分
  6. H5录音音频可视化-实时波形频谱绘制、频率直方图
  7. 6441. 【GDOI2020模拟01.17】小 ω 维护序列
  8. Python思维导图(二)—— 数据类型
  9. 2.Scanner的进阶使用
  10. Serilog高级玩法之用Serilog记录所选终结点附加属性