[模板]树状数组(以hdu2352 Stars(一维)/hdu1195 Mobile phones(二维)为例

本文深入讲解了一维和二维树状数组的应用,通过具体实例解析树状数组如何高效地解决区间和查询问题。针对不同场景,文章提供了清晰的代码实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

树状数组就是巧妙的利用了二分,虽然代码很短,但是很奇妙。

一维树状数组

学习链接,模板以hdu2352 Stars为例

题意:
每颗星星都有坐标,若一个星星坐标为(x,y),记录坐标为(_x,_y)(_x<=x && _y<=y && (_x!=x && _y!=y))的星星数目,也就是除了它本身的左下角的星星啦,这个数目就是星星的等级,求(0-(N-1))等级星星的数目;
星星数目 N (1<=N<=15000),接下来N行坐标(X,Y)( 0<=X,Y<=32000);
坐标被以Y坐标的上升序列给出,Y坐标相等的以X坐标的上升序列给出。
思路:
题目其实很明显的提供了思路,“坐标被以Y坐标的上升序列给出,Y坐标相等的以X坐标的上升序列给出”,根据这句话我们可以“降维”,tree[]数组只需要记录X坐标就好,因为此刻的星星,Y坐标必然>=前面的星星,只需要记录X坐标>=前面的星星数目就好;
还要特别注意一点,坐标为(0,Y)的情况,树状数组不要出现(x=0)这种情况,我们可以人为使x+=1。
代码:

#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
/*
树状数组就是巧妙的利用了二分
k为i的末尾0的个数
*/
#define lowbit(x) (x&(-x))
const int N = 32000+5;
int tree[N];//树状数组
int f[N];//记录星星的X坐标
//向上更新
void update(int k, int num){
    while(k <= N){
        tree[k] += num;
        k += lowbit(k);
    }
}
//向下统计
int sum(int k){
    int re = 0;
    while(k>0){
        re += tree[k];
        k -= lowbit(k);
    }
    return re;
}
int main(){
    int n, a, b;
    while(~scanf("%d", &n)){
        memset(tree, 0, sizeof(tree));
        memset(f, 0, sizeof(f));
        for(int i=0; i<n; i++){
            scanf("%d%d", &a, &b);
            f[sum(a+1)]++;//记录不同等级的星星个数
            update(a+1,1);//向上更新
        }
        for(int i=0; i<n; i++)
            printf("%d\n", f[i]);
    }
    return 0;
}


二维树状数组

模板以hdu1195 Mobile phones为例

题意:
开arr[]代表基站,有0、1、2、3四个命令,分别代表的意思是
指令 参数 说明
0 S initialize初始化规格为S*S的矩阵为0,只会出现在第一行,且只出现一次
1 X Y A add arr[X][Y]+=A(A可以为负数,但题目有说我们可以默认此操作不会使arr数组出现负值)
2 L B R T query 以arr[L][B]为左上角arr[R][T]为右上角的矩阵的数值和
3 end 结束指令,只会出现在最后一行,且只出现一次
矩阵规格S*S: 1 * 1 <= S * S <= 1024 * 1024 
各基站的手机数目A: 0 <= V <= 32767 
变化的手机数目: -32768 <= A <= 32767 
指令数目: 3 <= U <= 60002 
整个矩阵所有基站的手机数目: M= 2^30
思路:
二维树状数组,是一维的扩展,与一维不同的是,一维管的是一段区间,二维管的是一块矩阵;
数组a[x][y]的x、y是负数不要紧,默认a[负值][负值]=0
代码:

#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;

#define lowbit(x) (x&(-x))
const int N = 1500;
int tree[N][N];
int s;
//向上更新
void update(int x, int y, int val){
    for(int i=x; i<=s; i+=lowbit(i))
    for(int j=y; j<=s; j+=lowbit(j))
        tree[i][j] += val;
}
int query_sum(int x, int y){
    int re = 0;
    for(int i=x; i>0; i-=lowbit(i))
    for(int j=y; j>0; j-=lowbit(j))
        re += tree[i][j];
    return re;
}
int main(){
    int order;
    int x, y, a;
    int l, b, r, t;

    scanf("%d%d", &order, &s);
    memset(tree, 0, sizeof(tree));
    while(scanf("%d", &order) && order!=3)
    switch(order){
        case 1:{
            scanf("%d%d%d", &x, &y, &a);
            update(x+1,y+1,a);
            break;
        }
        case 2:{
            scanf("%d%d%d%d", &l, &b, &r, &t);
            printf("%d\n", query_sum(r+1,t+1)+query_sum(l,b)-query_sum(r+1,b)-query_sum(l,t+1));
            break;
        }
    }
    return 0;
}


反思:

这道题有毒,还是注意坐标有可能有0,因为这个TLE了一次,附上discuss找的测试数据,谢谢这位叫 lanjiongjiang 的同学

/*
0 1
1 0 0 28766
2 0 0 0 0
1 0 0 -2826
2 0 0 0 0
1 0 0 -587
2 0 0 0 0
1 0 0 1017
2 0 0 0 0
1 0 0 -17234
2 0 0 0 0
1 0 0 19618
2 0 0 0 0
1 0 0 -2237
2 0 0 0 0
1 0 0 -9192
2 0 0 0 0
1 0 0 -11132
2 0 0 0 0
1 0 0 1842
2 0 0 0 0
3
--------------
28766
25940
25353
26370
9136
28754
26517
17325
6193
8035
*/


<think>我们注意到用户要求解决HDU1081问题,但实际提供的引用中并没有直接关于HDU1081的内容。HDU1081题目名是"ToTheMax",是一个最大子矩阵和问题,而不是二维最长递增子序列问题。这里可能存在用户描述的错误。然而,根据我们的知识,HDU1081(ToTheMax)的题意:给定一个二维矩阵(可能包含正数和负数),求子矩阵的最大和。这是一个经典的最大子矩阵和问题,通常使用动态规划来解决,可以通过转换为一维的最大子段和问题来求解。解决思路:1.将二维问题转化为一维:对于每一列,计算从第i行到第j行之间的每一列的元素和,形成一个一维数组。2.对形成的一维数组,求最大子段和(即连续子数组的最大和)。3.遍历所有可能的i和j(行),取最大子段和的最大值。因此,我们需要实现一个最大子矩阵和算法。但是,用户提到了“二维最长递增子序列”,而实际上HDU1081是最大子矩阵和问题。这里我们按照正确题意给出代码。为了确认,我们参考引用:引用[1]:排列型枚举,与问题无关。引用[2]:最长上升子序列(LIS),一维情况,采用二分优化。引用[3]:最长公共子序列(LCS)。用户的问题描述为“二维最长递增子序列”,但根据题号HDU1081,实际应为最大子矩阵和问题。我们按照最大子矩阵和来给出Java实现。注意:如果用户坚持要求“二维最长递增子序列”的代码,那么需要澄清,因为通常“最长递增子序列”是一维的,二维情况下可以转化为偏序问题(如按一维排序,另一维求最长上升子序列),但题目HDU1081并不是。因此,我们按照HDU1081的正确题意(最大子矩阵和)给出Java代码。代码步骤:1.读取一个n*n的矩阵。2.用两层循环枚举起始行i和终止行j。3.对于固定的i和j,计算每一列k从第i行到第j行的元素和,保存到一个一维数组colSum中(注意:对于每一列,我们只需要累加从i到j行的元素)。4.对colSum数组求最大子段和。5.更新最大子矩阵和。最大子段和算法(Kadane算法):intmaxEndingHere=0;intmaxSoFar=Integer.MIN_VALUE;for(intk=0;k< n;k++){maxEndingHere=Math.max(0,maxEndingHere+colSum[k]);//如果允许负数,则可能使用Math.max(colSum[k],maxEndingHere+colSum[k]),但题目要求子矩阵可以是任意子矩阵,包括负数,所以正确写法如下://但注意,如果矩阵中有负数,那么初始值不能设为0,因为可能整个都是负数。所以Kadane算法标准写法://初始化://intmaxEndingHere=colSum[0];//intmaxSoFar=colSum[0];//for(intk=1;k<n;k++){//maxEndingHere=Math.max(colSum[k],maxEndingHere+colSum[k]);//maxSoFar=Math.max(maxSoFar,maxEndingHere);//}}但是,我们的colSum数组是通过多行累加得到的,所以可能包含负数。因此,必须用标准的Kadane算法。具体Java代码实现如下:</think> ### HDU 1081 "To The Max" 最大子矩阵和问题(Java实现) 虽然您提到“二维最长递增子序列”,但根据题号HDU 1081,实际是求最大子矩阵和问题(给定包含正负整数的矩阵,求元素和最大的子矩阵)。以下是完整的Java实现方案: ```java import java.util.Scanner; public class Main { public static void main(String[] args) { Scanner sc = new Scanner(System.in); while (sc.hasNext()) { int n = sc.nextInt(); int[][] matrix = new int[n][n]; // 输入矩阵 for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { matrix[i][j] = sc.nextInt(); } } System.out.println(maxSubmatrixSum(matrix)); } sc.close(); } // 动态规划求最大子矩阵和 private static int maxSubmatrixSum(int[][] matrix) { int n = matrix.length; int maxSum = Integer.MIN_VALUE; // 枚举起始行 for (int top = 0; top < n; top++) { int[] colSum = new int[n]; // 存储列累加值 // 枚举结束行 for (int bottom = top; bottom < n; bottom++) { // 累加当前行各列值 for (int j = 0; j < n; j++) { colSum[j] += matrix[bottom][j]; } // 对当前列累加数组求最大子序和 int current = colSum[0]; int maxSeq = colSum[0]; for (int j = 1; j < n; j++) { current = Math.max(colSum[j], current + colSum[j]); maxSeq = Math.max(maxSeq, current); } // 更新全局最大值 maxSum = Math.max(maxSum, maxSeq); } } return maxSum; } } ``` ### 算法说明 1. **核心思路**: - 将二维问题转换为一维:枚举所有行组合`(top, bottom)` - 将每列在`top`到`bottom`行间的元素压缩成一维数组 - 对压缩后的一维数组使用Kadane算法求最大子序和[^1] 2. **时间复杂度**: - 压缩行组合:$O(n^2)$ - 一维Kadane算法:$O(n)$ - 总时间复杂度:$O(n^3)$(适合$n \leq 100$的题目规模) 3. **空间优化**: - 使用固定长度数组`colSum[n]`存储列压缩结果 - 空间复杂度$O(n)$ ### 测试示 输入矩阵: ``` 4 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2 ``` 输出结果: ``` 15 ``` (最大子矩阵为左上角4×1矩阵:0 + 9 + (-4) + (-1) 计算有误,实际最大子矩阵是第二行的[9, 2]和第四行的[-1, 8]部分区域,但示输出应为15,来自子矩阵: ``` 9 2 -4 1 -1 8 ``` 的和:9+2+(-4)+1+(-1)+8=15) ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值