BZOJ4569 [Scoi2016]萌萌哒

一眼觉得是并查集，yy了一下线段树优化发现不可sing，orz了ljss神犇被告知是倍增优化

倍增，一个点[i][j]代表某个i往后2^j个数，把点按倍增次数分层，两个同层的点在一个连通块说明两个点所代表的区间相等

把每个操作拆到logn个点上，两两用并查集合并

然后从高层向低层下传，设当前层为j，假设x1,x2,x3……xt在一个连通块中，则将x2到xt的左儿子(即第j-1层的x)都与x1的左儿子合并，x2到xt的右儿子(即j-1层的x+2^(j-1))都与x1的右儿子合并，传到最底层的时候统计最底层有多少个连通块，设有tot个，答案即为9*10^(tot-1)

拆是m log n a(n)的，下传是n log n a(n) 的，总复杂度O((n+m)log n a(n))

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 100010
#define MAXM 20
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
int f[MAXN*20];
int n,m;
ll ans=1;
bool vis[MAXN*20];
int mi[MAXM];
int fa(int x){
	return f[x]==x?x:f[x]=fa(f[x]);
}
int pt(int x,int y){
	return (x-1)*20+y+1;
}
int fa(int x,int y){
	return fa(pt(x,y));
}
void uni(int x,int y,int l){
	int i;
	for(i=19;~i;i--){
		if(mi[i]<=l){
			f[fa(x,i)]=fa(y,i); 
			x+=mi[i];
			y+=mi[i];
			l-=mi[i];
		}
	}
}
int main(){
	//freopen("moe1.in","r",stdin);
	//freopen("moe.out","w",stdout);
	int i,j,k;
	int x,y,xx,yy;
	scanf("%d%d",&n,&m);
	mi[0]=1;
	for(i=1;i<20;i++){
		mi[i]=mi[i-1]<<1;
	}
	for(i=1;i<=n;i++){
		for(j=0;j<=19;j++){
			f[pt(i,j)]=pt(i,j);
		}
	}
	for(i=1;i<=m;i++){
		scanf("%d%d%d%d",&x,&y,&xx,&yy);
		uni(x,xx,yy-xx+1);
	}
	for(j=19;j;j--){
		for(i=1;i<=n;i++){
			if(i+mi[j]-1<=n&&fa(i,j)!=pt(i,j)){
				x=fa(i,j);
				x=(x-1)/20+1;
				f[fa(i,j-1)]=fa(x,j-1);
				f[fa(i+mi[j-1],j-1)]=fa(x+mi[j-1],j-1);
			}
		}
	}
	int tot=0;
	for(i=1;i<=n;i++){
		if(!vis[fa(i,0)]){
			//cout<<i<<' '<<j<<fa(i,0)<<endl;
			tot++;
			vis[fa(i,0)]=1;
		}
	}
	ans=9;
	for(i=1;i<tot;i++){
		(ans*=10)%=MOD;
	}
	printf("%lld\n",ans);
	return 0;
}

/*
8 3
1 2 3 4
5 6 7 8
3 6 5 8

4 2
1 2 3 4
3 3 4 4

2 1
1 1 2 2

*/