woj4764 子矩阵

最新推荐文章于 2021-08-13 01:52:32 发布

「已注销」

最新推荐文章于 2021-08-13 01:52:32 发布

阅读量163

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/g21wcr/article/details/102700989

无传送门

枚举上下界，上下界确定后在区间中统计答案。
对于这些上下界相同，左右端不同的矩形，可能的贡献对象在每一列的最小值中。
记一下每一列的最小值，然后找左边、右边比它小的第一个位置，就可以统计出该列的贡献了。用单调栈维护即可。

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