oldssoj1074糖果盒（悬线法）

题目描述

一个被分为 n*m　个格子的糖果盒，第 i 行第 j 列位置的格子里面有 a [ i ][ j ] 颗糖。 本来 tenshi 打算送这盒糖果给某 PPMM 的，但是就在要送出糖果盒的前一天晚上，一只极其可恶的老鼠夜袭糖果盒，有部分格子被洗劫并且穿了洞。 tenshi 必须尽快从这个糖果盒里面切割出一个矩形糖果盒，新的糖果盒不能有洞，并且 tenshi 希望保留在新糖果盒内的糖的总数尽量多。 计算一下新糖果盒最多能够保留多少糖果.

输入

第一行有两个整数 n、m。第 i + 1 行的第 j 个数表示 a [ i ][ j ]，如果这个数为 0 ，则表示这个位置的格子被洗劫过。 其中： 1 ≤ n，m ≤ 1000 0 ≤ a [ i ][ j ]≤ 255

输出

输出最大糖果数

样例输入

3 41 2 3 45 0 6 310 3 4 0

样例输出

17

提示

10 3 4
这个矩形的糖果数最大

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
using namespace std;
const int maxn=1003;
int n,m,H[maxn][maxn],minl[maxn][maxn],maxr[maxn][maxn],sum[maxn][maxn],a[maxn][maxn],h[maxn][maxn],l[maxn][maxn],r[maxn][maxn],ans=0;
inline int get(){
    char c;while(!isdigit(c=getchar()));
    int v=c-48;while(isdigit(c=getchar()))v=v*10+c-48;
    return v;
}
inline void geth(int x,int y){
    if(a[x][y]==0)h[x][y]=x;
    else h[x][y]=h[x-1][y];
}
inline void getl(int x,int y){
    if(a[x][y]==0)l[x][y]=y;
    else l[x][y]=l[x][y-1];
}
inline void getr(int x,int y){
    if(a[x][y]==0)r[x][y]=y;
    else r[x][y]=r[x][y+1];
}
int main(){
    memset(h,0,sizeof(h));
    memset(l,0,sizeof(l));
    memset(r,0,sizeof(r));
    memset(sum,0,sizeof(sum));
    n=get();m=get();
    for(int i=1;i<=n;++i){
        r[i][m+1]=m+1;
        for(int j=1;j<=m;++j){
            a[i][j]=get();
            sum[i][j]=sum[i-1][j]+sum[i][j-1]+a[i][j]-sum[i-1][j-1];
            geth(i,j);
            getl(i,j);
        }
    }
    for(int i=1;i<=n;++i)
        for(int j=m;j>=1;--j)getr(i,j);
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j){
            minl[i][j]=maxn;maxr[i][j]=0;           
            if(i==1){
                minl[i][j]=l[i][j];
                maxr[i][j]=r[i][j];
            }
            else{
                if(h[i-1][j]!=i-1){
                    minl[i][j]=max(minl[i-1][j],l[i][j]);
                    maxr[i][j]=min(maxr[i-1][j],r[i][j]);
                }
                else{                            //碰到上一个点是0的时候要特殊处理一下
                    minl[i][j]=l[i][j];
                    maxr[i][j]=r[i][j];
                }
            }
            int rr=maxr[i][j]-1,ll=minl[i][j],hh=h[i][j]; 
            ans=max(ans,sum[i][rr]-sum[i][ll]-sum[hh][rr]+sum[hh][ll]);
        }
    printf("%d\n",ans);
    return 0;
}

思路：悬线法。  这三条线即为悬线。每次在保证高度的情况下，向左向右寻找其最远边界（碰到0点或到边界），计算出该矩阵中的和。其中一点的最左、最右和高度可以预处理。矩阵的值也可以用前缀和预处理。正确性容易证明。