C++解决最大大于0子矩阵和

Description
_AP神牛准备给自己盖一座很华丽的宫殿。于是，他看中了一块N*M的矩形空地。空地中每个格子都有自己的海拔高度。AP想让他的宫殿的平均海拔在海平面之上（假设海平面的高度是0，平均数都会算吧？）。而且，AP希望他的宫殿尽量大，能够容纳更多的人来膜拜他。请问AP的宫殿最后会有多大？
Input Format
第一行为N和M。之后N行，每行M个数，描述的空地的海拔。
Output Format
输出一行,表示宫殿最大面积。
_
Sample Input
3 2
4 0
-10 8
-2 -2

Sample Output
4

Data Limit
对于30%的数据，N,M≤50;
对于100%的数据，N,M≤200；
数据很大，O（n^4）肯定过不了，希望大神能在不改变源程序的思路的基础上将此题改好，代码如下

 #include<iostream>
using namespace std;
int n,m,s[201][201],sum[201],f[201];
int main()
{
    int i,j,ans=0,j1,j2,k,p,t;
    cin>>n>>m; 
    for(i=1;i<=n;i++) 
    for(j=1;j<=m;j++) 
    { 
        cin>>t;
        s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+t; 
    }
    for(j1=1;j1<=n;j1++)
    for(j2=j1;j2<=n;j2++)
    {
        for(k=1;k<=n;k++)
        sum[k]=s[k][j2]-s[k][j1-1];
        for(k=1;k<=n;k++)
        if(sum[k]>sum[k-1])
        {
            f[k]=k;
            for(p=0;p<=f[k-1];p++)
            if(sum[k]-sum[p]>0)
            {
                f[k]=p;
                if((j2-j1+1)*(k-p)>ans)
                ans=(j2-j1+1)*(k-p);
            }
            else
            {
                f[k]=k;
                for(p=f[k-1];p<=k-1;p++)
                if(sum[k]-sum[p]>0)
                {
                    f[k]=p;
                    if((j2-j1+1)*(k-p)>ans)
                    ans=(j2-j1+1)*(k-p);
                }
            }
        }
    }
    cout<<ans;
    return 0;
}

http://blog.csdn.net/pzler/article/details/40116845

for(j1=1;j1<=n;j1++)
for(j2=j1;j2<=n;j2++)
{
for(k=1;k<=n;k++)
sum[k]=s[k][j2]-s[k][j1-1];

k应从1 to m而非 n
我原来写过牛宫 不知道还对不对
var
s,a:array[0..200,0..200]of longint;
sum,x:array[0..200]of longint;
i,j,n,m,i1,i2,z,p,ans:longint;
begin
readln(n,m);
for i:=1 to n do begin
for j:=1 to m do read(a[i,j]);
readln;
end;
for i:=1 to n do
for j:=1 to m do s[i,j]:=s[i-1,j]+s[i,j-1]-s[i-1,j-1]+a[i,j];
for i1:=1 to m do
for i2:=i1 to m do begin
for i:=1 to n do
sum[i]:=s[i,i2]-s[i,i1-1];
for i:=1 to n do
if sum[i]>sum[i-1] then begin
x[i]:=i;
for p:=0 to x[i-1] do
if sum[i]>sum[p] then begin
x[i]:=p;
z:=(i2-i1+1)*(i-p);
if z>ans then begin writeln(z);ans:=z;end;

break;
end;
end else begin
x[i]:=i;
for p:=x[i-1] to i-1 do
if sum[i]>sum[p] then begin
x[i]:=p;
z:=(i2-i1+1)*(i-p);
if z>ans then begin writeln(z);ans:=z;end;
break;
end;
end;
end;
writeln(ans);
end.

pascal的