[SCOI2016]萌萌哒

题目描述

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

直观的想法

并查集，对应位直接并一下，然后最后看有多少个联通块，恶搞一下。

机智的正解

我们可以用ST表(就比如a[i][j]表示i位置开始的2^j长度的串)
首先对于每个条件，我们直接将对应的块并一下。
输入完然后将联通关系在ST表中下传，例如两个长度为2^i的串为同一联通块，那么各自的前后2^(i-1)分别都是一个联通块的。

贴代码：

/**************************************************************
    Problem: 4569
    User: samjia2000
    Language: C++
    Result: Accepted
    Time:4048 ms
    Memory:40444 kb
****************************************************************/

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>

#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)

using namespace std;

typedef long long LL;

const int N = 100010;

int fa[N*20];
int num[N][20];
int n,m,k;
int h[N*20],vis[N*20],next[N*20];

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

void prepare(){
    fo(i,0,log(n)/log(2))
        fo(j,1,n)fa[num[j][i]=++k]=k;
}

void merge(int x,int y,int len){
    fd(i,log(n)/log(2),0)
        if (len>=(1<<i)){
            len-=(1<<i);
            int fx=getfather(num[x][i]),fy=getfather(num[y][i]);
            fa[fx]=fy;
            x+=(1<<i);
            y+=(1<<i);
        }
}

int main(){
    scanf("%d%d",&n,&m);
    prepare();
    fo(i,1,m){
        int l1,r1,l2,r2;
        scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
        merge(l1,l2,r1-l1+1);
    }
    fd(x,log(n)/log(2),1){
        fo(i,1,n){
            h[num[i][x]]=next[num[i][x]]=0;
            getfather(num[i][x]);
        }
        fo(i,1,n){
            next[i]=h[fa[num[i][x]]];
            h[fa[num[i][x]]]=i;
        }
        fo(i,1,n)
        if (vis[fa[num[i][x]]]!=x){
            vis[fa[num[i][x]]]=x;
            int y=h[fa[num[i][x]]];
            int f1=getfather(num[y][x-1]),f2=getfather(num[y+(1<<(x-1))][x-1]);
            for(y=next[y];y;y=next[y]){
                int x1=getfather(num[y][x-1]);
                f1=getfather(f1);
                fa[x1]=f1;
                x1=getfather(num[y+(1<<(x-1))][x-1]);
                f2=getfather(f2);
                fa[f2]=x1;
            }
        }
    }
    fo(i,1,n)getfather(num[i][0]);
    int tot=0,ans=9;
    fo(i,1,n)
    if (fa[num[i][0]]==num[i][0])tot++;
    fo(i,2,tot)ans=LL(ans)*10%1000000007;
    printf("%d\n",ans);
    return 0;
}