zoukankan      html  css  js  c++  java
  • 洛谷P3295 [SCOI2016]萌萌哒(倍增+并查集)

    传送门

    思路太妙了啊……

    容易才怪想到暴力,把区间内的每一个数字用并查集维护相等,然后设最后总共有$k$个并查集,那么答案就是$9*10^{k-1}$(因为第一位不能为0)

    考虑倍增。我们设$f[i][j]$表示区间$[i,i+2^j-1]$,那么我们可以把原区间给拆成$log$个区间,然后维护这些区间的连通性

    然而我们最后需要的是最底层的,也就是单独的节点的连通性。那么我们考虑如何将连通性向下传递。如果$f[i][j]$和$f[a][b]$连通,那么$f[i][j-1]$和$f[a][b-1]$一定连通(前半部分区间),$f[i+2^{j-1}][j-1]$和$f[a+2^{b-1}][b-1]$也一定连通

    ps:连通性肯定都在同一层,所以实际上上面的$j$和$b$一般都是相等的

    然后只要最后判最底层有几个并查集就好了

     1 //minamoto
     2 #include<iostream>
     3 #include<cstdio>
     4 using namespace std;
     5 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
     6 char buf[1<<21],*p1=buf,*p2=buf;
     7 int read(){
     8     #define num ch-'0'
     9     char ch;bool flag=0;int res;
    10     while(!isdigit(ch=getc()))
    11     (ch=='-')&&(flag=true);
    12     for(res=num;isdigit(ch=getc());res=res*10+num);
    13     (flag)&&(res=-res);
    14     #undef num
    15     return res;
    16 }
    17 const int N=1e5+5,mod=1e9+7;
    18 int fa[N*17],id[N][21],num[N*17],log[N*17],bin[21],is[N*17],h[N*17],tot=0,cnt;
    19 int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
    20 void merge(int x,int y){
    21     x=find(x),y=find(y);
    22     if(h[x]<h[y]) fa[x]=y;
    23     else if(h[x]>h[y]) fa[y]=x;
    24     else fa[x]=y,++h[y];
    25 }
    26 int ksm(int b){
    27     int res=9,a=10;
    28     while(b){
    29         if(b&1) res=1ll*res*a%mod;
    30         a=1ll*a*a%mod,b>>=1;
    31     }
    32     return res;
    33 }
    34 int main(){
    35 //    freopen("testdata.in","r",stdin);
    36     int n=read(),m=read();
    37     bin[0]=1;for(int i=1;i<=16;++i) bin[i]=bin[i-1]<<1;
    38     for(int j=0;j<=16;++j) for(int i=1;i<=n;++i) id[i][j]=++tot,num[tot]=i,fa[tot]=tot,h[tot]=1;
    39     while(m--){
    40         int l1=read(),r1=read(),l2=read(),r2=read();
    41         for(int i=16;i>=0;--i) if(l1+bin[i]-1<=r1){
    42             merge(id[l1][i],id[l2][i]),l1+=bin[i],l2+=bin[i];
    43         }
    44     }
    45     for(int j=16;j;--j) for(int i=1;i+bin[j]-1<=n;++i){
    46         int x=find(id[i][j]),a=num[x];
    47         merge(id[a][j-1],id[i][j-1]),merge(id[a+bin[j-1]][j-1],id[i+bin[j-1]][j-1]);
    48     }
    49     for(int i=1;i<=n;++i)
    50     if(find(id[i][0])==id[i][0]) ++cnt;
    51     printf("%d
    ",ksm(cnt-1));
    52     return 0;
    53 }
  • 相关阅读:
    第12组 Beta冲刺 (4/5)
    第12组 Beta冲刺 (3/5)
    代码用辅助多了 基础的读取config都忘记了
    wpf 动态添加控件 通过xmal实习 c#代码插入控件
    C#里调用非托管的Dll -z注意 net版本
    动态调用 类库
    c#时间的生成
    c# 第三方 修改版本号 MSBuildTasks, 解决 通配符不匹配问题
    c#.exe以管理员身份运行
    log4
  • 原文地址:https://www.cnblogs.com/bztMinamoto/p/9812565.html
Copyright © 2011-2022 走看看