具体讲解可看我的B站视频蒟蒻关于最大子矩阵的和详细讲解,从一维过渡到三维!_哔哩哔哩_bilibili
需要提前知道的知识就是前缀和,大佬别喷我这个蒟蒻就可以了,我有鱼鱼症;
一维子矩阵的最大子矩阵和
例题cf中的Problem - B - Codeforces
#include <vector>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=100l;
int x;
int n; //n代表数据的数量
int sum=-1e5; //存过程,因为数据可能小于0
int res=-1e5; //存结果
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
sum=max(sum+x,x);
res=max(sum,res);
}
printf("%d\n",res);
return 0;
}
二维子矩阵的最大子矩阵和
例题洛谷P1719 最大加权矩形 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)
#include <vector>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=100l;
int arr[N][N];
int arr1[N][N]; //取前缀之后
int n,res=0,sum=0; //res存结果,sum存过程中的值
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
scanf("%d",&arr[i][j]);
arr1[i][j]=arr1[i][j-1]+arr[i][j];
}
for(int i=0;i<=n-1;i++)
for(int j=i+1;j<=n;j++)
{
sum=0;
for(int k=1;k<=n;k++)
{
sum=max(arr1[k][j]-arr1[k][i]+sum,arr1[k][j]-arr1[k][i]); //如果sum<0就从之后开始存
res=max(res,sum);
}
}
printf("%d\n",res);
return 0;
}
三维矩阵的最大子矩阵和
#include <vector>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N=26;
int arr[N][N][N];
int qz[N][N][N];
int z,x,y;
int sum; //存过程
int res; //存结果
void init()
{
for(int h=1;h<=z;h++)
for(int i=1;i<=x;i++)
for(int j=1;j<=y;j++)
scanf("%d",&arr[h][i][j]);
for(int h=1;h<=z;h++)
for(int i=1;i<=x;i++)
for(int j=1;j<=y;j++)
qz[h][i][j]=qz[h][i-1][j]+qz[h][i][j-1]-qz[h][i-1][j-1]+arr[h][i][j];
}
int main()
{
scanf("%d %d %d",&z,&x,&y);
init(); //初始化
for(int x1=1;x1<=x;x1++)
{
for(int y1=1;y1<=y;y1++)
{
for(int x2=x1;x2<=x;x2++)
{
for(int y2=y1;y2<=y;y2++)
{
sum=0;
for(int k=1;k<=z;k++)
{
int jie=qz[k][x2][y2]-qz[k][x1-1][y2]-qz[k][x2][y1-1]+qz[k][x1-1][y1-1];
sum =max(jie+sum,jie);
res=max(sum,res);
}
}
}
}
}
printf("%d\n",res);
return 0;
}