题意
有一个01矩阵,q个操作,每次将一个子矩阵全变为1,问每次操作后1的联通块的个数
题解
并查集
将一个0变为1时,看他4个方向是否为1,然后进行合并,动态维护联通块的数量
问题是,如何快速找0呢
b
i
t
s
e
t
bitset
bitset 里的骚函数
s
t
.
_
F
i
n
d
_
n
e
x
t
(
x
)
st.\_Find\_next(x)
st._Find_next(x)
查找大于x开始的下一次1的位置
我们枚举行,在每一行使用哪个函数就可以了
复杂度
O
(
n
∗
m
∗
q
/
64
)
O(n*m*q/64)
O(n∗m∗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;
}