BZOJ 4569: [Scoi2016]萌萌哒 ST表并查集

最新推荐文章于 2021-08-11 11:43:06 发布

原创最新推荐文章于 2021-08-11 11:43:06 发布 · 459 阅读

0 ·

CC 4.0 BY-SA版权

—————————基础数据结构同时被 2 个专栏收录

43 篇文章

订阅专栏

并查集

12 篇文章

订阅专栏

本文介绍了一道名为“萌萌哒”的算法题目，通过使用ST表维护区间和区间操作的方法解决了一个关于寻找满足特定条件的大数数量的问题。该题需要考虑大数长度和限制条件，最终输出答案模10^9+7的结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

4569: [Scoi2016]萌萌哒

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 976 Solved: 473
[Submit][Status][Discuss]

Description

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

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

刚开始连并查集都没想

就想只要一个区间与别的联动就不对答案有贡献，然后乱搞就WA了

翻题解。。。

用ST表维护区间，区间操作

最后下传

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;

typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=1000100,mod=int(1e9)+7;

int f[N][20];

int find(int x,int y)
{return f[x][y]==x?x:f[x][y]=find(f[x][y],y);}

void update(int a,int b,int c)
{if(find(a,c)^find(b,c))f[f[a][c]][c]=f[b][c];}

int main()
{
	int n=read(),m=read();
	register int i,j,l,r,x,y,len;
	for(i=1;i<=n;++i)for(j=0;j<=17;++j)f[i][j]=i;
	while(m--)
	{
		l=read();r=read();x=read();y=read();
		len=r-l+1;
		for(i=0;(1<<i)<=len;++i)if((1<<i)&len)
		update(l,x,i),l+=(1<<i),x+=(1<<i);
	}
	for(j=17;j;j--)for(i=1;i<=n;++i)
	update(i,find(i,j),j-1),update(i+(1<<(j-1)),f[i][j]+(1<<(j-1)),j-1);
	x=0;
	for(i=1;i<=n;++i)if(find(i,0)==i)x++;
	y=9;
	for(i=1;i<x;++i)y=1ll*y*10%mod;
	print(y);puts("");
	return 0;
}
/*
4 2
1 2 3 4 
3 3 3 3

90
*/