考虑将$(2i-1,2j-1)$和$(2i,2j)$缩为一个点,记作$(i,j)$

对于每一个点,只能选$(2i-1,2j-1)$或$(2i,2j)$(显然不能都选),而这样恰好为$nm$个,因此必须要至少选择一个

对于每一个点,障碍的状态分为以下几类:

1.无障碍,这类点暂时不考虑

2.都有障碍,若存在此类点必然不合法

3.有一个障碍,分两类讨论:

(1)$(2i-1,2j-1)$上有障碍,即这个点只能选$(2i,2j)$,考虑$(i+1,j)$和$(i,j+1)$,也只能选这个点,因此即对于其右下角的点,都要选$(2i,2j)$;

(2)$(2i,2j)$上有障碍,类似的,对于其左上角的点,都要选$(2i-1,2j-1)$

矛盾在于存在一个点$(i_{1},j_{1})$在$(2i_{1}-1,2j_{1}-1)$上有障碍,$(i_{2},j_{2})$在$(2i_{2},2j_{2})$上有障碍,且满足$i_{1}\le i_{2}$、$j_{1}\le j_{2}$,那么对于$([i_{1},i_{2}],[j_{1},j_{2}])$这个矩形就无解了

(其实第2类也可以看作此类情况,即$i_{1}=i_{2}$且$j_{1}=j_{2}$)

通过set,可以求出对于所有$i$,在$(2i-1,2j-1)$上有障碍的最小的$j$和在$(2i,2j)$上有障碍的最大的$j$

对于插入可以很好的支持,即判断新来的点是否会产生即可(通过线段树维护第一个前缀最小值以及第二个后缀最大值),那么删除通过线段树分治即可

 1 #include<bits/stdc++.h>
2 using namespace std;
3 #define N 200005
4 #define L (k<<1)
5 #define R (L+1)
6 #define mid (l+r>>1)
7 #define pii pair<int,int>
8 #define mp make_pair
9 #define fi first
10 #define se second
11 map<int,int>mat[N<<1];
12 map<int,int>::iterator it;
13 multiset<int>s1[N],s2[N];
14 vector<pii>v[N<<2];
15 pii f[N<<2];
16 int n,m,q,x,y,ans[N];
17 void add(int k,int l,int r,int x,int y,pii z){
18 if ((l>y)||(x>r))return;
19 if ((x<=l)&&(r<=y)){
20 v[k].push_back(z);
21 return;
22 }
23 add(L,l,mid,x,y,z);
24 add(R,mid+1,r,x,y,z);
25 }
26 pii merge(pii x,pii y){
27 return mp(min(x.fi,y.fi),max(x.se,y.se));
28 }
29 pii get(int k){
30 int x=m+1,y=0;
31 if (s1[k].size())x=(*s1[k].begin());
32 if (s2[k].size())y=(*--s2[k].end());
33 return mp(x,y);
34 }
35 void update(int k,int l,int r,int x){
36 if (l==r){
37 f[k]=get(x);
38 return;
39 }
40 if (x<=mid)update(L,l,mid,x);
41 else update(R,mid+1,r,x);
42 f[k]=merge(f[L],f[R]);
43 }
44 pii query(int k,int l,int r,int x,int y){
45 if ((l>y)||(x>r))return mp(m+1,0);
46 if ((x<=l)&&(r<=y))return f[k];
47 return merge(query(L,l,mid,x,y),query(R,mid+1,r,x,y));
48 }
49 void dfs(int k,int l,int r){
50 if (!ans[0])return;
51 for(int i=0;i<v[k].size();i++){
52 x=v[k][i].fi,y=v[k][i].se;
53 if (y&1){
54 s1[(x+1)/2].insert((y+1)/2);
55 if (query(1,1,n,(x+1)/2,n).se>=(*s1[(x+1)/2].begin()))ans[0]=0;
56 }
57 else{
58 s2[x/2].insert(y/2);
59 if (query(1,1,n,1,x/2).fi<=(*--s2[x/2].end()))ans[0]=0;
60 }
61 update(1,1,n,(x+1)/2);
62 }
63 if (l==r)ans[l]=ans[0];
64 else{
65 dfs(L,l,mid);
66 dfs(R,mid+1,r);
67 }
68 for(int i=0;i<v[k].size();i++){
69 x=v[k][i].fi,y=v[k][i].se;
70 if (y&1)s1[(x+1)/2].erase(s1[(x+1)/2].find((y+1)/2));
71 else s2[x/2].erase(s2[x/2].find(y/2));
72 update(1,1,n,(x+1)/2);
73 }
74 ans[0]=1;
75 }
76 int main(){
77 scanf("%d%d%d",&n,&m,&q);
78 for(int i=1;i<=q;i++){
79 scanf("%d%d",&x,&y);
80 if (!mat[x][y])mat[x][y]=i;
81 else{
82 add(1,1,q,mat[x][y],i-1,mp(x,y));
83 mat[x][y]=0;
84 }
85 }
86 for(int i=1;i<=2*n;i++)
87 for(it=mat[i].begin();it!=mat[i].end();it++)
88 if ((*it).se)add(1,1,q,(*it).se,q,mp(i,(*it).fi));
89 for(int i=1;i<=n;i++)update(1,1,n,i);
90 ans[0]=1;
91 dfs(1,1,q);
92 for(int i=1;i<=q;i++)
93 if (ans[i])printf("YES\n");
94 else printf("NO\n");
95 }

最新文章

  1. java反射详解
  2. 教你写能被舒服舒服又舒服地调用的iOS库
  3. polling 和 long polling 工作原理
  4. 恶心的Oracle的if else if...
  5. eclipse引用头文件报错问题-解决方法
  6. Java注解的简单了解
  7. jquery动态连接节点
  8. C# 委托和方法
  9. strut2服务器与android交互数据
  10. IS-IS完整笔记
  11. Fantacy团队周二站立会议
  12. 在Windows系统配置Jekyll
  13. [LeetCode&amp;Python] Problem 217. Contains Duplicate
  14. MapReduce文件配置和测试
  15. bzoj 3531 [Sdoi2014]旅行 (树剖+线段树 动态开点)
  16. vanilla
  17. Advanced Debugging and the Address Sanitizer
  18. python之旅:函数基础
  19. PHP isset和empty 详细比较
  20. JBOSS集群和安装

热门文章

  1. C语言日记② 初识 C的内存
  2. 测试开发【提测平台】分享12-掌握日期组件&amp;列表状态格式化最终实现提测管理多条件搜索展示功能
  3. 使用CEF(一)— 起步
  4. 永久修改alias
  5. Kali安装OWASP
  6. 初学Python “登录”案例 更新!!
  7. SharkCTF2021 fastcalc题记
  8. cs224n 2019
  9. Java:锁笔记
  10. MySQL:提高笔记-2