[FROM WOJ]#4764 子矩阵

#4764 子矩阵

传送门

SOL
30pts：暴力维护即可（由于数据很强，暴力没法骗更多分）
100pts：写了可得30pts暴力的人应该知道可以考虑枚举行，固定一行$i$去枚举后面每一行$j$，然后可以用单调栈维护列的信息，复杂度$O(N^3)$

代码：

#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);
}