矩阵区域和(Java)

给一个 m * n 的矩阵 arr 和一个整数 K ，返回一个矩阵 result ，其中每个 result[i][j] 是所有满足下述条件的元素 arr[r][c] 的和：
    i - K <= r <= i + K, j - K <= c <= j + K
    (r, c) 在矩阵内。
 

示例 1：
输入：arr = [[1,2,3],[4,5,6],[7,8,9]], K = 1
输出：[[12,21,16],[27,45,33],[24,39,28]]

示例 2：
输入：arr = [[1,2,3],[4,5,6],[7,8,9]], K = 2
输出：[[45,45,45],[45,45,45],[45,45,45]]

提示：
    m == arr.length
    n == arr[i].length
    1 <= m, n, K <= 100
    1 <= arr[i][j] <= 100

用数组 P 表示数组 arr 的二维前缀和，P 的维数为 (m + 1) * (n + 1)，其中 P[i][j] 表示数组 arr 中以 (0, 0) 为左上角，(i - 1, j - 1) 为右下角的矩形子数组的元素之和。
需要对数组 arr 中的每个位置，计算以 (i - K, j - K) 为左上角，(i + K, j + K) 为右下角的矩形子数组的元素之和，可以在前缀和数组的帮助下，通过：
sum = P[i + K + 1][j + K + 1] - P[i - K][j + K + 1] - P[i + K + 1][j - K] + P[i - K][j - K]
得到元素之和。注意 i + K + 1、j + K - 1、i - K 和 j - K 这些下标有可能不在矩阵内，因此对于所有的横坐标，需要将其规范在 [0, m] 的区间内；对于所有的纵坐标，我们需要将其规范在 [0, n] 的区间内。具体地：
    i + K + 1 和 j + K - 1 分别可能超过 m 和 n，因此我们需要对这两个坐标与 m 和 n 取较小值，忽略不在矩阵内的部分；
    i - K 和 j - K 可能小于 0，因此我们需要对这两个坐标与 0 取较大值，忽略不在矩阵内的部分。
一般的，对所有的横坐标与 m 取较小值，纵坐标与 n 取较小值，再将所有坐标与 0 取较大值，就可以将这些坐标规范在前缀和数组 P 的范围内。

package com.loo;

import java.util.Arrays;

public class MatrixAreaSum {

    public static void main(String[] args) {
        int[][] arr = {{1,2,3},{4,5,6},{7,8,9}};
        int k = 2;
        int[][] result = getMatrixAreaSum(arr , k);
        if (result!=null && result.length!=0) {
            for (int i=0;i<result.length;i++) {
                System.out.println(Arrays.toString(result[i]));
            }
        }
    }
    
    public static int[][] getMatrixAreaSum(int[][] arr , int k) {
        if (arr == null || arr.length == 0 || arr[0].length == 0) {
            return new int[0][];
        }
        int m = arr.length;
        int n = arr[0].length;
        int[][] result = new int[m][n];
        int[][] dp = new int[m+1][n+1];
        for (int i=1;i<=m;i++) {
            for (int j=1;j<=n;j++) {
                dp[i][j] = arr[i-1][j-1] + dp[i-1][j] + dp[i][j-1] - dp[i-1][j-1];
            }
        }
        for (int i=1;i<=m;i++) {
            for (int j=1;j<=n;j++) {
                int x0 = Math.max(i - k -1, 0);
                int y0 = Math.max(j - k -1, 0);
                int x1 = Math.min(i + k, m);
                int y1 = Math.min(j + k, n);
                result[i-1][j-1] = dp[x1][y1] - dp[x0][y1] - dp[x1][y0] + dp[x0][y0];
            }
        }
        return result;
    }

}