CF869E The Untended Antiquity（二维数状数组+差分+hash）

考虑什么情况会不可达，当覆盖两点的最小矩形不同时，一定不可达。因此我们需要快速的知道覆盖一个点的最小矩形是哪个。我们考虑每次把一个矩形染色，那么复杂度是不可接受的。联想到我们一维做区间加法，单点查询时，可以使用差分+前缀和的方法解决。我们此时同样可以使用差分+前缀和的思想。只不过是二维的。给(x1,y1,x2,y2)的矩阵加x，我们可以这样做：(x1,y1)+x,(x1,y2+1)-x,(x2+1,y1)-x,(x2+1,y2+1)+x。有点容斥原理的感觉吧。然后用二维树状数组来维护这个前缀和。复杂度O(qlogn^2)。考虑到可能会有冲突的问题，我们需要hash解决。这里用了一个小质数来实现。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 2510
#define pa pair<int,int>
#define mod 10037
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,q;
ll c[N][N],num=1;
map<pair<pa,pa>,int>mp;
void add(int x,int y,int val){
    for(int i=x;i<=n;i+=i&(-i))
        for(int j=y;j<=m;j+=j&(-j)) c[i][j]+=val;
}
ll ask(int x,int y){
    ll res=0;
    for(int i=x;i;i-=i&(-i))
        for(int j=y;j;j-=j&(-j)) res+=c[i][j];return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();q=read();
    while(q--){
        int op=read(),x1=read(),y1=read(),x2=read(),y2=read();
        if(op==1){
            num*=mod;
            mp[make_pair(make_pair(x1,y1),make_pair(x2,y2))]=num;
            add(x1,y1,num);add(x1,y2+1,-num);add(x2+1,y1,-num);add(x2+1,y2+1,num);
        }if(op==2){
            int x=mp[make_pair(make_pair(x1,y1),make_pair(x2,y2))];
            add(x1,y1,-x);add(x1,y2+1,x);add(x2+1,y1,x);add(x2+1,y2+1,-x);
        }if(op==3){
            if(ask(x1,y1)==ask(x2,y2)) puts("Yes");
            else puts("No");
        }
    }return 0;
}