题目描述
一个长度为n的大数，用S1S2S3...Sn表示，其中Si表示数的第i位,S1是数的最高位，告诉你一些限制条件，每个条件表示为四个数，l1，r1，l2，r2，即两个长度相同的区间，表示子串Sl1Sl1+1Sl1+2...Sr1与Sl2Sl2+1Sl2+2...Sr2完全相同。
比如n=6时，某限制条件l1=1，r1=3，l2=4，r2=6，那么123123，351351均满足条件，但是12012，131141不满足条件，前者数的长度不为6，后者第二位与第五位不同。问满足以上所有条件的数有多少个。
输入输出格式
输入格式：
第一行两个数n和m，分别表示大数的长度，以及限制条件的个数。接下来m行，对于第i行,有4个数li1，ri1，li2，ri2，分别表示该限制条件对应的两个区间。1<=n<=10^5，1<=m<=10^5，1<=li1,ri1,li2,ri2<=n；并且保证ri1-li1=ri2-li2。
输出格式：
一个数，表示满足所有条件且长度为n的大数的个数，答案可能很大，因此输出答案模10^9+7的结果即可。
题意：
给出长度为n，并且有m个限制l1,r1,l2,r2，指区间l1-r1和l2-r2的数完全相同，求最后数的可能性；
题解：
①笨拙的我首先考虑暴力算法，枚举i = l1 to l2，然后把l1-r1和l2-r2的元素一一合并（并查集），但是太慢了。。。。。。。。
②太慢加速就好了，引入倍增并查集：
考虑f[x][y]为第x层，第一个位置为y，f[x][y]代表的x层区间y,y+(1<<x)-1的合并情况，每次我们把l1-r1区间长度二进制拆分后与l2-r2的对应拆分区间在正确的x上合并，统计答案的时候，如果f[x][y]和另一个标号为fa的区间合并了，那么就让x-1层的y和fa合并，y+(1<<x-1)和fa+(1<<x-1)合并；
#include<cstdio>
#include<iostream>
using namespace std;
const int N = ,mod = 1e9+;
int n,m,f[N][],p[];
char gc(){
static char *p1,*p2,s[];
if(p1==p2) p2=(p1=s)+fread(s,,,stdin);
return(p1==p2)?EOF:*p1++;
}
int rd(){
int x = ; char c = gc();
while(c<''||c>'') c = gc();
while(c>=''&&c<='') x=x*+c-'',c=gc();
return x;
}
int find(int j,int k){
return(f[j][k]==j)?j:f[j][k]=find(f[j][k],k);
}
void Union(int l1,int l2,int k){
if(find(l1,k) != find(l2,k))
f[f[l1][k]][k] = f[l2][k];
}
int main()
{ freopen("bzoj4569.in","r",stdin);
freopen("bzoj4569.out","w",stdout);
n=rd();m=rd();
for(int i= p[] = ;i <= ;i++) p[i] = p[i-]<<;
for(int i = ;i <= ;i++)
for(int j = ;j <= n;j++){
if(j+p[i]->n) break;
f[j][i] = j;
}
for(int i = ,l1,l2,r1,r2,len;i <= m;i++){
l1=rd(),r1=rd(),l2=rd(),r2=rd();
for(int k = ;k >= ;k--) if(l1+p[k]-<=r1) {
Union(l1,l2,k);
l1=l1+p[k],l2=l2+p[k];
}
}
for(int i = ;i >= ;i--)
for(int j = ;j+p[i]-<=n;j++){
if(find(j,i)!=j)
{
Union(j,f[j][i],i-);
Union(j+p[i-],f[j][i]+p[i-],i-);
}
}
int tot = ,ans=,tmp = ;
for(int i = ;i <= n;i++)
if(find(i,)==i) tot++;
tot--;while(tot){if(tot&) ans=1ll*ans*tmp%mod; tot>>=; tmp = 1ll*tmp*tmp%mod;}
printf("%d\n",ans);
return ;
}//by tkys_Austin;