BZOJ4569 萌萌哒 [倍增][并查集]

没想到这两个东西还可以这么用，我们可以用$f[i][j]$表示$[i,i+2^j-1]$这一段区间属于哪一个集合,没有则等于0.
那么每次合并我们可以把这个限制拆成log个区间，依次合并起来。
若$f[i][j]$和$f[s][t]$同属一个集合，那么$f[i][j-1]$和$f[s][t-1]$，$f[i+2^{j-1}][j-1]$和$f[s+2^{t-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;
}