HDU - 6992 Lawn of the Dead 线段树 + 思维

传送门

题意：

给你一张$n\ast m$的图，其中有$k$个点不能走，你只能向下和向右走，问你能到达多少点。
$n,m,k\le 1e5$

思路：

可以发现每个点如果其左边和上面都有障碍，那么这个点不可达。考虑到障碍数量比较少，所以肯定是枚举障碍的数量，现在就变成怎么快速统计答案了。
我们按照每一列来考虑，下面是样例的图，其中三角形代表障碍。
对于每个位置$(x,y)$，我们查询其它上面障碍出现的位置(如果之前没有障碍我们定义为$0$)之后到$(x-1,y-1)$的位置中第一次空位置，即能走的位置，将这个位置到当前位置之前都赋值为$1$，代表这些点可达。
下图中问号即为询问的区间，箭头指向$(x-1,y-1)$。

区间查询以及区间覆盖的操作显然可以用线段树来实现，再加点边界限制即可。

// Problem: H - Lawn of the Dead
// Contest: Virtual Judge - 2021多校第四场补题
// URL: https://vjudge.net/contest/450106#problem/H
// Memory Limit: 262 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;

int n,m,k;
vector<int>v[N];
struct Seg {
	struct Node{
		int l,r;
		int sum;
		int lazy;
	}tr[N<<2];
    
    void pushup(int u) {
    	tr[u].sum=tr[L].sum+tr[R].sum;
    }
	
	void pushdown(int u) {
		if(tr[u].lazy!=-1) {
			int lazy=tr[u].lazy; tr[u].lazy=-1; 
			tr[L].sum=Len(L)*lazy; tr[L].lazy=lazy;
			tr[R].sum=Len(R)*lazy; tr[R].lazy=lazy;
		} 
	}
	    
    void build(int u,int l,int r) {
    	tr[u]={l,r,0,-1};
    	if(l==r) return;
    	build(L,l,Mid); build(R,Mid+1,r);
    }
    
    void modify(int u,int l,int r,int x) {
    	if(tr[u].l>=l&&tr[u].r<=r) {
    		tr[u].sum=Len(u)*x;
    		tr[u].lazy=x;
    		return;
    	}
    	pushdown(u);
    	if(l<=Mid) modify(L,l,r,x);
    	if(r>Mid) modify(R,l,r,x);
    	pushup(u); 
    }
    
    int query(int u,int l,int r) {
    	if(tr[u].sum==0) return INF;
    	if(tr[u].l==tr[u].r) return tr[u].l;
    	pushdown(u);
    	if(tr[u].l>=l&&tr[u].r<=r) {
    		if(tr[L].sum>0) return query(L,l,r);
    		else return query(R,l,r);
    	}
    	int ans=INF;
    	if(l<=Mid) ans=min(ans,query(L,l,r));
    	if(r>Mid) ans=min(ans,query(R,l,r));
    	return ans;
    }
    
}t[2];

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);
	
	int _; scanf("%d",&_);
	for(int __=1;__<=_;__++) {
		scanf("%d%d%d",&n,&m,&k);
		for(int i=1;i<=m;i++) v[i].clear();
		for(int i=1;i<=k;i++) {
			int x,y; scanf("%d%d",&x,&y);
			v[y].pb(x);
		}	
		t[0].build(1,1,n); t[1].build(1,1,n);
		int now=0; t[now].modify(1,1,1,1); now^=1;
		LL ans=0;
		for(int i=1;i<=m;i++) {
			sort(v[i].begin(),v[i].end());
			int l=0;
			for(auto x:v[i]) {
				if(l+1<=x-1) {
					int pos=t[now^1].query(1,l+1,x-1);
					if(pos!=INF) t[now].modify(1,pos,x-1,1);
				}
				l=x;
			}
			if(l+1<=n) {
				int pos=t[now^1].query(1,l+1,n);
				if(pos!=INF) t[now].modify(1,pos,n,1);
			}
			ans+=t[now].tr[1].sum;
			t[now^1].modify(1,1,n,0);
			now^=1;
		}
		printf("%lld\n",ans);
		
	}



	return 0;
}
/*

*/