bzoj 4569(并查集)

来源:互联网 发布:angularjs repeat源码 编辑:程序博客网 时间:2024/05/22 16:03

传送门
讲解来自:http://blog.csdn.net/werkeytom_ftd/article/details/51334030
我们用ST表,f[i,j]表示[i,i+2^j-1]这一段。
那么初始时每一段单独成一个集合。
对于一个限制可以拆成log 份,然后进行集合合并。
然后呢,如果任意f[s,t]和f[i,j]属于同一集合,那么f[s,t-1]与f[i,j-1]以及f[s+2^(t-1)-1,t-1]和f[i+2^(j-1)-1,j-1]都应该属于同一集合。
为了满足这个限制,只要最后再一层一层的做,把下一层的合并了即可。
合并注意必须让编号大的合进编号小的里。
统计答案就很简单啦,答案是9*10^(集合个数-1)。
莫名其妙的因为变量定义顺序没对?!RE了半小时,对于题目本身就不死死纠结了,方法学会就好:http://blog.csdn.net/KGV093/article/details/78039328

#include<bits/stdc++.h>using namespace std;#define fo(i,a,b) for(int i=a;i<=b;i++)#define fd(i,a,b) for(int i=a;i>=b;i--)const int MAXN=1e5+2,MOD=1e9+7;int n,m;int fa[MAXN*20],id[MAXN][20],mk[MAXN*20][2],pw[20],tot=0,cnt=0,ans;int find(int x){    return fa[x]?fa[x]=find(fa[x]):x;}int merge(int x,int y){    x=find(x);    y=find(y);    if (x!=y) fa[y]=x;}inline int read() {    int x=0;char c=getchar();    while (c<'0'||c>'9') c=getchar();    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();    return x;}int main(){    pw[0]=1;    fo(i,1,17) pw[i]=pw[i-1]*2;    n=read(),m=read();    if (n==1){        printf("10\n");        return 0;    }    fo(i,1,n)        fo(j,0,17)            id[i][j]=++tot,mk[tot][0]=i,mk[tot][1]=j;    fo(i,1,m){        int l1=read(),r1=read(),l2=read(),r2=read();        fd(j,17,0)            if (l1+pw[j]-1<=r1){                merge(id[l1][j],id[l2][j]);                l1+=pw[j];l2+=pw[j];            }    }    fd(j,17,1)        fo(i,1,n){            int k=find(id[i][j]);            int s=mk[k][0];            merge(id[s][j-1],id[i][j-1]);            merge(id[s+pw[j-1]][j-1],id[i+pw[j-1]][j-1]);        }    fo(i,1,n)        if (!fa[id[i][0]]) cnt++;    ans=9;    fo(i,1,cnt-1) ans=(long long)ans*10%MOD;    printf("%d\n",ans);}