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;
}
相关推荐
AI科技星8 分钟前
三维网格—素数对偶性及其严格证明(全域数学·统一基态演化版)
算法·数学建模·数据挖掘
诸葛务农20 分钟前
光电对抗:多模复合制导烟雾干扰外场试验及仿真(4)
人工智能·算法·光电对抗
WolfGang00732128 分钟前
代码随想录算法训练营 Day39 | 动态规划 part12
算法·动态规划
阿Y加油吧41 分钟前
动态规划经典题解:最长递增子序列 & 乘积最大子数组
算法·动态规划·代理模式
f3iiish1 小时前
3783. 整数的镜像距离 力扣
算法·leetcode
Not Dr.Wang4221 小时前
基于matlab的控制系统奈氏图及其稳定性分析
数据结构·算法·matlab
闻缺陷则喜何志丹1 小时前
【排序 离散化 二维前缀和】 P7149 [USACO20DEC] Rectangular Pasture S|普及+
c++·算法·排序·离散化·二维前缀和
rainbow7242441 小时前
AI学习路线分享:通用型认证与算法认证学习体验对比
人工智能·学习·算法
君义_noip1 小时前
信息学奥赛一本通 4163:【GESP2512七级】城市规划 | 洛谷 P14921 [GESP202512 七级] 城市规划
c++·算法·图论·gesp·信息学奥赛
Simon_lca1 小时前
验厂不翻车!Acushnet 11 项核心政策 + 自查要点,一文搞定
大数据·人工智能·经验分享·算法·制造