hdu 1081 最大子矩阵(二维前缀和优化）@

最新推荐文章于 2022-12-13 14:51:18 发布

mrcoderrev

最新推荐文章于 2022-12-13 14:51:18 发布

阅读量932

点赞数

分类专栏：思维

本文链接：https://blog.csdn.net/yjf3151731373/article/details/53325662

版权

思维专栏收录该内容

59 篇文章 0 订阅

订阅专栏

最大子矩阵

Time Limit:10000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u

Submit Status

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

暴力T了，前缀和优化。。。

输入挂快的一笔。。。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 1100;
int s[N][N], sum[N][N];
int m, n;
const int BufferSize=1<<16;
char buffer[BufferSize],*head,*tail;
inline char Getchar() {
    if(head==tail) {
        int l=fread(buffer,1,BufferSize,stdin);
        tail=(head=buffer)+l;
    }
    return *head++;
}
inline int read() {
    int x=0,f=1;char c=Getchar();
    for(;!isdigit(c);c=Getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=Getchar()) x=x*10+c-'0';
    return x*f;
}


LL judge(int l, int r,int x, int y)
{
    LL ans=0;
    for(int i=l; i<l+x; i++)
    {
        for(int j=r; j<r+y; j++)
        {
            ans+=s[i][j];
        }
    }
    return ans;
}
LL maxt(LL l,LL r)
{
    return l>r?l:r;
}


int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        int x, y;
        memset(sum,0,sizeof(sum));
        //scanf("%d %d %d %d", &m, &n, &x, &y);
        m=read(),n=read(),x=read(),y=read();
        for(int i=1; i<=m; i++)
        {
            for(int j=1; j<=n; j++)
            {
                s[i][j]=read();
                sum[i][j]=sum[i][j-1]+s[i][j];
            }
            for(int j=1; j<=n; j++) sum[i][j]+=sum[i-1][j];
        }
        LL ans=0;
        for(int i=1; i<=m-x+1; i++)
        {
            for(int j=1; j<=n-y+1; j++)
                ans=maxt(ans,sum[i+x-1][j+y-1]+sum[i-1][j-1]-sum[i+x-1][j-1]-sum[i-1][j+y-1]);
        }
        printf("%I64d\n",ans);
    }
    return 0;
}