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));
}
}
上一篇:GSON 简介 示例


下一篇:文件类型解析漏洞防御与攻击(PHP)