#4764 子矩阵
SOL
30pts:暴力维护即可(由于数据很强,暴力没法骗更多分)
100pts:写了可得30pts暴力的人应该知道可以考虑枚举行,固定一行
i
i
i去枚举后面每一行
j
j
j,然后可以用单调栈维护列的信息,复杂度
O
(
N
3
)
O(N^3)
O(N3)
代码:
#include<bits/stdc++.h>
#pragma GCC optimize(2)
using namespace std;
#define re register
inline int rd(){
int re data=0;static char ch=0;ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))data=(data<<1)+(data<<3)+(ch^48),ch=getchar();
return data;
}
const int N=301;
int n,m,a[N][N],top,s[N],t[N],r[N],l[N],mn[N],ans[N*N];
signed main(){
// freopen("submatrix.in","r",stdin);
// freopen("submatrix.out","w",stdout);
n=rd(),m=rd();
for(int re i=1;i<=n;++i)for(int re j=1;j<=m;++j)a[i][j]=rd();
for(int re i=1;i<=n;++i){
memset(mn,0x3f,sizeof(mn));
for(int re j=i;j<=n;++j){
for(int re k=1;k<=m;++k)mn[k]=min(mn[k],a[j][k]);
top=0;
for(int re k=1;k<=m;++k){
while(top&&mn[k]<s[top])r[t[top--]]=k-1;
s[++top]=mn[k],t[top]=k;
}while(top)r[t[top--]]=m;
for(int re k=m;k;--k){
while(top&&mn[k]<s[top])l[t[top--]]=k+1;
s[++top]=mn[k],t[top]=k;
}while(top)l[t[top--]]=1;
for(int re k=1;k<=m;++k)ans[mn[k]]+=(r[k]-k+1)*(k-l[k]+1);
}
}
for(int re i=1;i<=n*m;++i)printf("%d\n",ans[i]);
exit(0);
}