前缀和矩阵

前缀和矩阵(Prefix Sum Matrix)是一种预处理技术,用于快速计算二维矩阵中任意子矩阵的元素和。其核心思想是通过提前计算并存储每个位置左上角所有元素的和,将子矩阵和的查询时间从暴力计算的 (O(mn)) 优化到 (O(1))。以下是构建前缀和矩阵的详细步骤和示例:


文章目录

      • [1. 定义原矩阵与前缀和矩阵](#1. 定义原矩阵与前缀和矩阵)
      • [2. 构建前缀和矩阵的步骤](#2. 构建前缀和矩阵的步骤)
        • [步骤 1:初始化前缀和矩阵](#步骤 1:初始化前缀和矩阵)
        • [步骤 2:递推公式](#步骤 2:递推公式)
      • [3. 构建示例](#3. 构建示例)
      • [4. 查询子矩阵和](#4. 查询子矩阵和)
      • [5. 代码实现(C++)](#5. 代码实现(C++))
      • 关键点总结

1. 定义原矩阵与前缀和矩阵

  • 原矩阵 :一个 m X n 的二维数组 matrix,元素为整数或浮点数。
  • 前缀和矩阵 :一个与 matrix 同尺寸的二维数组 prefix,其中 prefix[i][j] 表示从左上角 (0,0)(i,j) 形成的子矩阵所有元素之和。

2. 构建前缀和矩阵的步骤

步骤 1:初始化前缀和矩阵
  • 创建一个与原矩阵大小相同的二维数组 prefix
  • 通常将 prefix 的索引从 (0,0) 开始(与 matrix 对齐)。
步骤 2:递推公式
  • 边界条件

    • i=0j=0 时:

      markdown 复制代码
      							prefix[0][0] = matrix[0][0]
    • i=0 时(第一行):

      markdown 复制代码
      					prefix[0][j] = prefix[0][j-1] + matrix[0][j]
    • j=0 时(第一列):

      markdown 复制代码
      					prefix[i][0] = prefix[i-1][0] + matrix[i][0]
  • 一般情况i>0j>0):

    markdown 复制代码
    				prefix[i][j] = matrix[i][j] + prefix[i-1][j] + prefix[i][j-1] - prefix[i-1][j-1]

    解释

    当前元素的值 matrix[i][j],加上上方子矩阵的和 prefix[i-1][j],加上左方子矩阵的和 prefix[i][j-1],再减去重复计算的左上角子矩阵 prefix[i-1][j-1]


3. 构建示例

假设原矩阵 matrix 如下:

markdown 复制代码
                                           [1  2  3]
                                           [4  5  6]
                                           [7  8  9]

逐行构建前缀和矩阵 prefix

  1. 初始化 prefix[0][0]

    markdown 复制代码
    							prefix[0][0] = matrix[0][0] = 1
  2. 第一行(i=0

    markdown 复制代码
    					prefix[0][1] &= prefix[0][0] + matrix[0][1] = 1 + 2 = 3 
                        prefix[0][2] &= prefix[0][1] + matrix[0][2] = 3 + 3 = 6 
  3. 第一列(j=0

    markdown 复制代码
    					prefix[1][0] &= prefix[0][0] + matrix[1][0] = 1 + 4 = 5 
                        prefix[2][0] &= prefix[1][0] + matrix[2][0] = 5 + 7 = 12 
  4. 一般位置(i>0j>0

    • prefix[1][1]

      markdown 复制代码
      								5 + 5 + 3 - 1 = 12
    • prefix[1][2]

      markdown 复制代码
      								6 + 12 + 6 - 3 = 21
    • prefix[2][1]

      markdown 复制代码
      								8 + 12 + 12 - 5 = 27
    • prefix[2][2]

      markdown 复制代码
      								9 + 27 + 21 - 12 = 45

最终前缀和矩阵

markdown 复制代码
                                            [1  3   6 ]
                                            [5  12  21]
                                            [12 27  45]

4. 查询子矩阵和

构建前缀和矩阵后,计算子矩阵 (x1, y1)(x2, y2) 的和公式为:

markdown 复制代码
	Sum = prefix[x2][y2] - prefix[x1-1][y2] - prefix[x2][y1-1] + prefix[x1-1][y1-1]

示例

计算子矩阵 (1,1)(2,2)(即原矩阵中的 5,6,8,9)的和:

markdown 复制代码
	Sum = 45 - 6 - 12 + 1 = 28

5. 代码实现(C++)

cpp 复制代码
#include <vector>
using namespace std;

vector<vector<int>> buildPrefixSum(vector<vector<int>>& matrix) {
    int m = matrix.size(), n = matrix[0].size();
    vector<vector<int>> prefix(m, vector<int>(n, 0));

    // 初始化第一行和第一列
    prefix[0][0] = matrix[0][0];
    for (int j = 1; j < n; j++) 
        prefix[0][j] = prefix[0][j-1] + matrix[0][j];
    for (int i = 1; i < m; i++) 
        prefix[i][0] = prefix[i-1][0] + matrix[i][0];

    // 填充其他位置
    for (int i = 1; i < m; i++) {
        for (int j = 1; j < n; j++) {
            prefix[i][j] = matrix[i][j] 
                          + prefix[i-1][j] 
                          + prefix[i][j-1] 
                          - prefix[i-1][j-1];
        }
    }
    return prefix;
}

关键点总结

  • 时间复杂度:构建前缀和矩阵需 (O(mn)),查询子矩阵和仅需 (O(1))。
  • 适用场景:频繁查询子矩阵和的场景(如动态规划、图像处理)。
  • 边界处理 :注意索引从 0 开始,避免越界访问。
相关推荐
努力学习的小廉21 小时前
深入了解linux系统—— 线程同步
linux·服务器·数据库·算法
数据爬坡ing21 小时前
从挑西瓜到树回归:用生活智慧理解机器学习算法
数据结构·深度学习·算法·决策树·机器学习
luoganttcc21 小时前
小鹏汽车 vla 算法最新进展和模型结构细节
人工智能·算法·汽车
骄傲的心别枯萎21 小时前
RV1126 NO.16:通过多线程同时获取H264和H265码流
linux·c++·音视频·rv1126
落羽的落羽21 小时前
【C++】特别的程序错误处理方式——异常机制
开发语言·c++
空山新雨(大队长)21 小时前
C 语言第一课:hello word c
c++·c·exe
春蕾夏荷_72829772521 小时前
c++ 第三方库与个人封装库
c++·三方库
牵牛老人1 天前
Qt C++ 复杂界面处理:巧用覆盖层突破复杂界面处理难题之一
数据库·c++·qt
wallflower20201 天前
滑动窗口算法在前端开发中的探索与应用
前端·算法
林木辛1 天前
LeetCode热题 42.接雨水
算法·leetcode