【BZOJ4569】萌萌哒

题目链接

BZOJ4569

解法：倍增+并查集

先想暴力的做法。每次读入$l_1,r_1,l_2,r_2$，令$i$从$0$枚举到$r_1-l_1$，然后用并查集合并$l_1+i$和$l_2+i$，最后统计并查集中联通块的个数$s$，答案即为$9\times 10^{s-1}$。
但是这样会TLE。
对暴力进行优化。更改并查集中数组$f$的含义，令$f_{x,k}$表示从$x$到$x+2^k-1$这两段的共同编号，亦即若$f_{x,k}=f_{y,k}$则表示$f_x,\dots ,f_{x+2^k-1}$和$f_y,\dots ,f_{y+2^k-1}$这两段数每项分别处于相同的联通块中。合并时用倍增的处理方式（参考ST表）。

代码

#include<iostream>
#include<cstdio>
#include<cmath>

using namespace std;

const long long p=1000000007ll;
long long ans=1ll;
int n,m,l1,l2,r1,r2,f[17][100001];
bool flag=false;

int find(int x,int k){
	return f[k][x]==x?x:f[k][x]=find(f[k][x],k);
}

void merge(int x,int y,int k){
	if(y>100000)return;
	int fx=find(x,k),fy=find(y,k);
	if(fx!=fy)f[k][fx]=fy;
}

void update(int a,int b,int c,int d){
	int k=(int)log2(b-a+1);
	merge(a,c,k),merge(b-(1<<k)+1,d-(1<<k)+1,k);
}

int main(){
	scanf("%d%d",&n,&m);int l=(int)log2(n);
	for(int i=1;i<=n;++i)for(int k=0;k<17;++k)f[k][i]=i;
	for(int i=1;i<=m;++i)scanf("%d%d%d%d",&l1,&r1,&l2,&r2),update(l1,r1,l2,r2);
	for(int i=l;i;--i)for(int j=1;j<=n-(1<<i)+1;++j){int f=find(j,i);merge(j,f,i-1),merge(j+(1<<i-1),f+(1<<i-1),i-1);}
	for(int i=1;i<=n;++i)if(find(i,0)==i)(ans*=(flag?10:((flag=true),9)))%=p;
	printf("%lld",ans);
}