ACM-ICPC 2018 南京赛区网络预赛 B题 - 单调队列

题目链接:点击这里

 

解题思路:

此问题相当于在01矩阵中找全1的子矩阵个数.

从下到上，从左到右枚举点作为子矩阵的右上角.那么就是去找合法的左下角了.

假设现在以(x,y)作为右上角,在(x,y)的左下角去找一个点作为子矩阵的左下角使得这个子矩阵不会包含任意的0点

那么(x,y)左下角的0点坐标会围成一个高度递减的楼梯型围线,它与x轴,y轴组成一个多边形，那么在这多边形内的点作为左下角一个会包含0点,所以我们就是去找多边形外的点.由于是构成高度递减的围线所以我们就想到了用单调队列来维护这个问题.

 

#include<bits/stdc++.h>
#define fi first
#define se second
using namespace std;
typedef long long ll;
const int mx = 1e5 + 5;
vector <int> vec[mx];
int hig[120];
int n,m,q,top;
pair <int,int> sta[120];
int main(){
    int a,b,t,ca = 1;
    scanf("%d",&t);
    while(t--){
	    scanf("%d%d%d",&n,&m,&q);
	    for(int i=1;i<=n;i++) vec[i].clear();
	    memset(hig,0,sizeof(hig));
	    for(int i=0;i<q;i++){
	    	scanf("%d%d",&a,&b);
			vec[a].push_back(b); 
		}
		ll ans = 0;
		for(int i=1;i<=n;i++){
			for(int j=0;j<vec[i].size();j++) 
				hig[vec[i][j]] = i;
			int top = 0;
			ll ret = 0;
			for(int j=1;j<=m;j++){
				int y = i - hig[j];
				if(!y) top = 0,ret = 0;
				else{
					int x = 1;
					while(top&&sta[top].se>y){
						ret -= 1ll*sta[top].fi*sta[top].se;
						x += sta[top].fi;
						top--;
					}
					sta[++top] = make_pair(x,y);
					ret += 1ll*x*y;
					ans += ret; 
				}
			}
		}
		printf("Case #%d: %lld\n",ca++,ans);
	}
    return 0;
}