题意:给定一个n*m的01矩阵,求不严格次大子矩阵的大小。(子矩阵中元素必须都为1)
就是在求最大子矩阵的过程中加入一些判断,来求次大子矩阵。
分为三种情况:
- 如果该子矩阵的值小于最大值,那么用这个子矩阵的值去更新次大值
- 如果该子矩阵的值等于最大值,首先判断这个矩阵是否被使用过,若没有使用过则次大值等于最大值
- 如果该子矩阵的值大于最大值,那么次大值等于最大值,最大值等于该子矩阵的值。
而且次大值可能是最大值去掉一行一列得到的,所以也要考虑到这个。
标记一个矩阵可以使用矩阵的左下角左边以及矩阵的高度。
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e3+7;
int mapp[maxn][maxn];
char s1[maxn];
int l[maxn],r[maxn];
stack<int> s;
map<pair<pair<int,int>,int>,int> mmm;//标记出现过的矩阵;
int main(){
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i){
scanf("%s",s1+1);
for(int j=1;j<=m;++j){
mapp[i][j]=s1[j]-'0';
if(mapp[i-1][j]&&mapp[i][j]) mapp[i][j]+=mapp[i-1][j];
}
}
int maxx=0;
int mm=0;
for(int i=1;i<=n;++i){
while(!s.empty()) s.pop();
for(int j=1;j<=m;++j){
while(s.size()&&mapp[i][s.top()]>=mapp[i][j]) s.pop();
if(s.empty()) l[j]=0;
else l[j]=s.top();
s.push(j);
}
while(!s.empty()) s.pop();
for(int j=m;j>=1;--j){
while(s.size()&&mapp[i][s.top()]>=mapp[i][j]) s.pop();
if(s.empty()) r[j]=m+1;
else r[j]=s.top();
s.push(j);
if(maxx==mapp[i][j]*(r[j]-1-l[j])){
if(mmm.count(make_pair(make_pair(l[j],i),mapp[i][j]))){
}
else
mm=maxx;
}
else {
if(maxx>mapp[i][j]*(r[j]-1-l[j])) mm=max(mm,mapp[i][j]*(r[j]-1-l[j]));
else{
mm=maxx;
maxx=mapp[i][j]*(r[j]-1-l[j]);
}
mm=max(mm,mapp[i][j]*(r[j]-2-l[j]));
if(mapp[i][j])mm=max(mm,(mapp[i][j]-1)*(r[j]-1-l[j]));
}
mmm[make_pair(make_pair(l[j],i),mapp[i][j])]++;
}
}
printf("%d\n",mm);
return 0;
}