bzoj 4569: [Scoi2016]萌萌哒 ST表+并查集

题意

给出若干个限制l1 r1 l2 r2表示区间[l1,r1]和[l2,r2]相同，问有多少个不含前导0的n位数可以满足所有限制，答案模10^9+7.
n<=100000

分析

一开始全往大数据结构上面想了，根本没有想到用ST表来做。

题解这个想法非常的棒棒，用f[i,j]表示以i开头长度为2^j的区间
若f[i,j]和f[k,j]属于统一集合则表示区间$[i,i+2^j-1]和[k,k+2^j-1]$相同。那么对于一个限制，我们可以将其分为logn个部分分别合并。
若$[i,i+2^j-1]$和$[k,k+2^k-1]$相等，则必有$[i,i+2^{j-1}-1]$和$[k,k+2^{j-1}-1]$相等，所以就一直往下合并即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100005
#define LL long long
#define MOD 1000000007
using namespace std;

int n,m,tot,num[N][25],fa[N*25],t[N*25];

int find(int x)
{
    if (fa[x]==x) return x;
    fa[x]=find(fa[x]);
    return fa[x];
}

void merge(int x,int y,int j)
{
    int fx=find(num[x][j]),fy=find(num[y][j]);
    if (fx!=fy) fa[fx]=fy;
    else return;
    merge(x,y,j-1);merge(x+(1<<(j-1)),y+(1<<(j-1)),j-1);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        for (int j=0;i+(1<<j)-1<=n;j++)
            num[i][j]=++tot;
    for (int i=1;i<=tot;i++) fa[i]=i;
    for (int i=1;i<=m;i++)
    {
        int l1,r1,l2,r2;
        scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
        int len=r1-l1+1,lg=log(len)/log(2),now=0;
        for (int j=lg;j>=0;j--)
            if (now+(1<<j)<=len)
            {
                merge(l1+now,l2+now,j);
                now+=(1<<j);
            }
    }
    int ans=0;
    for (int i=1;i<=n;i++)
    {
        int x=find(num[i][0]);
        if (!t[x]) ans++,t[x]=1;
    }
    int w=9;
    for (int i=1;i<ans;i++) w=(LL)w*10%MOD;
    printf("%d",w);
    return 0;
}