codevs 1373 射命丸文（矩阵前缀和）

题目描述 Description
在幻想乡，射命丸文是以偷拍闻名的鸦天狗。当然，文文的照相机可不止能够照相，还能够消除取景框里面所有的弹幕。假设现在文文面前有一块N行M列的弹幕群，每一个单位面积内有分值有num[i][j]的弹幕。相机的取景框可以将一块R行C列的弹幕消除，并且得到这一块区域内所有弹幕的分值(累加)。现在文文想要取得尽可能多的分值，请你计算出她最多能够得到的分值。

输入描述 Input Description
第1行：4个正整数N,M,R,C 第2..N+1行：每行M个正整数，第i+1行第j个数表示num[i][j]

输出描述 Output Description
第1行：1个整数，表示文文能够取得的最大得分

样例输入 Sample Input
3 5 2 3
5 2 7 1 1
5 9 5 1 5
3 5 1 5 3

样例输出 Sample Output
33

数据范围及提示 Data Size & Hint
对于60%的数据：1 <= N,M <= 200
对于100%的数据：1 <= N,M <= 1,000 1 <= R <= N, 1 <= C <= M 1 <= num[i][j] <= 1000
保证结果不超过2,000,000,000

题解：最暴力的做法就是枚举每一个子矩阵求和找最大值，这样显然会超时。我们可以用矩阵前缀和来做这道题。我们读入的时候顺便求一下矩阵前缀和，然后就可以直接用。那么我们怎么通过矩阵前缀和求子矩阵的和呢？其实跟前缀和差不多，详见下图：

矩阵4的和等于最大的矩阵减去（（矩阵1+矩阵2）+（矩阵1+矩阵3）- 矩阵1），矩阵前缀和就是这样233，然后我们就可以O(1)询问某个子矩阵的和，是不是很神奇╮(╯▽╰)╭。

代码如下：

#include<cstdio>#include<cstring>#include<iostream>using namespace std;int sum[1010][1010],map[1010][1010];int main(){    int n,m,r,c,maxx=0;    scanf("%d%d%d%d",&n,&m,&r,&c);    for(int i=1;i<=n;i++)        for(int j=1;j<=m;j++)        {            scanf("%d",&map[i][j]);            sum[i][j]=map[i][j]+sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];//求矩阵前缀和        }    for(int i=1;i<=n-r+1;i++)        for(int j=1;j<=m-c+1;j++)            if(sum[i+r-1][j+c-1]-sum[i+r-1][j-1]-sum[i-1][j+c-1]+sum[i-1][j-1]>maxx)//找子矩阵和的最大值                 maxx=sum[i+r-1][j+c-1]-sum[i+r-1][j-1]-sum[i-1][j+c-1]+sum[i-1][j-1];    printf("%d",maxx);    return 0;}