poj 1050 最大子矩阵和

题意很明确了，给定矩阵求最大子矩阵使得和最大

基本方法：枚举，枚举x, y的上下边界需要O(n^4)时间复杂度，计算矩阵和需要O(n^2)的复杂度，总共要O(n^6)无法承受

改进： 

         采用二维前缀和（类似一维），sub[i][j]表示以i,j坐标为右下角的矩阵的和

         递推计算公式sub[i][j] = sub[x-1][y] + sub[x][y-1] - sub[x-1][y-1] + array[x][y] 需要O(n^2)的复杂度计算sub矩阵

        枚举：枚举纵坐标上下边界O(n^2)，固定纵坐标上下边界的情况下，从左到右扫描O(n)，一共O(n^3)

                    这里用到了一维子序列和的O(n)算法，见http://www.cnblogs.com/CCBB/archive/2009/04/25/1443455.html

                    因为在固定y1, y2的情况下，扫描x的时候可以在O(1)的时间里得到sigma(array[x][y1]...array[x][y2])的和。将二维

                    问题降为一维解决（sub数组发挥作用使得值能够在O(1)时间得到）

                    总时间复杂度为O(n^3)，相对于O（n^6)是一个改进

poj 1050

         

#include <iostream>
#include <vector>
#include <map>
#include <list>
#include <set>
#include <deque>
#include <stack>
#include <queue>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <cstdio>
#include <iomanip>
#include <cmath>
#include <cstdio>
#include <iostream>
#include <string>
#include <sstream>
#include <cstring>
#include <queue>
using namespace std;

///宏定义
const int  INF = 10000000;
const int MAXN = 110;
const int maxn = MAXN;
///全局变量 和 函数
///

int n;
int arrayy[maxn][maxn];
int sub[maxn][maxn];
int main()
{
	///变量定义
	int i, j;
    while(cin >> n)
    {
        memset(arrayy, 0, sizeof(arrayy));
        memset(sub, 0, sizeof(sub));
        for(i = 1; i <= n; i++)
        {
            for(j = 1;  j <= n; j++)
            {
                cin >> arrayy[i][j];
            }
        }

        //计算前缀和
        for(i = 1; i <= n; i++)
        {
            for(j = 1; j <= n; j++)
            {
                sub[i][j] = sub[i - 1][j] + sub[i][j - 1] - sub[i - 1][j - 1] + arrayy[i][j];
            }
        }

        //枚举最大值
        int ans = -INF;
        for(i = 1; i <= n; i++)
        {
            for(j = i; j <= n; j++) 
            {
               int maxSum = 0, thisSum = 0;
               for (int k = 1; k <= n; k++)//一维最大序列和联机算法O(n)

               {
                      thisSum += sub[j][k] + sub[i][k - 1] - sub[i][k] - sub[j][k - 1];  //本质上就是一维问题
                      if (thisSum > maxSum)
                             maxSum = thisSum;
                      else if (thisSum < 0)
                             thisSum = 0;
               }
               ans = max(ans, maxSum);
            }
        }
        cout << ans << endl;
    }



	///结束
	return 0;
}