BZOJ 4569 萌萌哒 （并查集）

4569: [Scoi2016]萌萌哒

Time Limit: 10 Sec Memory Limit: 256 MB
Description

一个长度为n的大数，用S1 S2 S3… Sn表示，其中Si表示数的第i位,S1是数的最高位，告诉你一些限制条件，每个条
件表示为四个数，l1，r1，l2，r2，即两个长度相同的区间，表示子串Sl1 Sl1+1 Sl1+2… Sr1与Sl2 Sl2+1 Sl2+2 … Sr2完全相同。比如n=6时，某限制条件l1=1，r1=3，l2=4，r2=6，那么123123，351351均满足条件，但是12012，131141不满足条件，前者数的长度不为6，后者第二位与第五位不同。问满足以上所有条件的数有多少个。

Input

第一行两个数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。

Output

一个数，表示满足所有条件且长度为n的大数的个数，答案可能很大，因此输出答案模10^9+7的结果即可。

Sample Input

4 2
1 2 3 4
3 3 3 3

Sample Output

90

思路：
因此使用ST表维护并查集，从高到低进行类似ST表合并的操作，最后判定有多少个代表元素，则每个并查集内的数字一定相同，简单的乘法原理即可得到答案。

#include<cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
using namespace std;

const int N = 1e5 + 5, MOD = 1e9+7, P = 17;
int n,m;
int fa[N*20],f[N][20],mk[N*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 fpow(int a, int b, int p) {
    int ret = 1;
    while ( b ) {
        if(b & 1) ret = (LL)ret * a  %p;
        b >>= 1, a = (LL)a * a % p;
    }
    return ret;
}

int main() {
    pw[0] = 1;
    for(int i=1; i<=P; ++i) pw[i] = pw[i-1] << 1;
    n=read() , m=read();
    if(n == 1) ans = 10;
    else{
        for(register int i=1; i<=n; ++i)
            for(int j=0; j<=P; ++j)
                f[i][j] = ++tot, mk[tot][0] = i, mk[tot][1] = j;//f[i,j]表示[i,i+2^j-1]这一段 
        for(int i=1; i<=m; ++i) {//那么初始时每一段单独成一个集合 
            int l1=read(), r1=read(), l2=read(), r2=read();
            for(int j=P; j>=0; --j)
                if(l1 + pw[j] - 1 <= r1) {
                    merge(f[l1][j], f[l2][j]);
                    l1 += pw[j]; l2 += pw[j];//对于一个限制可以拆成log 份，然后进行集合合并 
            }
        }
        for(int j=P; j>0; --j)
            for(register int i=1; i<=n; ++i) {
                int k = find(f[i][j]);
                int s = mk[k][0];
                merge(f[s][j-1], f[i][j-1]);
                merge(f[s+pw[j-1]][j-1], f[i+pw[j-1]][j-1]);
            }//如果任意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]都应该属于同一集合。
        for(register int i=1; i<=n; ++i)
            if( !fa[f[i][0]] ) ++cnt;
        ans = 9;//不含前导0 
        for(register int i=1; i<cnt; ++i) ans = (LL)ans * 10 % MOD;
    }
    printf("%lld\n", (LL)ans);
    return 0;
}