hdu 1559【最大子矩阵和DP】

来源：互联网发布：php权限修饰符编辑：程序博客网时间：2024/06/05 05:47

这题比poj那道多了限制，就是规定了子矩阵大小，开始看到时间10s，果断暴力，接着又果断TLE，加一点点优化还是TLE，坑爹！

接着就用DP，怎么DP？很简单，令a[i][j]为1<=s<=i，1<=t<=j所有元素的和，可以在线处理，那么就可以通过加减来求出当前子矩阵和，就好像计算几个矩阵面积一样（某部分覆盖在一起，所以相加后减去重复的面积）。

#include <vector>#include <list>#include <map>#include <set>#include <queue>#include <string.h>#include <deque>#include <stack>#include <bitset>#include <algorithm>#include <functional>#include <numeric>#include <utility>#include <sstream>#include <iostream>#include <iomanip>#include <cstdio>#include <cmath>#include <cstdlib>#include <limits.h>using namespace std;int lowbit(int t){return t&(-t);}int countbit(int t){return (t==0)?0:(1+countbit(t&(t-1)));}int gcd(int a,int b){return (b==0)?a:gcd(b,a%b);}int max(int a,int b){return a>b?a:b;}int min(int a,int b){return a>b?b:a;}#define LL long long#define PI acos(-1.0)#define N  1010#define MAX INT_MAX#define MIN INT_MIN#define eps 1e-8#define FRE freopen("a.txt","r",stdin)int n,m,x,y;int a[N][N];int main(){    int t;    scanf("%d",&t);    while(t--)    {        scanf("%d%d%d%d",&n,&m,&x,&y);        int i,j;        int ans=-1;        memset(a,0,sizeof(a));        for(i=1;i<=n;i++)        {            for(j=1;j<=m;j++)            {                scanf("%d",&a[i][j]);                a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];                if(i>=x && j>=y)                {                    int tmp=a[i][j]-a[i-x][j]-a[i][j-y]+a[i-x][j-y];                    if(tmp>ans)ans=tmp;                }            }        }        printf("%d\n",ans);    }    return 0;}