最大子矩阵——杨子曰算法

最大子矩阵——杨子曰算法

来，今天曰个模板——最大子矩阵
简单的说就是给你一个矩阵（n*m，有负数），让你在这个矩阵中取出一个子矩阵，是得子矩阵的和最大。

1.O(n^6)巨型大暴力
n^ 2枚举子矩阵左上角，n^ 2枚举子矩阵右下角，再n^ 2求个和，O(n^6)搞定，如果你想到的是这个，杨子只想曰：呵呵

2.O(n^4)机智大暴力
依然用n^ 2枚举子矩阵左上角，n^ 2枚举子矩阵右下角，但求和只用O(1),Why？——割补法（那是神马，戳→二位前缀和(矩阵割补法)）,于是O(n^4)搞定，能想到这个，不要得瑟，只是正常罢了，与正解的思路完全不同

3.O(n^3)正解走起

---

我相信99%打开这个Blog的人都已经是掌握了最大子序列了，如果你是那1%——请自觉关掉这个页面（我估计好奇心驱使着你也会往下看）

如果你，知道了最大子序列的做法，那它的复杂的是O(n)，那我们能不能，把一个二维数组压成一维能？，举个例子：

这个矩阵，我们可以把它压成3个一维序列，也就是分别求出（a[1][1],a[1]][2],a[1][3]）,(a[2][1],a[2][2],a[2][3]),(a[1][1]+a[2][1],a[1][2]+a[2][2],a[1][3]+a[2][3])的最大子序列，再取一个max，完事
So，总结一下，我们可以枚举上限和下限，把a[i][k]到a[j][k]的和存在f[k]中，在对f数组求一个最大子序列，欸欸欸，求和不又要一个O(n)吗?
So,我们是用前缀和，把序列的列存个前缀和
O(n^ 2)枚举上下限，O(n)求最大子序列。O(n^3)!!!
OK,完事

---

破天荒pascal代码：

var
    a:array[0..300,0..300]of longint;
    i,j,n,res,k,max,m,maxx:longint;
begin
    max:=-maxlongint;
    read(n,m);
    for i:=1 to n do
        for j:=1 to m do
            read(a[i,j]);
    for i:=1 to n do
        a[i,0]:=0;
    for j:=1 to m do
        a[0,i]:=0;
    for i:=1 to n do
        for j:=1 to m do
            a[i,j]:=a[i,j]+a[i-1,j];
    for i:=0 to n-1 do
        for j:=i+1 to n do begin
            res:=0;
            maxx:=-maxlongint;
            for k:=1 to m do  begin
                if res<=0
                    then res:=a[j,k]-a[i,k]
                    else res:=res+a[j,k]-a[i,k];
                if res>maxx then maxx:=res;
        end;
        if maxx>max
            then max:=maxx;
        end;
    write(max);
end.

于TJQ高层小区
未经作者允许，严禁转载：https://blog.csdn.net/HenryYang2018/article/details/79768658