[BZOJ4569][Scoi2016]萌萌哒（并查集+ST表）

=== ===

这里放传送门

=== ===

题解

显然这玩意儿它给出的两个区间前面的区间会限制后面的区间。但是有时候后面的区间也会反过来限制前面的区间，也就是说这个限制关系是有传递性的，可以想到用并查集来维护，对于每个点维护一下它和别的点的连通关系，读入一组约束的时候用O(n)的时间合并。但是这种做法显然很暴力啊是过不去的。。可以发现读入约束的时间是必须的，那要优化只能优化合并的时间。于是这里就有了一个神奇的思路就是把并查集套上ST表，倍增地来用！

用f[x][i]表示x这个点开头往后的$2^i$个位置，比如说如果f[x][i]和f[y][i]合并了就说明x开头的$2^i$个位置和y开头的$2^i$个位置是对应相同的。这里要区分“对应相同”和“完全一样”，不是说这$2*2^i$个字符都是一模一样的而是说它们对应位置相等，比如x和y位置上的字符相等，x+1和y+1位置上的字符相等那样。

但是这样的话怎么最后统计有多少个可以独立填数字的位置呢？也就是说，用这样类似ST表的形式合并了以后怎么查询单个位置的信息呢？这就要求每次合并以后要下放信息，就是说要把长度为$2^i$的大区间下面管辖的所有小区间也一个一个合并。为了保证时间复杂度，遇到已经合并过的区间就要及时退出，这样可以保证每个点只被合并一次，时间复杂度仍然是$O(nlogn)$的。

代码

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define Pow 17
using namespace std;
const int Mod=1e9+7;
int n,m,fa[20][100010],cnt;
long long ans;
int find(int p,int x){
    if (fa[p][x]==x) return fa[p][x];
    fa[p][x]=find(p,fa[p][x]);
    return fa[p][x];
}
void merge(int p,int x,int y){
    int r1,r2;
    r1=find(p,x);r2=find(p,y);
    if (r1==r2) return;
    fa[p][r2]=r1;
    if (p==0) return;
    merge(p-1,x,y);//递归合并较小的区间
    merge(p-1,x+(1<<(p-1)),y+(1<<(p-1)));
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=Pow;i>=0;i--)
      for (int j=1;j<=n;j++)
        fa[i][j]=j;
    for (int i=1;i<=m;i++){
        int x1,y1,x2,y2,j;
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        if (x1>x2){swap(x1,x2);swap(y1,y2);}
        j=floor(log2(y1-x1+1));
        merge(j,x1,x2);merge(j,y1-(1<<j)+1,y2-(1<<j)+1);//确保区间恰好被覆盖
    }
    for (int i=1;i<=n;i++)
      if (find(0,i)==i) ++cnt;
    ans=9;
    for (int i=2;i<=cnt;i++) ans=(ans*10)%Mod;
    printf("%d\n",ans);
    return 0;
}

偏偏在最后出现的补充说明

这玩意儿该叫什么呢？倍增并查集？好新鲜啊以前从来没见过诶。
话说回来，有时候均摊复杂度的算法还是要好好考虑一下诶。看起来很不科学实际上是很不错的对吧。
倍增的思想还是需要注意一下。