
本文介绍如何通过预计算二维前缀和(或局部区域和)优化子矩阵搜索,将暴力匹配的 o(rl·cl·rs·cs) 时间复杂度显著降低,在 hackerrank 等限时场景下稳定通过 1000×1000 规模测试用例。
本文介绍如何通过预计算二维前缀和(或局部区域和)优化子矩阵搜索,将暴力匹配的 o(rl·cl·rs·cs) 时间复杂度显著降低,在 hackerrank 等限时场景下稳定通过 1000×1000 规模测试用例。
在解决类似 HackerRank《The Grid Search》这类二维子矩阵匹配问题时,原始的四重嵌套暴力扫描(对每个合法左上角位置完整比对子矩阵)极易超时——尤其当大矩阵达 1000×1000、小矩阵达 100×100 时,最坏时间复杂度高达 $O(10^6 \times 10^4) = O(10^{10})$,远超 4 秒限制。
核心优化思路是:引入“区域和过滤”机制,将全量比对降为稀疏验证。具体而言,我们预先构建一个 sumGrid 二维数组,其中 sumGrid[i][j] 表示以 (i, j) 为右下角、尺寸与小矩阵完全一致(即 RS × CS)的子区域元素和。若该和不等于小矩阵总和,则该位置必然不匹配;仅当和相等时,才启动精确比对。这能大幅剪枝无效候选位置。
✅ 关键实现步骤
- 高效读取输入:避免 strings.Split(s, "") 拆分单字符字符串(生成大量短字符串对象),改用 for _, r := range s 遍历 rune,并直接减 '0' 转数字;
- 预计算小矩阵总和:一次性求和 targetSum;
-
构建区域和网格:
- 使用动态规划思想,利用二维滑动窗口技巧:先计算首行/首列的初始块和,再通过“移出左列 + 移入右列”方式逐列更新,避免重复求和;
- 或更简洁地:对每个合法右下角 (i, j)(满足 i ≥ RS-1 && j ≥ CS-1),用嵌套循环计算 RS×CS 子区域和(虽为 O(RS·CS),但仅执行 O(RL·CL) 次,总体仍优于暴力);
- 过滤 + 精确验证:遍历所有 sumGrid[i][j] == targetSum 的位置,调用 matchAt(i, j) 进行逐元素校验。
以下是优化后的 Go 实现(关键部分):
func solve() {
var t int
fmt.Scanf("%d", &t)
for ; t > 0; t-- {
var rL, cL, rS, cS int
fmt.Scanf("%d %d", &rL, &cL)
// 读取大矩阵:避免 strings.Split,直接解析数字字符
large := make([][]int, rL)
for i := range large {
large[i] = make([]int, cL)
var s string
fmt.Scanf("%s", &s)
for j, ch := range s {
large[i][j] = int(ch - '0')
}
}
fmt.Scanf("%d %d", &rS, &cS)
small := make([][]int, rS)
targetSum := 0
for i := range small {
small[i] = make([]int, cS)
var s string
fmt.Scanf("%s", &s)
for j, ch := range s {
val := int(ch - '0')
small[i][j] = val
targetSum += val
}
}
// 构建区域和网格 sumGrid[i][j] = 以 (i,j) 为右下角的 rS×cS 子矩阵和
sumGrid := make([][]int, rL)
for i := range sumGrid {
sumGrid[i] = make([]int, cL)
}
// 初始化:计算左上角首个 rS×cS 块
if rL >= rS && cL >= cS {
for i := 0; i <h3>⚠️ 注意事项与进阶建议</h3>
- 内存友好性:上述滑动更新法仅需 O(1) 额外空间(除存储矩阵本身),比维护完整二维前缀和数组更省内存;
- 边界处理:务必检查 rL >= rS && cL >= cS,否则跳过区域和计算;
- 精度保障:区域和相等是必要非充分条件(如 [[1,0],[0,1]] 和 [[0,1],[1,0]] 和相同但内容不同),必须保留最终逐元素验证;
- 进一步优化:对极端情况(如小矩阵极小),可结合 Rabin-Karp 思想对每行做哈希,再对行哈希序列做二维哈希,实现接近 O(RL·CL) 的期望时间复杂度;
- Go 特定陷阱:避免 strconv.Atoi 在高频循环中分配,字符转数字请用 ch - '0';切片预分配容量(如 make([]int, cL))减少扩容开销。
通过此方案,原超时代码在最大规模数据下运行时间可从 >10s 降至











