bzoj4569【SCOI2016】萌萌哒

4569: [Scoi2016]萌萌哒

Time Limit: 10 Sec   Memory Limit: 256 MB
Submit: 403   Solved: 176
[ Submit][ Status][ Discuss]

Description

一个长度为n的大数，用S1S2S3...Sn表示，其中Si表示数的第i位,S1是数的最高位，告诉你一些限制条件，每个条

件表示为四个数，l1，r1，l2，r2，即两个长度相同的区间，表示子串Sl1Sl1+1Sl1+2...Sr1与Sl2Sl2+1Sl2+2...S

r2完全相同。比如n=6时，某限制条件l1=1，r1=3，l2=4，r2=6，那么123123，351351均满足条件，但是12012，13

1141不满足条件，前者数的长度不为6，后者第二位与第五位不同。问满足以上所有条件的数有多少个。

Input

第一行两个数n和m，分别表示大数的长度，以及限制条件的个数。接下来m行，对于第i行,有4个数li1，ri1，li2

，ri2，分别表示该限制条件对应的两个区间。

1≤n≤10^5，1≤m≤10^5，1≤li1,ri1,li2,ri2≤n；并且保证ri1-li1=ri2-li2。

Output

 一个数，表示满足所有条件且长度为n的大数的个数，答案可能很大，因此输出答案模10^9+7的结果即可。

Sample Input

4 2
1 2 3 4
3 3 3 3

Sample Output

90

并查集+ST表

首先有一个直观的做法，每一个限制可以拆成O(n)个限制，然后并查集合并，最后统计答案。这样的复杂度是O(n^2)，显然是不可过接受的。

这里速度慢主要因为很多限制是重复的，所以我们要考虑减少重复计算本质上相同的限制。

用ST表实现这一过程，f[i][j]表示[i,i+2^j-1]这一区间的情况。

对于每一个限制，可以拆成两个大的限制，然后逐层向下传递，直到某一层的两个区间处于同一集合，也就说明这次限制之前已经又过了，就不需要再向下处理。

每一层最坏合并n次，一共logn层，复杂度是O(nlogn)。

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define N 100005
#define mod 1000000007
using namespace std;
int n,m,f[N][20],lg[N];
ll ans;
bool vst[N];
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x;
}
inline int find(int x,int t)
{
	return f[x][t]==x?x:f[x][t]=find(f[x][t],t);
}
inline void merge(int x,int y,int t)
{
	int fx=find(x,t),fy=find(y,t);
	if (fx==fy) return;
	f[fx][t]=fy;
	if (!t) return;
	merge(x,y,t-1);merge(x+(1<<(t-1)),y+(1<<(t-1)),t-1);
}
int main()
{
	n=read();m=read();
	for(int i=1;(1<<i)<=n;i++) lg[1<<i]=1;
	F(i,1,n) lg[i]+=lg[i-1];
	F(j,0,lg[n]) F(i,1,n) f[i][j]=i;
	F(i,1,m)
	{
		int a=read(),b=read(),c=read(),d=read();
		int len=b-a+1,t=lg[len];
		merge(a,c,t);merge(b-(1<<t)+1,d-(1<<t)+1,t);
	}
	ans=9;vst[find(1,0)]=true;
	F(i,2,n)
	{
		int x=find(i,0);
		if (!vst[x]) vst[x]=true,ans=ans*10%mod;
	}
	printf("%lld\n",ans);
	return 0;
}