Comet OJ - Contest #13 「佛御石之钵 -不碎的意志-」(困难版)

题意

有一个01矩阵，q个操作，每次将一个子矩阵全变为1，问每次操作后1的联通块的个数

题解

并查集
将一个0变为1时，看他4个方向是否为1，然后进行合并，动态维护联通块的数量
问题是，如何快速找0呢
$bitset$ 里的骚函数 $st.\_Find\_next(x)$
查找大于x开始的下一次1的位置
我们枚举行，在每一行使用哪个函数就可以了
复杂度 $O(n\ast m\ast q/64)$ 约等于 5e8 的样子
玄学复杂度，我也不知道该除32还是64
反正 bitset 牛逼就对了

代码

#include<bits/stdc++.h>
#define N 1003
#define INF 0x3f3f3f3f
#define eps 1e-5
#define pi 3.141592653589793
#define mod 998244353
#define P 1000000007
#define LL long long
#define pb push_back
#define fi first
#define se second
#define cl clear
#define si size
#define lb lower_bound
#define ub upper_bound
#define bug(x) cerr<<#x<<"      :   "<<x<<endl
#define mem(x,y) memset(x,0,sizeof(int)*(y+3))
#define sc(x) scanf("%d",&x)
#define scc(x,y) scanf("%d%d",&x,&y)
#define sccc(x,y,z) scanf("%d%d%d",&x,&y,&z)
using namespace std;
typedef  pair<int,int> pp;

char s[N][N];
bitset<N> st[N];
bool f[N][N];
int n,m,op;
int fa[N*N];
inline int id(int a,int b){return (a-1)*m+b; }

pp ls[N*N];

int getfa(int x){
    return fa[x]==x?x:fa[x]=getfa(fa[x]);
}

int X[]={0,1,0,-1},Y[]={1,0,-1,0};

void dfs(int x,int y){
    int h=0,t=1;
    ls[1]={x,y};
    while(h!=t){
        h++; pp st=ls[h];
        for(int i=0;i<4;i++){
            int x=st.fi+X[i],y=st.se+Y[i];
            if (x<1||y<1||x>n||y>m||s[x][y]=='1'||f[x][y]==1) continue;
            f[x][y]=1;
            fa[id(x,y)]=op;
            ls[++t]={x,y};
        }
    }
}

int main(int argc, char const *argv[])
{
    scc(n,m);
    for(int i=1;i<=n;i++){
        scanf("%s",s[i]+1);
        for(int j=1;j<=m;j++) s[i][j]='1'-s[i][j]+'0',st[i][j]=s[i][j]-'0';
    }
    for(int i=1;i<=n*m;i++) fa[i]=i;
    int ans=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) if (s[i][j]=='0'&&f[i][j]==0){
            ans++;
            op=id(i,j);
            dfs(i,j);
        }
    for(int i=1;i<=n;i++) st[i][m+1]=1;
    int q; sc(q);
    while(q--){
        int x1,x2,y1,y2;
        scc(x1,y1); scc(x2,y2);
        for(int i=x1;i<=x2;i++){
            for(int j=st[i]._Find_next(y1-1);j<=y2;j=st[i]._Find_next(j)){
                s[i][j]='0'; st[i][j]=0;
                ans++; 
                for(int k=0;k<4;k++){
                    int x=i+X[k],y=j+Y[k];
                    if (x<1||y<1||x>n||y>m||s[x][y]=='1') continue;
                    int b=getfa(id(x,y)),a=getfa(id(i,j));
                    if (a!=b){
                        ans--;
                        fa[b]=a;
                    }
                }
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}