[OIBH] 糖果盒（Candy Box）——又一个最大子矩形

http://codewaysky.sinaapp.com/problem.php?id=1056

这题和奶牛浴场略有区别，奶牛浴场只需要求出最大子矩形，而这题要求的是最大权重子矩形，不一定要最大的面积，但要最大的权重和

思路是先求出每个最大子矩形，然后求出每个矩形的左上点和右下点，然后用二维数转数组进行求和，保存最优解

#include<iostream>
#include<string>
#include<algorithm>
#include<stdio.h>
#include<memory.h>
using namespace std;

int sum[1010][1010];
int v[1010][1010];
int l[1010],r[1010],h[1010];
int n,m;

//****************************//树状数组
int lowbit(int x)
{
    return x&-x;
}

void add(int x,int y,int w)
{
    int i,j;
    for(i=x;i<=n;i+=lowbit(i))
    {
        for(j=y;j<=m;j+=lowbit(j))
            sum[i][j]+=w;
    }
}

int get_sum(int x,int y)
{
    int i,j,ans=0;
    for(i=x;i>0;i-=lowbit(i))
    {
        for(j=y;j>0;j-=lowbit(j))
        {
            ans+=sum[i][j];
        }
    }
    return ans;
}
//*****************************//

int find(int x1,int y1,int x2,int y2)
{
    return get_sum(x2,y2)-get_sum(x1-1,y2)-get_sum(x2,y1-1)+get_sum(x1-1,y1-1);
}

int main()
{
    int i,j,w,x1,x2,y1,y2;
    freopen("D:\\in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)==2)
    {
        memset(sum,0,sizeof(sum));
        memset(v,0,sizeof(v));
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=m;j++)
            {
                scanf("%d",&w);
                if(!w)
                    v[i][j]=1;
                else
                    add(i,j,w);
            }
        }
        for(i=0;i<=m;i++)
        {
            h[i]=0;l[i]=1;r[i]=m;
        }
        int lm,rm,ans=0,temp;
        for(i=1;i<=n;i++)
        {
            lm=1;
            for(j=1;j<=m;j++)
            {
                if(!v[i][j])
                {
                    h[j]++;
                    if(lm>l[j])
                        l[j]=lm;
                }
                else
                {
                    h[j]=0; //边界不能有洞，所以障碍点高度是0，而不是1
                    l[j]=1;
                    r[j]=m;
                    lm=j+1; //由于边界不能有洞，所以加1
                }
            }
            rm=m;
            for(j=m;j>0;j--)
            {
                if(r[j]>rm)
                    r[j]=rm;
                if(h[j])
                {
                    y1=l[j];y2=r[j]; //处理出举行的左上点和右下点
                    x1=i-h[j]+1;x2=i;
                    temp=find(x1,y1,x2,y2);  //利用二维数转数组进行求和
                    if(temp>ans)
                        ans=temp;
                }
                else
                    rm=j-1; //同理减一
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}