【XSY2741】网格 分治 LCT 并查集
题目描述
有一个\(n\times m\)的网格，线框的交点可以扭动，边不可伸缩。网格中有一些格子里面放了'x'形的支架，这些格子不会变形，但可以整体转动。如果所有格子都不能变形，那么称这个网格稳固。
有\(q\)个操作，每次改变一个格子的状态，即有支架给为无支架，无支架改为有支架。
请你判断初始及每次操作后这个网格是否稳固。
比如说下面这个网格就不稳固。
\(n,m\leq 3000,q\leq 100000\)
题解
先看看怎么判断一个网格是否稳固。
先给整个网格的左边和上边加上一行一列，然后把第一行的最左边两个格子设为有支架的格子，其他的设为没支架的格子。
样例那个图就会变成这样
可以发现这样操作是不会改变整个图形的稳定性的。
设格子\((i,j)\)右下角的角度为\(a_{i,j}+90\)。
因为一个交点四个角的度数和为\(360\)，所以可以列出以下方程：
然后通过一些简单变换可以得到
因为\(a_{0,0}=0\)，所以方程简化为
当\((i,j)\)有支架时\(a_{i,0}+a_{0,j}=a_{i,j}=0\)，即\(a_{i,0}=-a_{0,j}\)，那么我们就在图\(G\)的\(i\)和\(j+n\)两个点之间连一条边。
显然这个图是二分图。
因为边界上只有\(a_{0,1}=0\)，所以一个点只有和\(n+1\)号点\((0,1)\)属于同一个联通块，这个点对应的角的角度才是确定的。
当\((i,j)\)无支架时\(a_{i.j}=a_{i,0}+a_{0,j}\)。如果\(a_{i,0}\)和\(a_{0,j}\)之间有一个没有确定，那么\(a_{i,j}\)也没有确定。
所以说，这个网格是稳定的\(\Longleftrightarrow\)图\(G\)只有一个联通块。
现在问题就变成了：有一个\(n+m\)个点的图，有\(nm\)条边，还有\(q\)个加边删边的操作。问操作前和每一次操作完后联通块个数是不是\(1\)。
用分治+并查集和LCT都可以做。
可以把一定存在的边先用路径压缩的并查集处理完。
分治+并查集：\(O(nm\alpha+q\log^2(n+m))\)
LCT：\(O(nm\alpha+q\log (n+m))\)
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<utility>
#include<iostream>
#include<vector>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
void open(const char *s)
{
#ifndef ONLINE_JUDGE
char str[100];
sprintf(str,"%s.in",s);
freopen(str,"r",stdin);
sprintf(str,"%s.out",s);
freopen(str,"w",stdout);
#endif
}
vector<pii> a[400010];
int f[10010];
int r[10010];
int ans[100010];
int s1[10010];//x
int s2[10010];//f[x]
int s3[10010];//r[f[x]]
int top;
int n,m,q;
char s[10010];
int c[3010][3010];
int find(int x)
{
return f[x]==x?x:find(f[x]);
}
int find2(int x)
{
return f[x]==x?x:f[x]=find(f[x]);
}
int num=0;
int merge(int x,int y)
{
x=find(x);
y=find(y);
if(x==y)
return 0;
if(r[x]>r[y])
swap(x,y);
top++;
s1[top]=x;
s2[top]=y;
s3[top]=r[y];
if(r[x]==r[y])
r[y]++;
f[x]=y;
return 1;
}
int merge2(int x,int y)
{
x=find(x);
y=find(y);
if(x==y)
return 0;
if(r[x]>r[y])
swap(x,y);
f[x]=y;
if(r[x]==r[y])
r[y]++;
return 1;
}
void back()
{
f[s1[top]]=s1[top];
r[s2[top]]=s3[top];
top--;
}
void add(int p,int l,int r,int x,int y,int L,int R)
{
if(l<=L&&r>=R)
{
a[p].push_back(pii(x,y));
return;
}
int mid=(L+R)>>1;
if(l<=mid)
add(p<<1,l,r,x,y,L,mid);
if(r>mid)
add((p<<1)|1,l,r,x,y,mid+1,R);
}
void solve(int l,int r,int p)
{
int now=top;
for(auto v:a[p])
if(merge(v.first,v.second))
num++;
if(l==r)
ans[l]=(num==n+m-1);
else
{
int mid=(l+r)>>1;
solve(l,mid,p<<1);
solve(mid+1,r,(p<<1)|1);
}
while(top>now)
{
back();
num--;
}
}
int main()
{
open("grid");
scanf("%d%d%d",&n,&m,&q);
memset(c,-1,sizeof c);
for(int i=1;i<=n;i++)
{
scanf("%s",s+1);
for(int j=1;j<=m;j++)
if(s[j]=='x')
c[i][j]=0;
}
for(int i=1;i<=n+m;i++)
{
f[i]=i;
r[i]=1;
}
int x,y;
for(int i=1;i<=q;i++)
{
scanf("%d%d",&x,&y);
if(~c[x][y])
{
add(1,c[x][y],i-1,x,y+n,0,q);
c[x][y]=-1;
}
else
c[x][y]=i;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(~c[i][j])
{
if(c[i][j])
add(1,c[i][j],q,i,j+n,0,q);
else
if(merge2(i,j+n))
num++;
}
for(int i=1;i<=n+m;i++)
find2(i);
solve(0,q,1);
for(int i=0;i<=q;i++)
if(ans[i])
printf("S\n");
else
printf("U\n");
return 0;
}
