poj 2029 get many persimmon trees

题目大意：

有一个01矩阵，求其所有a*b的子矩阵中总和最大的子矩阵的权值和

思路：

二维树状数组

n²log²n

虽然poj数据太弱了，二维前缀和就过了

好多天以后发现二位前缀和比树状数组快。。。

经过实验：一样快。。。

附两种代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<vector>
#include<set>
#include<stack>
#define inf 2147483611
#define ll long long
#define MAXN 510
using namespace std;
inline int read()
{
    int x=0,f=1;
    char ch;ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int map[MAXN][MAXN],x,y,s,t,n,ans;
int main()
{
    int a,b;
    while(scanf("%d",&n)&&n)
    {
        memset(map,0,sizeof(map));
        x=read(),y=read();
        for(int i=1;i<=n;i++)
        {
            a=read(),b=read();
            map[b][a]=1;
        }
        s=read(),t=read();
        ans=-1;
        for(int i=1;i<=y;i++)
            for(int j=1;j<=x;j++)
            {
                map[i][j]+=map[i-1][j]+map[i][j-1]-map[i-1][j-1];
            } 
        for(int i=t;i<=y;i++)
            for(int j=s;j<=x;j++)
            {
                ans=max(ans,map[i][j]-map[i-t][j]-map[i][j-s]+map[i-t][j-s]);
            } 
        printf("%d
",ans);
    }
}

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<vector>
#include<set>
#include<stack>
#define inf 2147483611
#define ll long long
#define MAXN 510
using namespace std;
inline int read()
{
    int x=0,f=1;
    char ch;ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int c[MAXN][MAXN],x,y,s,t,n,ans;
int lowbit(int x) {return x&(-x);}
void update(int a,int b)
{
    for(int i=a;i<=y;i+=lowbit(i))
        for(int j=b;j<=x;j+=lowbit(j))
        {
            c[i][j]++;
        }
}
int sum(int a,int b)
{
    int ans=0;
    for(int i=a;i>=1;i-=lowbit(i))
        for(int j=b;j>=1;j-=lowbit(j))
        {
            ans+=c[i][j];
        }
    return ans;
}
int main()
{
    int a,b;
    while(scanf("%d",&n)&&n)
    {
        memset(c,0,sizeof(c));
        x=read(),y=read();
        for(int i=1;i<=n;i++)
        {
            a=read(),b=read();
            update(b,a);
        }
        s=read(),t=read();
        ans=-1;
        for(int i=t;i<=y;i++)
            for(int j=s;j<=x;j++)
            {
                ans=max(ans,sum(i,j)-sum(i-t,j)-sum(i,j-s)+sum(i-t,j-s));
            } 
        printf("%d
",ans);
    }
}