力扣304. 二维区域和检索 - 矩阵不可变（动态规划，二维前缀和）

力扣304. 二维区域和检索 - 矩阵不可变（动态规划，二维前缀和）

给定一个二维矩阵，计算其子矩形范围内元素的总和，该子矩阵的左上角为 (row1, col1) ，右下角为 (row2, col2) 。
上图子矩阵左上角 (row1, col1) = (2, 1) ，右下角(row2, col2) = (4, 3)，该子矩形内元素的总和为 8。

示例：

给定 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

提示：

你可以假设矩阵不可变。
会多次调用 sumRegion 方法。
你可以假设 row1 ≤ row2 且 col1 ≤ col2 。

动态规划，二维前缀和

状态方程

mymatrix[i + 1][j + 1] = mymatrix[i][j + 1] + mymatrix[i + 1][j] - mymatrix[i][j] + matrix[i][j];//状态方程

复杂度分析

时间复杂度：初始化 O(mn)，每次检索 O(1)，其中 m 和 n 分别是矩阵 matrix 的行数和列数。
初始化需要遍历矩阵 matrix 计算二维前缀和，时间复杂度是 O(mn)。
每次检索的时间复杂度是 O(1)。

空间复杂度：O(mn)，其中 m 和 n 分别是矩阵 matrix 的行数和列数。需要创建一个 m+1 行 n+1 列的二维前缀和数组 sums。

// 304NumMatrix.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。
//

#include <iostream>
#include<vector>
using namespace std;

class NumMatrix {
public:
    vector<vector<int>> mymatrix;
    NumMatrix(vector<vector<int>>& matrix) {
        
        if (matrix.size() >0)//鲁棒性，防止非空
        {
            int rowlen = matrix.size();
            int collen = matrix[0].size();
            vector<int> tempmatrix(max(rowlen, collen)+1, 0);
            mymatrix.resize(max(rowlen, collen)+1, tempmatrix);
            for (int i = 0; i < rowlen; i++)
            {
                for (int j = 0; j < collen; j++)
                {
                    mymatrix[i + 1][j + 1] = mymatrix[i][j + 1] + mymatrix[i + 1][j] - mymatrix[i][j] + matrix[i][j];//状态方程
                }
            }
        }
        for (int i = 0; i < mymatrix.size(); i++)
        {
            for (int j = 0; j < mymatrix[0].size(); j++)
            {
                cout << mymatrix[i][j] << " ";
            }
            cout << endl;
        }
    }

    int sumRegion(int row1, int col1, int row2, int col2) {
        //if (row1 == 0 && col1 != 0)return mymatrix[row2 + 1][col2 + 1] - mymatrix[row2 + 1][col1];
        //if (row1 != 0 && col1 == 0)return mymatrix[row2 + 1][col2 + 1] - mymatrix[row1][col2 + 1];
        //if (row1 == 0 && col1 == 0)return mymatrix[row2 + 1][col2 + 1];
        //如果mymatrix数组没有多一行的话，就要对row1和col1分情况了，1顶上边row1 == 0，2顶左边col1 == 0，3顶左上角row1 == 0 && col1 == 0，4不顶
        return mymatrix[row2 + 1][col2 + 1] - mymatrix[row1][col2 + 1] - mymatrix[row2 + 1][col1] + mymatrix[row1][col1];
    }
};

int main()
{
    vector<vector<int>> 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} };
    for (int i = 0; i < matrix.size(); i++)
    {
        for (int j = 0; j < matrix[0].size(); j++)
        {
            cout << matrix[i][j] << " ";
        }
        cout << endl;
    }
    NumMatrix* obj = new NumMatrix(matrix);
    cout << obj->sumRegion(0, 0, 0, 3) << endl;
    std::cout << "Hello World!\n";
}