枚举上下界,上下界确定后在区间中统计答案。
对于这些上下界相同,左右端不同的矩形,可能的贡献对象在每一列的最小值中。
记一下每一列的最小值,然后找左边、右边比它小的第一个位置,就可以统计出该列的贡献了。用单调栈维护即可。
#include<bits/stdc++.h>
#define cs const
#define re register
cs int N=305;
namespace IO{
cs int Rlen=1<<22|1;
char buf[Rlen],*p1,*p2;
inline char gc(){return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;}
template<typename T>
inline T get(){
char ch=gc();T x=0;
while(!isdigit(ch)) ch=gc();
while(isdigit(ch)) x=((x+(x<<2))<<1)+(ch^48),ch=gc();
return x;
}
inline int gi(){return get<int>();}
}
using IO::gi;
int n,m,a[N][N],ans[N*N],mn[N];
int L[N],R[N],st[N],pos[N],top=0;
inline void Min(int &x,int y){if(x>y)x=y;}
inline void solve(int up,int dn){
for(int re i=1;i<=m;++i) Min(mn[i],a[dn][i]);
st[top=1]=0,pos[top]=0;
for(int re i=1;i<=m;++i){
while(top&&mn[i]<st[top]) --top;
L[i]=pos[top],st[++top]=mn[i],pos[top]=i;
}
st[top=1]=0,pos[top]=m+1;
for(int re i=m;i>=1;--i){
while(top&&mn[i]<st[top]) --top;
R[i]=pos[top],st[++top]=mn[i],pos[top]=i;
}
for(int re i=1;i<=m;++i)
ans[mn[i]]+=(i-L[i])*(R[i]-i);
}
int main(){
// freopen("submatrix.in","r",stdin);
n=gi(),m=gi();
for(int re i=1;i<=n;++i)
for(int re j=1;j<=m;++j)
a[i][j]=gi();
for(int re i=1;i<=n;++i){
memset(mn,0x3f,sizeof mn);
for(int re j=i;j<=n;++j) solve(i,j);
}
for(int re i=1;i<=n*m;++i)
printf("%d\n",ans[i]);
}