最大矩形問題總結

昨天晚上刷水題寫到一道最大全1正方形，於是興致勃勃將系列相關都復習了一遍qwq

1.最大子段和 

測試地址：https://www.luogu.org/problem/show?pid=1115

Problem：給出一段序列，選出其中連續且非空的一段使得這段和最大。

Key：f[i]表示以i為終點時的最大子段和，f[i]=max(f[i-1],0)+a[i]；ans=max(ans,f[i])；

#include<cmath>
#include<algorithm>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn=200009;
int a[maxn],f[maxn];
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
    }
    for(int i=1;i<=n;i++)
    {
        f[i]=max(f[i-1],0)+a[i];
    }
    int ans=-0x3f3f3f3f;
    for(int i=1;i<=n;i++)
    {
        ans=max(ans,f[i]);
    }
    cout<<ans<<endl;
    return 0;
}

 

2.最大子矩形和（最大加權矩形） 

測試地址：https://www.luogu.org/problem/show?pid=1719

Problem：給定一個整數矩陣，求其中一個任意大小的子矩形，使該子矩形內值的和最大。

Key：將二維的矩陣按行壓成一維，枚舉列來做最大子段和，復雜度O(n^3)。

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=150;
int a[maxn][maxn],n,m,sum[maxn];
int ans=-0x3f3f3f3f;
int f[maxn];
int calc_max(int val[],int n)
{
    int ret=-0x3f3f3f3f;
    f[1]=val[1];
    for(int i=2;i<=n;i++)
        f[i]=max(f[i-1],0)+val[i],ret=max(ret,f[i]);
    return ret;
}

int main()
{
    scanf("%d",&n);
    m=n;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int k=1;k<=m;k++)
        {
            sum[k]=0;
        }
        for(int j=i;j<=n;j++)
        {
            int k;
            for(k=1;k<=m;k++)
            {
                sum[k]+=a[j][k];
            }
            int mx=calc_max(sum,k);
            ans=max(ans,mx);
        }
    }
    cout<<ans<<endl;
    return 0;
}

 

3.最大全1正方形

測試地址：https://www.luogu.org/problem/show?pid=1387

Problem：在一個n*m的只包含0和1的矩陣里找出一個不包含0的最大正方形，輸出邊長。

Key：

　　解法1：暴力維護二維前綴和，枚舉邊長來更新ans，其中求二維前綴和是容斥原理。

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=109;
int n,m,sum[maxn][maxn],ans=0;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int x;
            scanf("%d",&x);
            sum[i][j]=sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1]+x;
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            for(int l=1;l+i<=n&&l+j<=m;l++)
            {
                int ok=sum[i+l][j+l]-sum[i+l][j]-sum[i][j+l]+sum[i][j];
                if(ok==l*l)
                {
                    ans=max(ans,l);
                }
            }
        }
    }
    printf("%d",ans);
    return 0;
}

　　解法2：f[i][j]表示以（i,j）為右下角的時候，能延伸左上角到最大的的正方形。如果（i,j）這個點為0時直接不考慮，不為0時才計算轉移方程f[i,j]=min(f[i-1,j-1],f[i-1,j],f[i,j-1])，轉移方程的實際意義是右上方、左方、右方能擴展的最大邊長+1 （加1的操作就是計算（i,j）本身）。

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=109;
int f[maxn][maxn],n,m,ans=0;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++)
    {
        int x;
        scanf("%d",&x);
        if(x)//如果這一個格子是1
        {//f[i][j]表示以(i,j)為右下角的，能延伸左上角最大的（正方形）邊長
            f[i][j]=min(f[i-1][j-1],min(f[i-1][j],f[i][j-1]))+1;
            //取min是為了驗證是不是三個格子都能達到擴展要求
            //如果其中一個值為0，說明那個角沒法擴展正方形
            //如果其中最小值為x，說明這個點最多擴展x，加他自己就是+1
            ans=max(ans,f[i][j]);
        }
    }
    cout<<ans<<endl;
    return 0;
}

 

4.最大全1矩形

測試地址：https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1158

Problem：給出1個M*N的矩陣M1，里面的元素只有0或1，找出M1的一個子矩陣M2，M2中的元素只有1，並且M2的面積是最大的。輸出M2的面積。

Key：

　　解法1：存儲矩陣時，將非1的元素存成-inf，直接做一遍最大子矩形和。

　　解法2：處理每一行的每一個元素能向下延伸的最大全1長度，這樣得到每一行都是一個可以“histogram求最大矩形”的序列。主要做法是對於每一行，從左到右遍歷一遍剛剛得到的序列，用f[i]表示以i為起點時，向右能擴展的最大矩形面積，最后在所有的f數組中取max。

　　解法3：遞推處理，f[i][j]表示第i行，第j列時，從i能延伸的最大全1長度。記錄l和r兩個數組，分別表示能向左或右擴展的最大位置。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
const int maxn=1009;
int l[maxn],r[maxn];
int f[maxn][maxn],n,m;
int main()
{
    
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j)
        {
            int x;
            scanf("%d",&x);
            if(x==1)f[i][j]=f[i-1][j]+1;
            else f[i][j]=0;
        }
    int ans=0;
    for(int i=1; i<=n; ++i)
    {
        for(int j=1; j<=m; ++j)
            l[j]=r[j]=j;
        f[i][0]=f[i][m+1]=-1;
        for(int j=1; j<=m; ++j)
        {
            while(f[i][j]<=f[i][l[j]-1])
                l[j]=l[l[j]-1];
        }
        for(int j=m; j>=1; --j)
        {
            while(f[i][j]<=f[i][r[j]+1])
                r[j]=r[r[j]+1];
        }
        for(int j=1; j<=m; ++j)
        {
            if(f[i][j])
            {
                int t=(r[j]-l[j]+1)*f[i][j];
                if(t>ans)ans=t;
            }
        }
    }
    printf("%d\n",ans);
    return 0;
}