一、 一维最大连续子序列问题
例题:KY141 最大连续子序列
描述
给定K个整数的序列{ N1, N2, ..., NK },其任意连续子序列可表示为{ Ni, Ni+1, ..., Nj },其中 1 <= i <= j <= K。最大连续子序列是所有连续子序列中元素和最大的一个,例如给定序列{ -2, 11, -4, 13, -5, -2 },其最大连续子序列为{ 11, -4, 13 },最大和为20。现在增加一个要求,即还需要输出该子序列的第一个和最后一个元素。
输入描述:
测试输入包含若干测试用例,每个测试用例占2行,第1行给出正整数K( K< 10000 ),第2行给出K个整数,中间用空格分隔。当K为0时,输入结束,该用例不被处理。
输出描述:
对每个测试用例,在1行里输出最大和、最大连续子序列的第一个和最后一个元素,中间用空格分隔。如果最大连续子序列不唯一,则输出序号i和j最小的那个(如输入样例的第2、3组)。若所有K个元素都是负数,则定义其最大和为0,输出整个序列的首尾元素。
示例1
输入:
6 -2 11 -4 13 -5 -2 10 -10 1 2 3 4 -5 -23 3 7 -21 6 5 -8 3 2 5 0 1 10 3 -1 -5 -2 3 -1 0 -2 0
输出:
20 11 13 10 1 4 10 3 5 10 10 10 0 -1 -2 0 0 0
AC 代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn = 10010;
int a[maxn];
void maxSum1D(int a[],int n)
{
int max_sum,cur_sum,maxst,maxed,st;
max_sum=a[0];
cur_sum=0;
st=maxst=maxed=0;
for(int i=0;i<n;i++)
{
if(cur_sum<0)
{
cur_sum=a[i];
st=i;
}
else
{
cur_sum+=a[i];
}
if(cur_sum>max_sum)
{
maxst=st;
maxed=i;
max_sum=cur_sum;
}
}
//若所有K个元素都是负数,则定义其最大和为0,输出整个序列的首尾元素。
if(max_sum<0)
{
cout<<0<<" "<<a[0]<<" "<<a[n-1]<<endl;
}
else
{
cout<<max_sum<<" "<<a[maxst]<<" "<<a[maxed]<<endl;
}
}
int main()
{
int k;
while(cin>>k)
{
if(k==0) break;
for(int i=0;i<k;i++)
{
cin>>a[i];
}
maxSum1D(a , k);
}
return 0;
}
二、 二维最大连续子序列问题
例题:KY63 最大子矩阵
描述
已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵,你的任务是找到最大的非空(大小至少是1 * 1)子矩阵。 比如,如下4 * 4的矩阵 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2 的最大子矩阵是 9 2 -4 1 -1 8 这个子矩阵的大小是15。
输入描述:
输入是一个N * N的矩阵。输入的第一行给出N (0 < N <= 100)。 再后面的若干行中,依次(首先从左到右给出第一行的N个整数,再从左到右给出第二行的N个整数……)给出矩阵中的N2个整数,整数之间由空白字符分隔(空格或者空行)。 已知矩阵中整数的范围都在[-127, 127]。
输出描述:
输出最大子矩阵的大小。
示例1
输入:
4 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2
输出:
15
AC 代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn=110;
int matrix[maxn][maxn];
int maxSum1D(int a[],int n)
{
int maxSum,curSum,st,maxSt,maxEd;
maxSum = INT_MIN;
curSum=0;
for(int i=0;i<n;i++)
{
if(curSum<0)
{
curSum=a[i];
st=i;
}
else
{
curSum+=a[i];
}
if(curSum>maxSum)
{
maxSum=curSum;
maxSt=st;
maxEd=i;
}
}
return maxSum;
}
int maxSum2D(int a[][maxn],int n,int m)
{
int maxSum = INT_MIN;
int *colSum = new int[m+1];
int curSum = 0;
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
colSum[j]=0;
}
for(int j=i;j<n;j++)
{
for(int k=0;k<m;k++)
{
colSum[k]+=a[j][k];
}
curSum = maxSum1D(colSum, m);
if(curSum>maxSum)
{
maxSum=curSum;
}
}
}
return maxSum;
}
int main()
{
int n,ans;
cin>>n;
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
cin>>matrix[i][j];
}
}
ans=maxSum2D(matrix,n,n);
cout<<ans<<endl;
return 0;
}
三、 三维最大连续子序列问题
例题:最大长方体问题
题目链接:最大长方体问题 | NOJ (njupt.edu.cn)https://acm.njupt.edu.cn/problem/NOJ1232/editor
最大长方体问题
Description:
一个长,宽,高分别为m,n,p的长方体被分割成个m*n*p个小立方体。每个小立方体内有一个整数。试设计一个算法,计算出所给长方体的最大子长方体。子长方体的大小由它所含所有整数之和确定。
对于给定的长,宽,高分别为m,n,p的长方体,编程计算最大子长方体的大小。
Input:
输入的第1行是3个正整数m,n,p,1≤m,n,p≤50。接下来m*n行每行p个正整数,表示小立方体中的数。
Output:
输出的第1行中的数是计算出的最大子长方体的大小。
Sample Input:
3 3 3 0 -1 2 1 2 2 1 1 -2 -2 -1 -1 -3 3 -2 -2 -3 1 -2 3 3 0 1 3 2 1 -3
Sample Output:
14
AC 代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn=55;
int matrix[maxn][maxn][maxn];
int maxSum1D(int a[],int n)
{
int maxSum,curSum,st,maxSt,maxEd;
maxSum = INT_MIN;
curSum = 0;
maxSt = maxEd = st = 0;
for(int i=0;i<n;i++)
{
if(curSum<0)
{
curSum=a[i];
st=i;
}
else
{
curSum+=a[i];
}
if(curSum>maxSum)
{
maxSum=curSum;
maxSt = st;
maxEd = i;
}
}
return maxSum;
}
int maxSum2D(int a[][maxn],int n,int p)
{
int maxSum,curSum;
int *colSum = new int[p+1];
maxSum=INT_MIN;
curSum=0;
for(int i=0;i<n;i++)
{
for(int j=0;j<p;j++)
{
colSum[j]=0;
}
for(int j=i;j<n;j++)
{
for(int k=0;k<p;k++)
{
colSum[k]+=a[j][k];
}
curSum = maxSum1D(colSum,p);
if(curSum>maxSum)
{
maxSum=curSum;
}
}
}
return maxSum;
}
int maxSum3D(int a[][maxn][maxn],int m,int n,int p)
{
int maxSum=INT_MIN;
int levelSum[maxn][maxn];
int curSum = 0;
for(int i=0;i<m;i++)
{
for(int j=0;j<n;j++)
{
for(int k=0;k<p;k++)
{
levelSum[j][k]=0;
}
}
for(int j=i;j<m;j++)
{
for(int k=0;k<n;k++)
{
for(int l=0;l<p;l++)
{
levelSum[k][l]+=a[j][k][l];
}
}
curSum = maxSum2D(levelSum,n,p);
if(curSum>maxSum)
{
maxSum=curSum;
}
}
}
return maxSum;
}
int main()
{
int m,n,p,ans;
cin>>m>>n>>p;
for(int i=0;i<m;i++)
{
for(int j=0;j<n;j++)
{
for(int k=0;k<p;k++)
{
cin>>matrix[i][j][k];
}
}
}
ans = maxSum3D(matrix,m,n,p);
cout<<ans<<endl;
return 0;
}