【BZOJ4537】[Hnoi2016]最小公倍数分块

【BZOJ4537】[Hnoi2016]最小公倍数

Description

　　给定一张N个顶点M条边的无向图(顶点编号为1,2,…,n)，每条边上带有权值。所有权值都可以分解成2^a*3^b的形式。现在有q个询问，每次询问给定四个参数u、v、a和b，请你求出是否存在一条顶点u到v之间的路径，使得路径依次经过的边上的权值的最小公倍数为2^a*3^b。注意：路径可以不是简单路径。下面是一些可能有用的定义：最小公倍数：K个数a₁,a₂,…,a_k的最小公倍数是能被每个a_i整除的最小正整数。路径：路径P:P₁,P₂,…,P_k是顶点序列，满足对于任意1<=i<k，节点P_i和P_i+1之间都有边相连。简单路径：如果路径P:P₁,P₂,…,P_k中，对于任意1<=s≠t<=k都有Ps≠Pt，那么称路径为简单路径。

Input

　　输入文件的第一行包含两个整数N和M，分别代表图的顶点数和边数。接下来M行，每行包含四个整数u、v、a、b代表一条顶点u和v之间、权值为2^a*3^b的边。接下来一行包含一个整数q，代表询问数。接下来q行，每行包含四个整数u、v、a和b，代表一次询问。询问内容请参见问题描述。1<=n,q<=50000、1<=m<=100000、0<=a,b<=10^9

Output

　　对于每次询问，如果存在满足条件的路径，则输出一行Yes，否则输出一行 No（注意：第一个字母大写，其余字母小写）。

Sample Input

4 5
1 2 1 3
1 3 1 2
1 4 2 1
2 4 3 2
3 4 2 2
5
1 4 3 3
4 2 2 3
1 3 2 2
2 3 2 2
1 3 4 4

Sample Output

Yes
Yes
Yes
No
No

题解：最小公倍数=2^a*3^b等价于存在一个连通块使得边权的2的次数最大值=a且3的次数最大值=b，先给出一种暴力方法：

将所有边和询问按a从小到大排序，然后对于每个询问，将a不超过它的所有边按b排序，然后从小到大扔到并查集中，并用并查集维护连通块中b的最大值，如果最后b的最大值=询问的最大值，则输出Yes。

如何优化这个过程呢？考虑按a分块，先将所有询问按(a值所在块-1)排序，然后我们处理第i块中的询问。先将前i块的所有边以及询问按b排序，然后一个一个处理询问和连边。但问题是，对当前询问产生贡献的某些边可能在第i+1块中，但是这些边不超过sqrt(n)个，于是用按秩合并的并查集暴力加入再暴力还原即可。（我写的是按siz合并）。

#include <iostream>
#include <cstring>
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=50010;
int n,m,Q,B,ltp,tp,tq,ltq,top;
struct edge
{
	int a,b,A,B,org;
}p[maxn<<1],q[maxn];
int ans[maxn],f[maxn],sz[maxn],sa[maxn],sb[maxn],px[maxn],psa[maxn],psb[maxn],psz[maxn],pf[maxn],vis[maxn];
bool cmpa(const edge &a,const edge &b)
{
	return (a.A==b.A)?(a.B<b.B):(a.A<b.A);
}
bool cmpb(const edge &a,const edge &b)
{
	return (a.B==b.B)?(a.A<b.A):(a.B<b.B);
}
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<‘0‘||gc>‘9‘)	{if(gc==‘-‘)f=-f;	gc=getchar();}
	while(gc>=‘0‘&&gc<=‘9‘)	ret=ret*10+gc-‘0‘,gc=getchar();
	return ret*f;
}
int find(int x)
{
	return (f[x]==x)?x:find(f[x]);
}
void link(int x)
{
	int a=find(p[x].a),b=find(p[x].b);
	if(a==b)
	{
		sa[a]=max(sa[a],p[x].A),sb[a]=max(sb[a],p[x].B);
		return ;
	}
	if(sz[a]<sz[b])	swap(a,b);
	f[b]=f[a],sz[a]+=sz[b],sa[a]=max(sa[a],max(sa[b],p[x].A)),sb[a]=max(sb[a],max(sb[b],p[x].B));
}
void add(int x)
{
	if(!vis[x])	vis[x]=1,px[++top]=x;
}
int main()
{
	n=rd(),m=rd(),B=int(sqrt(double(m)));
	int i,j,k,l,ra,rb;
	for(i=1;i<=m;i++)	p[i].a=rd(),p[i].b=rd(),p[i].A=rd(),p[i].B=rd();
	sort(p+1,p+m+1,cmpa);
	Q=rd();
	for(i=1;i<=Q;i++)	q[i].a=rd(),q[i].b=rd(),q[i].A=rd(),q[i].B=rd(),q[i].org=i;
	sort(q+1,q+Q+1,cmpa);
	p[0].A=p[0].B=-1,p[m+1].A=p[m+1].B=1<<30;
	for(i=0;i<=m;ltq=tq,i=tp+B)
	{
		tp=min(i,m);
		for(;tq<Q&&q[tq+1].A>=p[tp].A&&q[tq+1].A<p[min(tp+B,m+1)].A;tq++);
		sort(p+1,p+tp+1,cmpb),sort(q+ltq+1,q+tq+1,cmpb);
		for(j=1;j<=n;j++)	f[j]=j,sa[j]=sb[j]=-1,sz[j]=1;
		for(l=1,j=ltq+1;j<=tq;j++)
		{
			for(;l<=tp&&p[l].B<=q[j].B;link(l++));
			top=0;
			for(k=tp+1;p[k].A<=q[j].A;k++)	if(p[k].B<=q[j].B)
				ra=find(p[k].a),rb=find(p[k].b),add(ra),add(rb);
			for(k=1;k<=top;k++)	psa[k]=sa[px[k]],psb[k]=sb[px[k]],psz[k]=sz[px[k]];
			for(k=tp+1;p[k].A<=q[j].A;k++)	if(p[k].B<=q[j].B)	link(k);
			ra=find(q[j].a),rb=find(q[j].b);
			if(ra==rb&&sa[ra]==q[j].A&&sb[ra]==q[j].B)	ans[q[j].org]=1;
			for(k=1;k<=top;k++)	f[px[k]]=px[k],sa[px[k]]=psa[k],sb[px[k]]=psb[k],sz[px[k]]=psz[k],vis[px[k]]=0;
		}
	}
	for(i=1;i<=Q;i++)
	{
		if(ans[i])	printf("Yes\n");
		else	printf("No\n");
	}
	return 0;
}//