[BZOJ4569] [Luogu 3295] [SCOI2016]萌萌哒(并查集+倍增)
[BZOJ4569] [Luogu 3295] [SCOI2016]萌萌哒(并查集+倍增)
题面
有一个n位的十进制数a(无前导0),给出m条限制,每条限制\((l_1,r_1,l_2,r_2)(保证r_1-l_1=r_2-l_2)\)表示这个数的第\([l_1,r_1]\)位与\([l_2,r_2]\)位相同。问有多少个这样的数满足条件,答案取模\(10^9+7\),
\(n \leq 10^5\)
分析
约定:我们把十进制数a的每一位从高到低称为第1位,第2位...,记为\(a_i\)
首先很妙的一点是,把区间限制转化成元素相等,即\(a_{l_1+k}=a_{r_1+k}(k \in [0,r_1-l_1+1])\)。然后把所有相等的元素看作一个联通块,对于每条限制暴力合并,最后我们就会得到\(cnt\)个联通块。对于每个联通块,我们要给它的所有元素赋一个相同的值。包含\(a_1\)的联通块里的元素只能在\([1,9]\)内取值(因为第一位不能是0),其他联通块的元素可以在\([0,9]\)内取值。根据乘法原理,答案就是\(9 \times 10^{cnt-1}\)
那么问题就是如何快速维护连通性。直接做是\(O(n^2\alpha(n))\)的。注意到我们一个个合并\([l_1,r_1]\)内的所有元素,很浪费时间,可以用倍增优化。
我们将并查集的1维改成2维,\(f[i][j]\)表示区间\([j,j+2^i-1]\)的代表元。合并的时候我们只需要把\(r_1-l_1+1\)二进制拆分,如13拆成\(2^3+2^2+2^0\).对于每一位\(j\),合并\(f[j]\)上对应的一段长度为\(2^j\)的区间。
int len=r1-l1+1;
for(int j=log2n;j>=0;j--){//二进制拆分
if(len&(1<<j)){
S.merge(j,l1,l2);
l1+=(1<<j);
l2+=(1<<j);
}
}
[图片来自mhy博客]
如下图,3的二进制拆分为\(2^1+2^0\)先合并长度为\(2^1\)的区间,再合并长度为\(2^0\)的区间
但是我们查询的时候需要查询的是每个点的信息,即\(f[0][i]\),所以我们要下推合并。合并长度为\(2^j\),起点分别为\(x,y\)的区间时,我们把区间断成两半,合并\(f[j-1][x+2^{j-1}],f[j-1][y+2^{j-1}]\),还要合并\(f[j-1][x],f[j-1][y]\)这样做\(O(\log n)\)次,就可以得到\(f[0][i]\)
for(int j=log2n;j>=1;j--){
for(int i=1;i+(1<<j)-1<=n;i++){//下推关系,长度为2^j的推到2^(j-1)
int f=S.find(j,i);
S.merge(j-1,i,f);
S.merge(j-1,i+(1<<(j-1)),f+(1<<(j-1)));
}
}
时间复杂度为\(O(n \alpha(n)\log n )\),但是这里我偷了个懒,只写了路径压缩的并查集,时间复杂度\(O(n \log^2 n)\),由于常数很小,可以通过。
代码
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#define maxn 100000
#define maxlogn 21
#define mod 1000000007
using namespace std;
typedef long long ll;
inline ll fast_pow(ll x,ll k){
ll ans=1;
while(k){
if(k&1) ans=ans*x%mod;
x=x*x%mod;
k>>=1;
}
return ans;
}
struct disjoint_set{
int fa[maxlogn+5][maxn+5];//f[i][j]表示[j,j+2^i-1]的代表元
void ini(int n){
int k=log2(n)+1;
for(int j=0;j<=k;j++){
for(int i=1;i<=n;i++){
fa[j][i]=i;
}
}
}
int find(int b,int x){
if(fa[b][x]==x) return x;
else return fa[b][x]=find(b,fa[b][x]);
}
void merge(int b,int x,int y){
int fx=find(b,x);
int fy=find(b,y);
if(fx!=fy) fa[b][fx]=fy;
}
}S;
int log2n;
int n,m;
int main(){
int l1,l2,r1,r2;
scanf("%d %d",&n,&m);
log2n=log2(n)+1;
S.ini(n);
for(int i=1;i<=m;i++){
scanf("%d %d %d %d",&l1,&r1,&l2,&r2);
int len=r1-l1+1;
for(int j=log2n;j>=0;j--){//二进制拆分
if(len&(1<<j)){
S.merge(j,l1,l2);
l1+=(1<<j);
l2+=(1<<j);
}
}
}
for(int j=log2n;j>=1;j--){
for(int i=1;i+(1<<j)-1<=n;i++){//下推关系,长度为2^j的推到2^(j-1)
int f=S.find(j,i);
S.merge(j-1,i,f);
S.merge(j-1,i+(1<<(j-1)),f+(1<<(j-1)));
}
}
int cnt=0;
for(int i=1;i<=n;i++){
if(S.find(0,i)==i) cnt++;
}
printf("%lld\n",9*fast_pow(10,cnt-1)%mod);
}
[BZOJ4569] [Luogu 3295] [SCOI2016]萌萌哒(并查集+倍增)的更多相关文章
- Luogu P3295 [SCOI2016]萌萌哒(并查集+倍增)
P3295 [SCOI2016]萌萌哒 题面 题目描述 一个长度为 \(n\) 的大数,用 \(S_1S_2S_3 \cdots S_n\) 表示,其中 \(S_i\) 表示数的第 \(i\) 位, ...
- 洛谷 3295 [SCOI2016]萌萌哒——并查集优化连边
题目:https://www.luogu.org/problemnew/show/P3295 当要连的边形如 “一段区间内都是 i 向 i+L 连边” 的时候,用并查集优化连边. 在连边的时候,如果要 ...
- [bzoj4569][SCOI2016]萌萌哒-并查集+倍增
Brief Description 一个长度为n的大数,用S1S2S3...Sn表示,其中Si表示数的第i位,S1是数的最高位,告诉你一些限制条件,每个条 件表示为四个数,l1,r1,l2,r2,即两 ...
- BZOJ 4569: [Scoi2016]萌萌哒 [并查集 倍增]
传送门 题意:长为$n \le 10^5$的数字,给出$m \le 10^5$个限制$[l1,r1]\ [l2,r2]$两个子串完全相等,求方案数 把所有要求相等的位置连起来,不就是$9*10^{连通 ...
- LOJ2014 SCOI2016 萌萌哒 并查集、ST表优化连边
传送门 一个朴素的做法就是暴力连边并查集,可是这是\(O(n^2)\)的.发现每一次连边可以看成两个区间覆盖,这两个区间之间一一对应地连边.可线段树对应的两个节点的size可能不同,这会导致" ...
- bzoj 4569 [Scoi2016]萌萌哒 并查集 + ST表
题目链接 Description 一个长度为\(n\)的大数,用\(S_1S_2S_3...S_n\)表示,其中\(S_i\)表示数的第\(i\)位,\(S_1\)是数的最高位,告诉你一些限制条件,每 ...
- BZOJ4569 [SCOI2016]萌萌哒 【并查集 + 倍增】
题目链接 BZOJ4569 题解 倍增的思想很棒 题目实际上就是每次让我们合并两个区间对应位置的数,最后的答案\(ans = 9 \times 10^{tot - 1}\),\(tot\)是联通块数, ...
- 2018.07.31 bzoj4569: [Scoi2016]萌萌哒(并查集+倍增)
传送门 对于每个限制,使用倍增的二进制拆分思想,用并查集数组fa[i][j]" role="presentation" style="position: rel ...
- luogu P3295 [SCOI2016]萌萌哒
传送门 题目条件"两个子串\(S[l_1,r_1],S[l_2,r_2]\)完全相同"等价于\(\forall i \in[0,r_1-l_1+1],S_{l1+i}=S_{l_2 ...
随机推荐
- scipy.sparse 稀疏矩阵
from 博客园(华夏35度)http://www.cnblogs.com/zhangchaoyang 作者:Orisun 本文主要围绕scipy中的稀疏矩阵展开,也会介绍几种scipy之外的稀疏矩阵 ...
- 大数乘法(A * B Problem Plus)问题
大数乘法问题一般可以通过将大数转换为数组来解决. 解题思路 第1步 第2步 第3步 第4步 样例输入1 56 744 样例输出1 800 样例输入2 -10 678 样例输出2 -6780 样例输入3 ...
- MTBF测试
本文转载自:https://blog.csdn.net/liuhaoemail/article/details/50531489 MTBF测试 目前,终端侧的可靠性测试基本上是采用称为”MTBF测试” ...
- 【gym102394B】Binary Numbers(DP)
题意:From https://blog.csdn.net/m0_37809890/article/details/102886956 思路: 可以发现转移就是右上角的一个区间前缀和 std只要开1倍 ...
- Spring Boot教程(十六)属性配置文件详解(1)
相信很多人选择Spring Boot主要是考虑到它既能兼顾Spring的强大功能,还能实现快速开发的便捷.我们在Spring Boot使用过程中,最直观的感受就是没有了原来自己整合Spring应用时繁 ...
- sqli-lab(16)
现实证明 英语好才能学渗透 基于bool类型的时间盲注 双引号 0X01爱之初体验 说实话我还没有找到 盲注的时候怎么判断闭合的方法 so 这里我直接看源码 这里的语句应该是 (“”)这种的闭合形式 ...
- JDK中String类的源码分析(二)
1.startsWith(String prefix, int toffset)方法 包括startsWith(*),endsWith(*)方法,都是调用上述一个方法 public boolean s ...
- 64位 Qt5.12 MySql 连接问题
关于怎么检查Qt是否带MySql驱动 ,到Qt安装目录下 plugins\sqldrivers下寻找是否有qsqlmysql.dll文件 例如:F:\Qt\Qt5.9.6\5.9.6\msv ...
- 【机器学习速成宝典】模型篇02线性回归【LR】(Python版)
目录 什么是线性回归 最小二乘法 一元线性回归 多元线性回归 什么是规范化 Python代码(sklearn库) 什么是线性回归(Linear regression) 引例 假设某地区租房价格只与房屋 ...
- linux查看端口被那个进程占用
linux下遇到端口被暂用了 需要知道是哪个进程 比如80端口 可以这样 netstat -tunlp|