bzoj 4569: [Scoi2016]萌萌哒

最新推荐文章于 2019-03-26 14:12:34 发布

heheda_is_an_OIer

最新推荐文章于 2019-03-26 14:12:34 发布

阅读量978

点赞数

本文链接：https://blog.csdn.net/heheda_is_an_oier/article/details/51372726

版权

一种会T的做法：

发现并查集加边只会发生n-1次，所以只要一个较小的复杂度完成这n-1次加边即可

启发式合并维护并查集，每次暴力修改所有属于小区间的fa值，用vector记录每个连通块内的数，暴力在线段树内修改hash值，复杂度logn,需要进行nlogn次

线段树上hash把所有的fa hash起来，然后二分找位置，每次寻找log^2,需要进行On次

    
    
     
     
      
      #include<cstdio>
     
     
     
     
      
      #include<cstdlib>
     
     
     
     
      
      #include<cstring>
     
     
     
     
      
      #include<algorithm>
     
     
     
     
      
      #include<vector> 
     
     
     
     
      
      #define ll long long
     
     
     
     
      
      #define unl unsigned long long
     
     
     
     
      
      #define inf 1e9
     
     
     
     
      
      #define eps 1e-8
     
     
     
     
      
      #define md 1000000007
     
     
     
     
      
      #define N 100010
     
     
     
     
      
      using namespace std;
     
     
     
     
      
      unl tmp[4*N],mi[N],ha[4*N];
     
     
     
     
      
      int fa[N];
     
     
     
     
      
      int lx[4*N],rx[4*N];
     
     
     
     
      
      vector<int> vec[N];
     
     
     
     
      
      const unl cheng=23333333;
     
     
     
     
      
      void build(int i,int l,int r)
     
     
     
     
      
      {
     
     
     
     
      
      	tmp[i]=mi[r-l+1]; lx[i]=l; rx[i]=r;
     
     
     
     
      
      	if (l==r) { ha[i]=fa[l]; return;}
     
     
     
     
      
      	int mid=(l+r)>>1;
     
     
     
     
      
      	build(i<<1,l,mid); build(i<<1|1,mid+1,r);
     
     
     
     
      
      	ha[i]=ha[i<<1]*tmp[i<<1]+ha[i<<1|1];
     
     
     
     
      
      }
     
     
     
     
      
       
     
     
     
     
      
      void modify(int i,int x)
     
     
     
     
      
      {
     
     
     
     
      
      	if (lx[i]==rx[i]) { ha[i]=fa[x]; return;}
     
     
     
     
      
      	int mid=(lx[i]+rx[i])>>1;
     
     
     
     
      
      	if (x<=mid) modify(i<<1,x); else modify(i<<1|1,x);
     
     
     
     
      
      	ha[i]=ha[i<<1]*tmp[i<<1]+ha[i<<1|1];
     
     
     
     
      
      }
     
     
     
     
      
       
     
     
     
     
      
      unl query(int i,int ql,int qr)
     
     
     
     
      
      {
     
     
     
     
      
      	if (ql<=lx[i]&&rx[i]<=qr) return ha[i];
     
     
     
     
      
      	int mid=(lx[i]+rx[i])>>1;
     
     
     
     
      
      	if (qr<=mid) return query(i<<1,ql,qr);
     
     
     
     
      
      	if (mid+1<=ql) return query(i<<1|1,ql,qr);
     
     
     
     
      
      	return query(i<<1,ql,mid)*mi[mid-ql+1]+query(i<<1|1,mid+1,qr);
     
     
     
     
      
      }
     
     
     
     
      
       
     
     
     
     
      
      int main()
     
     
     
     
      
      {
     
     
     
     
      
      	int n,m;
     
     
     
     
      
      	scanf("%d%d",&n,&m);
     
     
     
     
      
      	for (int i=1;i<=n;i++) { fa[i]=i; vec[i].push_back(i);}
     
     
     
     
      
      	mi[0]=1; for (int i=1;i<=n;i++) mi[i]=mi[i-1]*cheng;
     
     
     
     
      
      	build(1,1,n);
     
     
     
     
      
      	for (int i=1;i<=m;i++)
     
     
     
     
      
      	{
     
     
     
     
      
      		int l1,r1,l2,r2;
     
     
     
     
      
      		scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
     
     
     
     
      
      		while (query(1,l1,r1)!=query(1,l2,r2))
     
     
     
     
      
      		{
     
     
     
     
      
      			int L=0,R=r1-l1;
     
     
     
     
      
      			while (L!=R)
     
     
     
     
      
      			{
     
     
     
     
      
      				int mid=(L+R)>>1;
     
     
     
     
      
      				if (query(1,l1,l1+mid)!=query(1,l2,l2+mid)) R=mid; else L=mid+1;
     
     
     
     
      
      			}
     
     
     
     
      
      			int f1=fa[l1+L],f2=fa[l2+L];
     
     
     
     
      
      			if (vec[f1].size()<vec[f2].size()) swap(f1,f2);
     
     
     
     
      
      			for (int j=vec[f2].size()-1;j>=0;j--)
     
     
     
     
      
      			{
     
     
     
     
      
      				int x=vec[f2][j];
     
     
     
     
      
      				fa[x]=f1; modify(1,x);
     
     
     
     
      
      				vec[f1].push_back(x);
     
     
     
     
      
      			}
     
     
     
     
      
      			vec[f2].clear();
     
     
     
     
      
      		}
     
     
     
     
      
      	}
     
     
     
     
      
      	int k=0;
     
     
     
     
      
      	for (int i=1;i<=n;i++) if (fa[i]==i) k++;
     
     
     
     
      
      	ll ans=9;
     
     
     
     
      
      	for (int i=2;i<=k;i++) ans=ans*10%md;
     
     
     
     
      
      	printf("%lld\n",ans);
     
     
     
     
      
      	return 0;
     
     
     
     
      
      }

