JZOJ 4820. 【NOIP2016提高A组模拟10.15】最大化

题目简介

给你一个矩阵，让你找到面积最大的一个，且这个矩阵的所有元素之和>0.

题解

首先，我们枚举上边界，左边界和右边界。
对于下边界，我们维护一个单调队列，然后按照$s_i-s_j(i,j为左边界和右边界)$从小到大二分。
如果$i>j$，那么说明这个状态可以更新答案。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#define N 305
#define LL long long
#define mo 1000000007
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
LL sum,i,j,k,l,r,mid,n,m,wz,ans,tot,a[N][N],s[N][N],d[N][2];
inline int read()
{
    int data=0,p=1;
    char ch=0;
    while (ch<'0' || ch>'9')
    {
        ch=getchar();
        if (ch=='-') p=-1; 
    }
    while (ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
    return data*p;
}
LL erfen(LL x)
{
    LL l=1,r=tot,mid;
    while (l<=r)
    {
        mid=(l+r)/2;
        if (d[mid][1]<x) r=mid-1;else l=mid+1;
    }
    return l;
}
int main()
{
    freopen("max.in","r",stdin);
    freopen("max.out","w",stdout);
    n=read();m=read();
    fo(i,1,n)
        fo(j,1,m)
        {
            a[i][j]=read();
            s[i][j]=s[i][j-1]+a[i][j];
        }
    fo(i,1,m)
        fo(j,i,m)
        {
            d[1][1]=2147483647;
            tot=1;sum=0;
            fo(k,1,n)
            {
                sum+=s[k][j]-s[k][i-1];
                wz=erfen(sum);
                if (wz<=tot) ans=max(ans,(k-d[wz][0])*(j-i+1));
                if (sum<d[tot][1])
                {
                    d[++tot][0]=k;
                    d[tot][1]=sum;
                }
            }
        }
    printf("%lld",ans);
}