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

没想到这两个东西还可以这么用，我们可以用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;}