文章目錄

• • 普通矩陣（單個(gè)矩陣值為任何數(shù)）
  • • 最大子段和
    • 擴(kuò)展到二維情況
  • 01矩陣（單個(gè)矩陣值為0或1）
  • 代碼：

普通矩陣（單個(gè)矩陣值為任何數(shù)）

例題：POJ 1074

求出其中最大的子矩陣
答案是：
9 2
-4 1
-1 8
最大和是15
我們先想想如果不是矩陣，是一個(gè)數(shù)組，求其中連續(xù)的最長一段，咋做？

最大子段和

我們用b[i]來表示a[0]…a[1]的最大子段和
那么b [ i ] =max ( b [i - 1 ] + a [ i ], a [ i ] )
當(dāng)a[n] 大于dp[n-1]（前n-1個(gè)元素序列的最大子段和）。此時(shí)舍棄前面的子段和，即斷開。
否則的話，把a(bǔ)[n]連接到前n-1個(gè)元素序列的最大和子段上，即dp[n] = dp[n-1] + a[n]

int MaxSum(int n,int *a) {int sum=0,b=0;for(int i=0;i<n;i++){if(b>0) sum=max(sum,b+=a[i]);//不要忘了把b更新else sum=max(b=a[i],sum);}return sum; }

擴(kuò)展到二維情況

二維矩陣排列方式如圖
參考題解

如果紅色是我們要找的最大子矩陣
那和就是：（a[q][i]+…+a[p][i],a[q][i+1]+…+a[p][i+1],…,a[q][j]+…a[p][j]）=（sum₁，sum₂，…,sum_n）
sum也就是每一列的相同行數(shù)相加
這不也是sun的一維數(shù)組，那其實(shí)就是一個(gè)一維數(shù)組的最大子段問題
如果把二維數(shù)組看成是縱向的一維數(shù)組和橫向的一維數(shù)組

#include <iostream> #include <cstring> using namespace std;int maxsub(int a[], int n) {int i, max = a[0], b = 0;for (i = 0; i < n; i++){if (b + a[i] >= a[i]) //當(dāng)n-1的最大字段和+a[i]>=原來的最大字段和b += a[i]; //則n的最大字段和為b+a[i]else //否則，就是a[i]b = a[i];if (b > max) //找到最大的最大字段和max = b;}return max; }int main() {int n, i, j, k, maxsubrec, maxsubarr, m;int dp[101][101], arr[101];cin >> n >> m;for (i = 0; i < n; i++)for (j = 0; j < m; j++)cin >> dp[i][j];maxsubrec = dp[0][0];for (i = 0; i < n; i++) //i用于標(biāo)識(shí)從哪行開始，依次選擇行{memset(arr, 0, sizeof(arr));for (j = i; j < n; j++) //從第i行開始，j為選擇幾行{for (k = 0; k < m; k++) //從1列選擇到m列，壓縮行arr[k] += dp[j][k];//arr表示每列的值maxsubarr = maxsub(arr, m); //從每一個(gè)壓縮行中選擇最大字段和if (maxsubarr > maxsubrec) maxsubrec = maxsubarr;}}cout << maxsubrec << endl;return 0; }

01矩陣（單個(gè)矩陣值為0或1）

參考題解
POJ - 3494

其實(shí)就是沒有負(fù)數(shù)情況
我一開始想的是將最大子矩陣的代碼中，如果一個(gè)數(shù)為0就賦值為-10000（這樣程序就肯定不會(huì)選擇這個(gè)塊）
但是這樣會(huì)超時(shí)
針對01矩陣我們可以這么想:
每一個(gè)單元格的值等于它所在列的連續(xù)1的數(shù)量，如果當(dāng)前行為1，則當(dāng)前行的值就等于上一行加1，如果當(dāng)前行為0，則他的值也是0
我們掃描每一行，求出這一行中可形成的矩形最大面積，然后去最大情況即可
然后我們要用到一個(gè)單調(diào)遞減的單調(diào)棧，
如果棧為空或者當(dāng)前元素大于等于棧頂元素，則入棧
若棧非空且當(dāng)前元素小于棧頂元素，則將棧頂彈出，更新面積最大值，直到棧空或者遇到第一個(gè)大于等于當(dāng)前元素
將最后一個(gè)出棧的棧頂元素（最后一個(gè)大于當(dāng)前元素的）向左右延伸，修改其值，并入棧。
核心都是對于每個(gè)高度為h的矩形，找到最左高度大于等于它的位置，和最右的位置。

我們看第三行：
0 1 3 2 0 0 0
0先入棧，然后1入棧，然后3入棧，然后到2時(shí)，2<3,計(jì)算此時(shí)的矩陣面積=（4-3） * 3=3，此時(shí)最大面積是3，然后將值3改為2，此時(shí)就是0 1 2 2 ，然后0<2,計(jì)算此時(shí)面積 = (5-3) * 2=4,最大面積就是4,一次類推，都是按照這個(gè)操作

代碼：

#include<stdio.h> #include<string.h> #include<iostream> #include<stack> using namespace std;int main() {//top指向棧頂；tmp為臨時(shí)變量，記錄面積的值；ans為答案，記錄面積的最大值 int i,j,m,n,x,top,tmp,ans,h[2020],a[2020];stack<int> st; //單調(diào)棧，記錄位置 while(~scanf("%d%d",&m,&n)){ans=0;memset(h,0,sizeof(h)); //用于第一行的處理for(i=0;i<m;i++){ //掃描每一行 for(j=1;j<=n;j++){scanf("%d",&x);if(x==1) h[j]=h[j]+1; //如果是1，則在上一行本列的基礎(chǔ)上+1 else h[j]=0; //否則為0 a[j]=h[j]; //a數(shù)組用來向左右擴(kuò)展}a[n+1]=-1; //設(shè)最后元素為最小值，以最后讓棧內(nèi)元素出棧 for(j=1;j<=n+1;j++){if(st.empty()||a[j]>=a[st.top()]){ //如果棧為空或入棧元素大于等于棧頂元素，則入棧 st.push(j);}else{while(!st.empty()&&a[j]<a[st.top()]){ //如果棧非空并且入棧元素小于棧頂元素，則將棧頂元素出棧 top=st.top();st.pop();tmp=(j-top)*a[top]; //計(jì)算面積值 if(tmp>ans) ans=tmp; //更新面積最大值 }st.push(top); //將最后一次出棧的棧頂元素延伸并入棧 a[top]=a[j]; //修改其對應(yīng)的值 }}}printf("%d\n",ans);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的最大子矩阵（普通和01）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。