bzoj 4569（并查集）

传送门
讲解来自：http://blog.csdn.net/werkeytom_ftd/article/details/51334030
我们用ST表，f[i,j]表示[i,i+2^j-1]这一段。
那么初始时每一段单独成一个集合。
对于一个限制可以拆成log 份，然后进行集合合并。
然后呢，如果任意f[s,t]和f[i,j]属于同一集合，那么f[s,t-1]与f[i,j-1]以及f[s+2^(t-1)-1,t-1]和f[i+2^(j-1)-1,j-1]都应该属于同一集合。
为了满足这个限制，只要最后再一层一层的做，把下一层的合并了即可。
合并注意必须让编号大的合进编号小的里。
统计答案就很简单啦，答案是9*10^(集合个数-1)。
莫名其妙的因为变量定义顺序没对？！RE了半小时，对于题目本身就不死死纠结了，方法学会就好：http://blog.csdn.net/KGV093/article/details/78039328

#include<bits/stdc++.h>
using namespace std;
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
const int MAXN=1e5+2,MOD=1e9+7;
int n,m;
int fa[MAXN*20],id[MAXN][20],mk[MAXN*20][2],pw[20],tot=0,cnt=0,ans;
int find(int x){
    return fa[x]?fa[x]=find(fa[x]):x;
}
int merge(int x,int y){
    x=find(x);
    y=find(y);
    if (x!=y) fa[y]=x;
}
inline int read() {
    int x=0;char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x;
}
int main(){
    pw[0]=1;
    fo(i,1,17) pw[i]=pw[i-1]*2;
    n=read(),m=read();
    if (n==1){
        printf("10\n");
        return 0;
    }
    fo(i,1,n)
        fo(j,0,17)
            id[i][j]=++tot,mk[tot][0]=i,mk[tot][1]=j;
    fo(i,1,m){
        int l1=read(),r1=read(),l2=read(),r2=read();
        fd(j,17,0)
            if (l1+pw[j]-1<=r1){
                merge(id[l1][j],id[l2][j]);
                l1+=pw[j];l2+=pw[j];
            }
    }
    fd(j,17,1)
        fo(i,1,n){
            int k=find(id[i][j]);
            int s=mk[k][0];
            merge(id[s][j-1],id[i][j-1]);
            merge(id[s+pw[j-1]][j-1],id[i+pw[j-1]][j-1]);
        }
    fo(i,1,n)
        if (!fa[id[i][0]]) cnt++;
    ans=9;
    fo(i,1,cnt-1) ans=(long long)ans*10%MOD;
    printf("%d\n",ans);
}