并查集奇妙应用-机房模拟赛某题-同BZOJ4569

黑白熊的密码

256MB / 1s ; monokuma.cpp / c / pas / in / out

题目描述

苗木诚发现了脱出装置，它可以帮助所有人逃脱希望之峰学园。但是要想使用它，必须输入一个十进制的n位数密码（没有前导0）。为了让游戏变得更加有趣，黑白熊提供了m条信息，其中第i条信息指出这个数的第Li位到Ri位和第li位到ri位完全相同。
苗木诚决定随机选择一个满足要求的n位数输入。请你告诉他他碰对密码的概率是多少。

输入格式

第一行两个整数n,m，表示数的位数和黑白熊提供的信息数。
接下来m行，每行四个正整数Li,Ri,li,ri。

输出格式

输出一行一个整数，表示苗木诚碰对密码的概率，对998244353取模。

样例数据

monokuma.in monokuma.out
4 2
1 2 3 4
3 3 3 3 144190851

样例解释

答案是1/90在模998244353意义下的值。

数据范围

对于30%的数据，n,m≤2000。
对于100%的数据，1≤n,m≤100000，1≤Li≤Ri≤n，1≤li≤ri≤n，且保证Ri-Li=ri-li。

---

这题，其实不难看出来要使用并查集
如果有两个区间相同，那么这两个区间中的每个对应位置都相同，相当于就少了那么多位（反正都是一样的嘛，那些关联起来的最后只留一个就可以了），于是只需要用并查集挨个挨个的合并就好了XD。
然而并没有那么水…但是这样会明显会GG！
数据规模那么大怎么可能过嘛qwq，naive！

这时候，一个谜一样的东西出现了！倍增思想
因为每次出现一个区间我们就要挨个挨个合并太麻烦，那我们不如用类似ST表的那种思想，维护log个并查集。表示从某位开始往后的几位都相同，这样的话每个条件就可以O(2)处理完，最后将高位的区间下放下去就好了（这个时候再统一的一个一个合并）

然后就可以轻松愉快的A题啦qwq

---

自带大常数的代码

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

const int mmod = 998244353 ;
int N , M , lg[210005] ;
struct Union_set{
    int b[100005] ;
    void clear(){
        for( int i = 1 ; i <= N ; i ++ )
            b[i] = i ;
    }
    int find( int x ){
        if( b[x] == x ) return x ;
        return b[x] = find( b[x] ) ;
    }
    bool Union( int x , int y ){
        int F_x = find( x ) , F_y = find( y ) ;
        if( F_x == F_y ) return false ;
        b[F_x] = F_y ;
        return true ;
    }
}U[17] ;

long long s_pow( long long x ){
    int b = mmod - 2 ;
    long long rt = 1 ;
    while( b ){
        if( b&1 )   rt = rt * x %mmod ;
        x = x * x %mmod ; b >>= 1 ;
    }
    return rt ;
}

void solve(){
    int cnt = 0 ;
    long long fm = 9 ;
    for( int i = 1 , L , R , lf , rg ; i <= M ; i ++ ){
        scanf( "%d%d%d%d" , &L , &R , &lf , &rg ) ;
        int LOG = lg[ R - L + 1 ] ;
        U[LOG].Union( L , lf ) ;
        U[LOG].Union( R - (1<<LOG) + 1 , rg - (1<<LOG) + 1 ) ;
    }
    //printf( "%d\n" , lg[N] ) ; 
    for( int i = lg[N] ; i >= 1 ; i -- ){
        //printf( "i = %d::\n" , i );
        for( int j = 1 ; j + ( 1<<i ) - 1 <= N ; j ++ ){
            int F_j = U[i].find( j ) ;
            U[i-1].Union( j , F_j ) ;
            U[i-1].Union( j + ( 1<<(i-1) ) , F_j + ( 1<<(i-1) ) ) ;
            //printf( "%d U %d\n%d U %d\n" , j , F_j , j + ( 1<<(i-1) ) , F_j + ( 1<<(i-1) ) ) ;
        }
    }
    for( int i = 1 ; i <= N ; i ++ )
        if( i == U[0].find( i ) ) cnt ++ ;
    for( int i = 1 ; i < cnt ; i ++ )   fm = fm * 10%mmod ;
    printf( "%I64d" , s_pow( fm ) ) ;
}

int main(){
    freopen( "monokuma.in" , "r" , stdin ) ;
    freopen( "monokuma.out", "w" , stdout) ;
    scanf( "%d%d" , &N , &M ) ;
    for( int i = 0 ; i <= 16 ; i ++ )
        U[i].clear() ;
    for( int i = 4 ; i <= ( N << 1 ) ; i <<= 1 )
        lg[i] = lg[i/2] + 1 ;
    for( int i = N * 2 ; i > 2 ; i -- )
        if( !lg[i] ) lg[i] = lg[i+1] ;
    /*for( int i = 1 ; i <= N ; i ++ )
        printf( "%d " , lg[i] ) ;
    printf( "\n" ) ;*/
    solve() ;
}