POJ 1050 to the max最大子矩阵

#include<iostream>
using namespace std;
int ko[105][105];
int b[105];
int dp[105];
int main()
{
int N,i,j,k,maxi,p,jo;
cin>>N;
for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
cin>>ko[i][j];
maxi=-1000000;
for(i=1;i<=N;i++)
{
memset(b,0,sizeof(b));
for(j=i;j<=N;j++)
{
for(k=1;k<=N;k++)
b[k]=b[k]+ko[j][k];
dp[1]=b[1];
jo=dp[1];
for(p=2;p<=N;p++)
{
if(dp[p-1]>0)
dp[p]=dp[p-1]+b[p];
else dp[p]=b[p];
if(dp[p]>jo)
jo=dp[p];
}
if(jo>maxi)
maxi=jo;
}
}
cout<<maxi<<endl;
return 0;
}