BZOJ1018 [SHOI2008]堵塞的交通traffic（分块+并查集）

这份题解没有在BZOJ上AC（狂汗不止），不过我觉得分块的思想挺不错的，就尝试了一下

不过相似的代码可以AC这道题LA 4082 Traffic Jam

如图，先给点编号，然后每一块中暴力维护并查集：



Open：在同一块中就维护并查集，在不同块互相连通就用一个link数组记录

Close：在同一块中就把块中元素的fa[]全部重新计算，在不同块中直接把link赋零

Ask：判断A及能到达的4个顶点，再递归4个顶点到其他块的顶点，看顶点与B是否连通

如果有哪位神犇能指点一下分块法如何能AC此题，或者想讨

论，欢迎随时与我交流^ ^

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define SIZE 512
int G[200005][5]={0},p[200005]={0},fa[200005]={0},jiao[200005/SIZE+5][5],link[200005]={0},vis[200005]={0};
int min(int a,int b)
{
	if(a<b) return a;
	return b;
}
int father(int x)
{
	if(fa[x]!=x) fa[x]=father(fa[x]);
	return fa[x];
}
int block(int x)
{
	return (x-1)/SIZE+1;
}
int cx(int a,int b)
{
	int i,x,BL=block(a);
	vis[a]==1;
	if(father(a)==father(b)) return 1;
	for(i=1;i<=4;i++)
		if(father(jiao[BL][i])==father(a))
		{
			vis[jiao[BL][i]]=1;
			x=link[jiao[BL][i]];
			if(x!=0&&vis[x]==0)
			if(cx(x,b)==1) return 1;
		}
	return 0;
}
int main()
{
	char s[20]={0};
	int n,i,j,r,c,a,b,tou,wei,BL;
	scanf("%d",&n);
	n*=2;
	for(i=1;i<=n;i++)
		fa[i]=i;
	for(i=1;i<=(n-1)/SIZE+1;i++)
	{
		jiao[i][1]=(i-1)*SIZE+1;
		jiao[i][2]=jiao[i][1]+1;
		jiao[i][4]=min(i*SIZE,n);
		jiao[i][3]=jiao[i][4]-1;
	}
	while(scanf("%s",s))
	{
		if(s[0]=='E') break;
		scanf("%d%d",&r,&c);
		a=2*(c-1)+r;
		scanf("%d%d",&r,&c);
		b=2*(c-1)+r;
		if(s[0]=='C')//关 
		{
			for(i=1;i<=p[a];i++)
				if(G[a][i]==b)
				{
					G[a][i]=G[a][p[a]--];
					break;
				}
			for(i=1;i<=p[b];i++)
				if(G[b][i]==a)
				{
					G[b][i]=G[b][p[b]--];
					break;
				}
			if(block(a)==block(b))
			{
				BL=block(a);
				tou=(BL-1)*SIZE+1;
				wei=min(BL*SIZE,n);
				for(i=tou;i<=wei;i++)
					fa[i]=i;
				for(i=tou;i<=wei;i++)
					for(j=1;j<=p[i];j++)
						if(block(G[i][j])==BL) fa[father(i)]=father(G[i][j]);
			}
			else link[a]=link[b]=0;
		}
		if(s[0]=='O')//开 
		{
			G[a][++p[a]]=b;
			G[b][++p[b]]=a;
			if(block(a)==block(b)) fa[father(a)]=father(b);
			else
			{
				link[a]=b;
				link[b]=a;
			}
		}
		if(s[0]=='A')//询问 
		{
			memset(vis,0,sizeof(vis));
			if(cx(a,b)==1) printf("Y\n");
			else printf("N\n");
		}
	}
	return 0;
}