求二维矩阵的最大子矩阵

1、使用积累数组cumarr[1...n][1...n]，使得求任意块儿子矩阵和的复杂度为O(1)，cumarr[i][j]为cumarr[1...i][1...j]子矩阵的和， 
    通过cumarr[i][j]=cumarr[i-1][j]+cumarr[i][j-1]-cumarr[i-1][j-1]计算积累数组，复杂度为O(N*N), 
2、枚举+dp 
    对矩阵的行枚举，确定矩阵区域的上下界imin和imax，此时就可以进行类似一维数组最大和的方法，每个枚举的矩阵区域负责度O(N)。 

 

总复杂度为O(N*N)+(N*N)*O(N)=O(N^3) 

 

 

   1: #include <stdio.h>
   2: #include <iostream>
   3:  
   4: using namespace std ;
   5:  
   6: //只用1...n，以便和积累数组对应
   7: int a[101][101];
   8: int cumarr[101][101];
   9:  
  10: void run1050()
  11: {
  12:     int n,val ;
  13:     int i,j ;
  14:     int imin,imax ;
  15:     int maxendinghere,maxsofar ;
  16:     int tmp ;
  17:     
  18:     scanf( "%d", &n );
  19:  
  20:     //积累数组边界值
  21:     for( i=0 ; i<=n ; ++i )
  22:     {
  23:         cumarr[0][i] = 0 ;
  24:         cumarr[i][0] = 0 ;
  25:     }
  26:  
  27:     //输入矩阵并求得积累数组
  28:     for( i=1 ; i<=n ; ++i )
  29:     {
  30:         for( j=1 ; j<=n ; ++j )
  31:         {
  32:             scanf( "%d", &(a[i][j]) ) ;
  33:             cumarr[i][j] = cumarr[i-1][j] + cumarr[i][j-1] - cumarr[i-1][j-1] + a[i][j] ;
  34:         }
  35:     }
  36:  
  37:     //对矩阵的行枚举，确定矩阵区域的上下界imin和imax，
  38:     maxsofar = a[1][1] ;
  39:     for( imin=1 ; imin<=n ; ++imin )
  40:     {
  41:         for( imax=imin ; imax<= n ; ++imax )
  42:         {
  43:             //进行类似一维数组最大和的方法
  44:             maxendinghere = cumarr[imax][1] - cumarr[imin-1][1] ;
  45:             for( j = 2 ; j<=n ; ++j )
  46:             {
  47:                 tmp = cumarr[imax][j] - cumarr[imin-1][j] - cumarr[imax][j-1] + cumarr[imin-1][j-1] ;
  48:                 maxendinghere = std::max( maxendinghere+tmp , tmp ) ;
  49:                 maxsofar = std::max( maxsofar , maxendinghere ) ;
  50:             }
  51:         }
  52:     }
  53:  
  54:     printf( "%d\n", maxsofar ) ;
  55: }