P1719 最大加权矩形 矩形压缩，前缀和

题目描述

为了更好的备战 NOIP2013 ，电脑组的几个女孩子 LYQ,ZSC,ZHQ 认为，我们不光需要机房，我们还需要运动，于是就决定找校长申请一块电脑组的课余运动场地，听说她们都是电脑组的高手，校长没有马上答应他们，而是先给她们出了一道数学题，并且告诉她们：你们能获得的运动场地的面积就是你们能找到的这个最大的数字。

校长先给他们一个 N×N矩阵。要求矩阵中最大加权矩形，即矩阵的每一个元素都有一权值，权值定义在整数集上。从中找一矩形，矩形大小无限制，是其中包含的所有元素的和最大 。矩阵的每个元素属于 [−127,127] ,例如

 0 –2 –7  0 
 9  2 –6  2
-4  1 –4  1 
-1  8  0 –2

在左下角：

9  2
-4  1
-1  8

和为15。

几个女孩子有点犯难了，于是就找到了电脑组精打细算的 HZH，TZY 小朋友帮忙计算，但是遗憾的是他们的答案都不一样，涉及土地的事情我们可不能含糊，你能帮忙计算出校长所给的矩形中加权和最大的矩形吗？

输入格式

第一行：n ，接下来是 n行 n 列的矩阵。

输出格式

最大矩形（子矩阵）的和。

输入输出样例

输入 #1

4
0 -2 -7 0
 9 2 -6 2
-4 1 -4  1 
-1 8  0 -2

输出 #1

15

说明/提示

n≤120

思路：

一维版本：给出一段序列，选出其中连续且非空的一段使序列和最大。

状态转移方程： dp[i]=max(temp, dp[i-1]+temp) temp表示这个数列的·第i项

此处为二维，那么我们想办法将二维压缩为一维，利用前缀和实现矩形压缩，

即在输入时加上a[i-1],那么a[ i ][ j ] - a[ i - k ][ j ]正好是以i为最下面一行，往上k行的压缩结果，这就很方便地表示了压缩后的矩形。然后我们用循环for(int i=1;i<=n;i++){

for(int k=1;k<=i;k++)

}就可以找出各行为最下一行，不同行数的矩阵最大值。

代码如下：

#include <bits/stdc++.h>
using namespace std;
int n,ans,a[200][200],dp[200];
int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            cin>>a[i][j];
            a[i][j]+=a[i-1][j];//前缀和 
        }
    }
    for(int i=1;i<=n;i++){
        for(int k=1;k<=i;k++){
            memset(dp,0,sizeof(dp));
            for(int j=1;j<=n;j++){
                dp[j]=a[i][j]-a[i-k][j];//相减求压缩后的矩形第j列的值
                dp[j]=max(dp[j-1]+dp[j],dp[j]);//更新dp值 ：权值和 
                // 一维版本：dp[i]=max(dp[i-1]+tmp,tmp)     ，此处是压缩成了一维 
                ans=max(dp[j],ans); //更新最大值 
            }
        }
    }
    cout<<ans<<endl;
    return 0;
}