AcWing 4405. 统计子矩阵(双指针,前缀和)

给定一个 N × M N×M N×M 的矩阵 A A A,请你统计有多少个子矩阵 (最小 1 × 1 1×1 1×1,最大 N × M N×M N×M) 满足子矩阵中所有数的和不超过给定的整数 K K K?

输入格式

第一行包含三个整数 N , M N,M N,M 和 K K K。

之后 N N N 行每行包含 M M M 个整数,代表矩阵 A A A。

输出格式

一个整数代表答案。

数据范围

对于 30 30% 30 的数据, N , M ≤ 20 N,M≤20 N,M≤20,

对于 70 70% 70 的数据, N , M ≤ 100 N,M≤100 N,M≤100,

对于 100 100% 100 的数据, 1 ≤ N , M ≤ 500 ; 0 ≤ A i j ≤ 1000 ; 1 ≤ K ≤ 2.5 × 1 0 8 1≤N,M≤500;0≤Aij≤1000;1≤K≤2.5×10^8 1≤N,M≤500;0≤Aij≤1000;1≤K≤2.5×108。

输入样例:

3 4 10
1 2 3 4
5 6 7 8
9 10 11 12

输出样例:

19

样例解释

满足条件的子矩阵一共有 19,包含:

大小为 1×1 的有 10 个。
大小为 1×2 的有 3 个。
大小为 1×3 的有 2 个。
大小为 1×4 的有 1 个。
大小为 2×1 的有 3 个。

没想到同一道题,仅仅是隔了一天,再一次遇见就又不会写了,所以复习很重要。。

对于此题,时间限制为 1 s 1s 1s,而矩阵的行和列的范围为 0 0 0到 500 500 500,假如我们暴力枚举每一行每一列,就会有 500 500 500的四次方的时间复杂度,肯定爆了,但是假如只枚举行,就会有 500 500 500的平方的复杂度,再加上运用双指针来扫描列,时间复杂度就不会爆了。

那么可以在处理前缀和的时候仅处理列前缀和,也就是不再处理出整个矩阵的前缀和,因为我们要枚举行。

之后再使用双指针扫描列,比如使用 l l l和 r r r。

  • 首先定义 l l l和 r r r为 1 1 1(初始下标),然后每次循环以 r r r为主线,让 r r r自增
  • 在循环中,首先处理出来 1 1 1 ~ r r r列(在我们枚举的行的区间中)的总和,如果这个总和大于了$k,就让总和减去 l l l列的总和,然后让 l l l往右走(因为这时候如果r在比较靠左的位置就已经大于k了,如果r再往右走,就也一定会大于 k k k,那么就必须让 l l l往右走),并且在一次加和大于 k k k之后,每次 r r r往右走, l l l都肯定要往右走,这里不需要进行判断或是什么。

代码:

cpp 复制代码
#include<iostream>
using namespace std;
const int N = 510;
#define ll long long
ll s[N][N];
int n, m, k;
int main() {
	cin >> n >> m >> k;
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			cin >> s[i][j];
			s[i][j] += s[i - 1][j]; //处理列前缀和
		}
	}

	ll res = 0;
	for (int i = 1; i <= n; i++) {      //枚举行
		for (int j = i; j <= n; j++) {
			//注意这里sum定义到for里,因为在双指针扫描的过程中sum的值不能被初始化
			for (int l = 1, r = 1, sum = 0; r <= m; r++) {  
				sum += s[j][r] - s[i-1][r]; //先加上r列的总和
				while (sum > k) {           //如果超过了
					sum -= s[j][l] - s[i - 1][l];
					l++;   //在sum > k之后,每次r走,j必定走
				}

				res += r - l + 1;   //矩阵数量
			}
		}
	}
	cout << res;
	return 0;
}
相关推荐
retaw_027 分钟前
74. 搜索二维矩阵
线性代数·矩阵
知识鱼丸38 分钟前
machine learning knn算法之使用KNN对鸢尾花数据集进行分类
算法·机器学习·分类
周杰伦_Jay40 分钟前
简洁明了:介绍大模型的基本概念(大模型和小模型、模型分类、发展历程、泛化和微调)
人工智能·算法·机器学习·生成对抗网络·分类·数据挖掘·transformer
凭君语未可1 小时前
豆包MarsCode:小C点菜问题
算法
C语言魔术师1 小时前
【小游戏篇】三子棋游戏
前端·算法·游戏
自由自在的小Bird1 小时前
简单排序算法
数据结构·算法·排序算法
王老师青少年编程7 小时前
gesp(C++五级)(14)洛谷:B4071:[GESP202412 五级] 武器强化
开发语言·c++·算法·gesp·csp·信奥赛
DogDaoDao8 小时前
leetcode 面试经典 150 题:有效的括号
c++·算法·leetcode·面试··stack·有效的括号
Coovally AI模型快速验证8 小时前
MMYOLO:打破单一模式限制,多模态目标检测的革命性突破!
人工智能·算法·yolo·目标检测·机器学习·计算机视觉·目标跟踪
可为测控9 小时前
图像处理基础(4):高斯滤波器详解
人工智能·算法·计算机视觉