4569: [Scoi2016]萌萌哒
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 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
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;
}