CF 376D. Maximum Submatrix 2 枚举+前缀和

题意:n*m(01)矩阵,操作:交换任意两行. 现在可以操作任意次.问能得到的最大全1子矩阵size为多少? n,m<=5e3.


操作任意次 也就是将所有行重新排列一次.也可以将最后的行看作(1,2..n)的一个置换.
则合法子矩阵h*w 其h不一定是连续的 但是每行的1~w都为1.


暴力的做法:枚举结尾列,在枚举宽度,在O(n)确定该宽度所能得到的高度. O(n^3).

同样枚举结尾列 不枚举宽度,O(n)枚举行,通过预处理出,d[i][j],第i行第j列结尾的最大宽度x 则第i行j列u结尾的宽度为[1..x]

设置f(j)[w]:表示第j列结尾宽度为w的有多少行. 将f[d[i][j]]++,最后从后往前累加一次即可. O(n*m).

#include <bits/stdc++.h>
using namespace std;
const int N=5e3+5;
char s[N];
int n,m,d[N][N],a[N][N],f[N];
int main()
{
    cin>>n>>m;

    for(int i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        for(int j=1;j<=m;j++)
        {
            a[i][j]=s[j]-'0';
            if(a[i][j])
                d[i][j]=d[i][j-1]+1;
            else
                d[i][j]=0;
        }
    }
    int ans=0;
    for(int j=1;j<=m;j++)
    {
        memset(f,0,sizeof(f));
        for(int i=1;i<=n;i++)
            f[d[i][j]]++;//f[w]:how many rows which finish at col j with width w
        for(int i=m;i>=1;i--)
            f[i]+=f[i+1];

        for(int i=1;i<=m;i++)
            ans=max(ans,f[i]*i);
    }
    printf("%d\n",ans);

    return 0;
}