HDU - 1559 - 最大子矩阵

先上题目：

最大子矩阵

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

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

 

　　中文题意不解释，这题和HDU的1003 1081是一类的题，方法一样，枚举区间的值，然后找出最大值。

　　用两个坐标i,j在其中一维上面移动，保持区间长度为要求的值，然后求这个区间里面要求大小的矩阵的和，这里需要先以另一维为方向预处理，预处理的效果是得出前k项元素的和。

　　详细做法可参考：

　　http://www.cnblogs.com/sineatos/p/3515309.html

　　http://www.cnblogs.com/sineatos/p/3560325.html

 

上代码：

 

#include <cstdio>
#include <cstring>
#define MAX 1002
#define max(x,y) (x > y ? x : y)
using namespace std;

int s[MAX][MAX];
int cy[MAX][MAX];

int main()
{
    int t;
    int n,m,x,y,maxn,sum;
    //freopen("data.txt","r",stdin);
    scanf("%d",&t);
    while(t--){
        scanf("%d %d %d %d",&n,&m,&y,&x);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&s[i][j]);
            }
        }
        memset(cy,0,sizeof(cy));
        for(int j=1;j<=m;j++){
            for(int i=1;i<=n;i++){
                cy[i][j]+=s[i][j]+cy[i-1][j];
            }
        }
        maxn=0;
        int i,j;
        for(i=1,j=y;j<=n;i++,j++){
            sum=0;
            for(int k=1;k<=m;k++){
                sum+=cy[j][k]-cy[i-1][k];
                if(k>x){
                    sum-=cy[j][k-x]-cy[i-1][k-x];
                }
                maxn=max(sum,maxn);
            }
        }
        printf("%d\n",maxn);
    }
    return 0;
}