没想到这两个东西还可以这么用,我们可以用
f[i][j]
表示
[i,i+2j−1]
这一段区间属于哪一个集合,没有则等于0.
那么每次合并我们可以把这个限制拆成log个区间,依次合并起来。
若
f[i][j]
和
f[s][t]
同属一个集合,那么
f[i][j−1]
和
f[s][t−1]
,
f[i+2j−1][j−1]
和
f[s+2t−1][t−1]
也同属一个集合。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 100005;
const long long mod = 1000000007;
template<class T>inline void read(T &res){
static char ch;T flag=1;
while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;
while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;
}
int tot,n,m,lx,ly,rx,ry,xt,yt,fa[N*17],Log[N*17],a[N*17],b[17],pos[N][17];
int find(int x){
return fa[x]?fa[x]=find(fa[x]):x;
}
int unionn(int x,int y){
if(find(x)!=find(y))fa[find(y)]=find(x);
}
int main(){
read(n),read(m),b[0]=1;
for(register int i=1;i<=16;++i)b[i]=b[i-1]*2;
for(register int i=1;i<=n;i++)
for(register int j=0;j<=16;j++)
pos[i][j]=++tot,a[tot]=i,Log[tot]=j;
for(register int i=1;i<=m;i++){
read(lx),read(ly),read(rx),read(ry);
for(register int j=16;j>=0;--j)
if(lx+b[j]-1<=ly){
unionn(pos[lx][j],pos[rx][j]);
rx+=b[j],lx+=b[j];
}
}
for(register int j=16;j>=1;--j)
for(register int i=1;i<=n;i++){
xt=a[find(pos[i][j])],yt=Log[find(pos[i][j])];
unionn(pos[xt][yt-1],pos[i][j-1]);
if(i+b[j-1]<=n)unionn(pos[xt+b[yt-1]][yt-1],pos[i+b[j-1]][j-1]);
}
long long ans=9,flag=0;
for(register int i=1;i<=n;i++)
if(!fa[pos[i][0]]){
if(flag)ans=ans*10%mod;
flag=1;
}
printf("%lld",ans);
return 0;
}