[原创]HDU 1559 最大子矩阵 [二维树状数组]【数据结构】

2016-11-09 15:09:25 Tabris_ 阅读数：265

博客爬取于2020-06-14 22:42:43
以下为正文

版权声明：本文为Tabris原创文章，未经博主允许不得私自转载。
https://blog.csdn.net/qq_33184171/article/details/53100103

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1559
---------------------------------------------------------------------------------------------.
最大子矩阵

Time Limit: 30000/10000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4361 Accepted Submission(s): 2273

Problem Description
给你一个m×n的整数矩阵，在上面找一个x×y的子矩阵，使子矩阵中所有元素的和最大。

Input
输入数据的第一行为一个正整数T，表示有T组测试数据。每一组测试数据的第一行为四个正整数m,n,x,y（0< m,n<1000 AND 0< x<=m AND 0< y<=n），表示给定的矩形有m行n列。接下来这个矩阵，有m行，每行有n个不大于1000的正整数。

Output
对于每组数据，输出一个整数，表示子矩阵的最大和。

Sample Input
1
4 5 2 2
3 361 649 676 588
992 762 156 993 169
662 34 638 89 543
525 165 254 809 280

Sample Output
2474

Author
lwg

Source
HDU 2006-12 Programming Contest

---------------------------------------------------------------------------------------------.

题目大意：中文题

解题思路 ：
根据这个时限 目测暴力都能过去 ，但是为了练习树状数组 还是用二维的树状数组过得

最最最最裸的二维树状数组 没什么可说的

附本题代码
----------------------------------------------------------------------------------------.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
# include <bits/stdc++.h>

# define abs(x)          (((x)>0)?(x):-(x))
# define lalal           puts("*********")
# define Rep(a,b,c)      for(int a=(b);a<=(c);a++)
# define Req(a,b,c)      for(int a=(b);a>=(c);a--)
# define Rop(a,b,c)      for(int a=(b);a<(c);a++)
# define s1(a)           scanf("%d",&a)
typedef long long int LL;
using namespace std;
/**************************************/
const int N = 1000+5;
# define lowbit(x)       (x&-x)
int cnt,sum[N][N],n,m;
void update(int xi,int yi,int val){
    for(int i=xi;i<=n;i+=lowbit(i))
        for(int j=yi;j<=m;j+=lowbit(j))
            sum[i][j]+=val;
    return;
}
int getSum(int xi,int yi){
    int ans = 0;
    for(int i=xi;i>0;i-=lowbit(i))
        for(int j=yi;j>0;j-=lowbit(j))
            ans+=sum[i][j];
    return ans;
}
int main(){
    int _;
    while(~s1(_)){
        while(_--){
            memset(sum,0,sizeof(sum));
            int x,y,temp;
            s1(n),s1(m),s1(x),s1(y);
            Rep(i,1,n)Rep(j,1,m) s1(temp),update(i,j,temp);
            int mx=-1;
            Rep(i,x,n)Rep(j,y,m){
                temp=getSum(i,j)-getSum(i,j-y)-getSum(i-x,j)+getSum(i-x,j-y);
                mx=max(mx,temp);
            }
            printf("%d\n",mx);
        }
    }
    return 0;
}