hdu 1081(最大子矩阵和)

　　题目很简单，就是个最大子矩阵和的裸题，看来算法课本的分析后也差不多会做了。利用最大子段和的O(n)算法，对矩阵的行（或列）进行 i和j的枚举，对于第 i到j行，把同一列的元素进行压缩，得到一整行的一维数组后直接调用O(n)算法即可。我一开始还想着同一列的元素压缩不是也要耗费O(n)的时间吗，看了书上的代码后才知道原来数组b[]的每个元素都可以利用上一次的结果在O(1)时间内算出（当 i固定，j向下枚举时），当 i移动时，b[]就要清零进行重新计算了（在这里很奇怪动态分配的数组竟然不能直接用memset来清零，必须手动开个for循环来清零的，为了先跳过这些细枝末节只好开个全局数组了），代码如下：

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 const int INF= 0x3fffffff;
 int a[][];

 int maxSum(int n, int c[]){
     int sum= -INF, b=;
     for(int i=; i<=n; ++i){
         if(b>=)    b+= c[i];
         else    b= c[i];
         if(b>sum)    sum=b;
     }
     return sum;
 }

 int b[];
 int maxMatrix(int n, int a[][]){
     int sum= -INF;
 //    int *b= new int[n+1];
 //    printf("%d\n",sizeof(b));
     for(int i=; i<=n; ++i){
         memset(b,,sizeof(b));
 //        for(int k=1; k<=n; ++k)        b[k]= 0;
         for(int j=i; j<=n; ++j){
             for(int k=; k<=n; ++k)      b[k]+= a[j][k];
             sum= max(sum, maxSum(n,b));
         }
     }
     return sum;
 }

 int main(){
     int n;
     while(~scanf("%d",&n)){
         for(int i=; i<=n; ++i)
             for(int j=; j<=n; ++j)
                 scanf("%d",a[i]+j);
         printf("%d\n",maxMatrix(n,a));
     }
 }