正解:

又是ST表。用ST表记录每一层的并查集 fa[j][i]的含义是从i开始的2^j个字符分别和从fa[j][i]开始的2^j个字符对应相等。

每次并查集合并fa[j][x],fa[j][y]的时候，要把fa[j-1][x],fa[j-1][y] 以及 fa[j-1][x+(1<<j-1)] fa[j-1][y+(1<<(j-1)] 合并

总共发生nlogn次合并，复杂度nlogn

%能想出来的人

    
    
     
     
      
      #include<cstdio>
     
     
     
     
      
      #include<cstdlib>
     
     
     
     
      
      #include<cstring>
     
     
     
     
      
      #include<algorithm>
     
     
     
     
      
      #include<cmath>
     
     
     
     
      
      #define ll long long
     
     
     
     
      
      #define unl unsigned long long
     
     
     
     
      
      #define inf 1e9
     
     
     
     
      
      #define eps 1e-8
     
     
     
     
      
      #define md 1000000007
     
     
     
     
      
      #define N 100010
     
     
     
     
      
      using namespace std;
     
     
     
     
      
      int fa[20][N];
     
     
     
     
      
      int find(int j,int i) { return fa[j][i]==i?i:fa[j][i]=find(j,fa[j][i]);}
     
     
     
     
      
       
     
     
     
     
      
      void merge(int j,int x,int y)
     
     
     
     
      
      {
     
     
     
     
      
      	int f1=find(j,x),f2=find(j,y);
     
     
     
     
      
      	if (f1==f2) return;
     
     
     
     
      
      	fa[j][f2]=f1;
     
     
     
     
      
      	if (!j) return;
     
     
     
     
      
      	merge(j-1,x,y); merge(j-1,x+(1<<(j-1)),y+(1<<(j-1)));
     
     
     
     
      
      }
     
     
     
     
      
       
     
     
     
     
      
      int main()
     
     
     
     
      
      {
     
     
     
     
      
      	int n,m;
     
     
     
     
      
      	scanf("%d%d",&n,&m);
     
     
     
     
      
      	if (n==1) { printf("10\n"); return 0;}
     
     
     
     
      
      	for (int j=0;j<18;j++)
     
     
     
     
      
      	  for (int i=1;i<=n;i++)
     
     
     
     
      
      		fa[j][i]=i;
     
     
     
     
      
      	for (int i=1;i<=m;i++)
     
     
     
     
      
      	{
     
     
     
     
      
      		int l1,r1,l2,r2;
     
     
     
     
      
      		scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
     
     
     
     
      
      		if (l1>l2) { swap(l1,l2); swap(r1,r2);}
     
     
     
     
      
      		int j=log2(r1-l1+1+eps);
     
     
     
     
      
      		merge(j,l1,l2); merge(j,r1-(1<<j)+1,r2-(1<<j)+1);
     
     
     
     
      
      	}
     
     
     
     
      
      	int ans=0;
     
     
     
     
      
      	for (int i=1;i<=n;i++)
     
     
     
     
      
      		if (find(0,i)==i) ans++;
     
     
     
     
      
      	ll ret=9;
     
     
     
     
      
      	for (int i=2;i<=ans;i++) ret=ret*10%md;
     
     
     
     
      
      	printf("%lld\n",ret);
     
     
     
     
      
      	return 0;
     
     
     
     
      
      }