[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,131141不满足条件,前者数的长度不为6,后者第二位与第五位不同。问满足以上所有条件的数有多少个。

Solution

涨姿势了。

不难想到用并查集维护数字之间的相等关系,最后用联通块个数统计答案。

但这样的复杂度是n^2的,需要去优化它,

考虑到每次合并都是两段等长的区间进行合并,所以我们考虑使用倍增。

我们开nlogn个并查集,num[i][j]表示从i开始的2^j个数,每次区间合并我们把它拆成logn个区间分别合并。

最后自顶向底合并儿子,就像线段树一样,

Code

#include<iostream>
#include<cstdio>
#define N 100002
using namespace std;
typedef long long ll;
const int mod=1e9+;
int num[N][],f[N*],n,m,tot,son[N*][],l1,r1,l2,r2;
int find(int x){return f[x]=f[x]==x?x:find(f[x]);}
long long power(ll x,int y){
ll ans=;
while(y){
if(y&)(ans*=x)%=mod;
(x*=x)%=mod;
y>>=;
}
return ans;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=;(<<i)<=n;++i)
for(int j=;j+(<<i)-<=n;++j){
num[j][i]=++tot;f[tot]=tot;
if(i){
son[tot][]=num[j][i-];
son[tot][]=num[j+(<<i-)][i-];
}
}
for(int i=;i<=m;++i){
scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
int len=r1-l1+;
for(int j=;j>=;--j)
if((<<j)<=len){
int x=find(num[l1][j]),y=find(num[l2][j]);
if(x!=y)f[x]=y;
len-=(<<j);l1+=(<<j);l2+=(<<j);
}
}
for(int i=;i>=;--i)
for(int j=;j+(<<i)-<=n;++j){
int root=num[j][i];
if(find(root)!=root){
int x=find(son[root][]),y=find(son[f[root]][]);
if(x!=y)f[x]=y;
x=find(son[root][]),y=find(son[f[root]][]);
if(x!=y)f[x]=y;
}
}
int ans=;
for(int i=;i<=n;++i)if(find(num[i][])==num[i][])ans++;
printf("%lld",*power(,ans-)%mod);
return ;
}
上一篇:json前后台传值


下一篇:XHR——XMLHttpRquest对象