【线段树分治】【P5227】 [AHOI2013]连通图
Description
给定一个无向连通图和若干个小集合,每个小集合包含一些边,对于每个集合,你需要确定将集合中的边删掉后改图是否保持联通。集合间的询问相互独立
定义一个图为联通的当且仅当对于任意的两个顶点,都存在一条路径连接它们
Input
第一行为两个整数 \(n,m\),代表无向图的点数和边数
下面 \(m\) 行,包含两个整数 \(u,v\),代表该边连接点 \(u,v\)。第 \(i + 1\) 行的边的编号为 \(i\)。保证不存在重边和自环
下面一行包含一个整数 \(k\),表示集合个数
下面 \(k\) 行每行描述一个集合,每行的第一个数为集合中边的个数 \(c\),后面 \(c\) 个数代表集合内的边
Output
对于每个集合,输出一行代表去掉该集合中的边后图是否联通,如果联通输出 Connected
,否则输出 Disconnected
Hint
\(1~\leq~n,k~\leq~10^5\)
\(1~\leq~m~\leq~2~\times~10^5\)
\(1~\leq~c~\leq~4\)
Solution
这个分治方法叫做线段树分治,而不是cdq分治。
我们考虑分治,设 \([l,r]\) 为区间 \([l,r]\) 中询问涉及边都没有被连上且其它边都连上时图是否联通,那么递归到底时即为一次询问的答案。
具体方法为先将没有被询问过的边连上。然后在递归区间 \([l,r]\) 时,将区间内 \([l,mid]\) 的询问中不涉及的边全部连上,递归处理左侧区间,然后撤销上次的连边,再将 \((mid, r]\) 中询问不涉及的边都连上,处理右边,回溯后再次撤销。在这个过程中用并查集维护一下联通块个数,在底部即可判断。存储需要撤销的操作可以用一个栈实现。
考虑时间复杂度:递归层数为 \(O(\log k)\),于是每个询问被处理 \(O(\log k)\) 次,对每个询问中的每条边,每被处理一次就会有一次并查集的删除和撤销操作,由于是按秩合并实现,单次操作复杂度 \(O(\log n)\)。所以单个询问的复杂度为 \(O(c~\log n)\)。共有 \(k\) 个操作,于是总的时间复杂度为 \(O(k~\times~c~\times~\log k~\log n)\)。
这个分治和cdq分治的区别大概在于……此分治时 mid 两侧的区间是被 “平等对待” 的,即记录右边区间对左边区间的贡献,但是cdq分治不是。
另外上面做线段树分治的时候并没有将树显式的建出来,但他们在本质上是一样的。
另外还有一种神仙做法在这里
Code
#include <cstdio>
#include <vector>
#ifdef ONLINE_JUDGE
#define freopen(a, b, c)
#endif
typedef long long int ll;
namespace IPT {
const int L = 1000000;
char buf[L], *front=buf, *end=buf;
char GetChar() {
if (front == end) {
end = buf + fread(front = buf, 1, L, stdin);
if (front == end) return -1;
}
return *(front++);
}
}
template <typename T>
inline void qr(T &x) {
char ch = IPT::GetChar(), lst = ' ';
while ((ch > '9') || (ch < '0')) lst = ch, ch=IPT::GetChar();
while ((ch >= '0') && (ch <= '9')) x = (x << 1) + (x << 3) + (ch ^ 48), ch = IPT::GetChar();
if (lst == '-') x = -x;
}
namespace OPT {
char buf[120];
}
template <typename T>
inline void qw(T x, const char aft, const bool pt) {
if (x < 0) {x = -x, putchar('-');}
int top=0;
do {OPT::buf[++top] = static_cast<char>(x % 10 + '0');} while (x /= 10);
while (top) putchar(OPT::buf[top--]);
if (pt) putchar(aft);
}
const int maxn = 100010;
const int maxm = 200010;
struct Edge {
int from, to;
};
Edge edge[maxm];
struct Ask {
int c;
int mu[5];
};
Ask ask[maxn];
struct OP {
int a, b, p;
};
OP stack[maxm * 20];
int n, m, k, cnt, top;
int ufs[maxn], rk[maxn];
bool vis[maxm];
void reset();
int find(int);
void cont(int, bool);
void unionn(int, int, bool);
void work(int, int, std::vector<int>&);
void divide(int, int, std::vector<int>&);
int main() {
freopen("1.in", "r", stdin);
qr(n); qr(m); cnt = n;
for (int i = 1; i <= m; ++i) {
qr(edge[i].from); qr(edge[i].to);
}
qr(k);
for (int i = 1; i <= k; ++i) {
qr(ask[i].c);
for (int j = 1; j <= ask[i].c; ++j) {qr(ask[i].mu[j]); vis[ask[i].mu[j]] = true;}
}
for (int i = 1; i <= n; ++i) ufs[i] = i, rk[i] = 1;
std::vector<int>dn;
for (int i = 1; i <= m; ++i)
if (!vis[i]) cont(i, false);
else dn.push_back(i);
work(1, k, dn);
return 0;
}
void work(int l, int r, std::vector<int>&e) {
for (int i = l; i <= r; ++i) {
for (int j = 1; j <= ask[i].c; ++j) vis[ask[i].mu[j]] = true;
}
std::vector<int>dn;
int _tp = top;
for (int i = 0, len = e.size(); i < len; ++i) {
if (!vis[i]) cont(i, true);
else dn.push_back(i);
}
for (int i = l; i <= r; ++i) {
for (int j = 1; j <= ask[i].c; ++j) vis[ask[i].mu[j]] = false;
}
divide(l, r, dn);
while (top != _tp) reset();
}
void divide(int l, int r, std::vector<int>&e) {
if (l == r) {
puts(cnt == 1 ? "Connected" : "Disconnected");
return;
}
int mid = (l + r) >> 1;
work(l, mid, e);
work(mid + 1, r, e);
}
int find(int x) {return ufs[x] == x ? x : find(ufs[x]);}
void cont(int x, bool rec) {
int fa = find(edge[x].from), fb = find(edge[x].to);
if (fa != fb) {
unionn(fa, fb, rec);
--cnt;
}
}
void unionn(int a, int b, bool rec) {
int c;
if (rk[a] < rk[b]) {c = 1; ufs[a] = b;}
else if (rk[a] > rk[b]) {c = 2; ufs[b] = a;}
else {c = 3; ufs[b] = a; ++rk[a];}
if (rec) stack[++top] = {a, b, c};
}
void reset() {
int a = stack[top].a, b = stack[top].b, p = stack[top--].p;
switch (p) {
case 1:
ufs[a] = a; break;
case 2:
ufs[b] = b; break;
case 3:
ufs[b] = b; --rk[a]; break;
}
++cnt;
}
【线段树分治】【P5227】 [AHOI2013]连通图的更多相关文章
- 线段树分治初步学习&洛谷P5227[AHOI2013]连通图
线段树分治 其实思想说起来是比较简单的,我们把这个题里的所有操作(比如连边删边查询balabala)全部拍到一棵线段树上,然后对着整棵树dfs一下求解答案,顺便把操作做一下,回溯的时候撤销一下即可.虽 ...
- BZOJ3237:[AHOI2013]连通图(线段树分治,并查集)
Description Input Output Sample Input 4 5 1 2 2 3 3 4 4 1 2 4 3 1 5 2 2 3 2 1 2 Sample Output Connec ...
- BZOJ3237 AHOI2013连通图(线段树分治+并查集)
把查询看做是在一条时间轴上.那么每条边都有几段存在时间.于是线段树分治就好了. 然而在bzoj上t掉了,不知道是常数大了还是写挂了. 以及brk不知道是啥做数组名过不了编译. #include< ...
- 3237: [Ahoi2013]连通图 线段树分治
题解: cf765f cf671e bzoj4184 bzoj4552 线段树分治裸题 还是介绍一下线段树分治 这个东西其实挺简单但也挺有用的 可以把删除+插入操作变成只有插入(倒着就是删除) 像这一 ...
- 线段树分治总结(线段树分治,线段树,并查集,树的dfn序,二分图染色)
闲话 stO猫锟学长,满脑子神仙DS 网上有不少Dalao把线段树分治也归入CDQ分治? 还是听听YCB巨佬的介绍: 狭义:只计算左边对右边的贡献. 广义:只计算外部对内部的贡献. 看来可以理解为广义 ...
- 【线段树分治 线性基】luoguP3733 [HAOI2017]八纵八横
不知道为什么bzoj没有HAOI2017 题目描述 Anihc国有n个城市,这n个城市从1~n编号,1号城市为首都.城市间初始时有m条高速公路,每条高速公路都有一个非负整数的经济影响因子,每条高速公路 ...
- LOJ2312 LUOGU-P3733「HAOI2017」八纵八横 (异或线性基、生成树、线段树分治)
八纵八横 题目描述 Anihc国有n个城市,这n个城市从1~n编号,1号城市为首都.城市间初始时有m条高速公路,每条高速公路都有一个非负整数的经济影响因子,每条高速公路的两端都是城市(可能两端是同一个 ...
- loj#2312. 「HAOI2017」八纵八横(线性基 线段树分治)
题意 题目链接 Sol 线性基+线段树分治板子题.. 调起来有点自闭.. #include<bits/stdc++.h> #define fi first #define se secon ...
- BZOJ.4184.shallot(线段树分治 线性基)
BZOJ 裸的线段树分治+线性基,就是跑的巨慢_(:з」∠)_ . 不知道他们都写的什么=-= //41652kb 11920ms #include <map> #include < ...
随机推荐
- JavaScript学习(2)call&apply&bind&eval用法
javascript学习(2)call&apply&bind&eval用法 在javascript中存在这样几种特别有用的函数,能方便我们实现各种奇技淫巧.其中,call.bi ...
- redis高级应用(集群搭建、集群分区原理、集群操作)
文章主目录 Redis集群简介 Redis集群搭建 Redis集群分区原理 集群操作 参考文档 本文是redis学习系列的第四篇,前面我们学习了redis的数据结构和一些高级特性,点击下面链接可回看 ...
- MVC 漫长之路(一)
1.新建项目 mvc 视图引擎选中 Razor 2.允许我们设置这个项目关于 MVC 的一些设置,确认选中了“空”项目模板 3. 4.打开 Global.asax 文件 配置路由名称等 ...
- 做一名合格的程序员(learning of a previous team)
合格程序员的基本素质 作一个真正合格的程序员,或者说就是可以真正合格完成一些代码工作的程序员,应该具 有的素质. 1 团队精神和协作能力 把它作为基本素质,并不是不重要,恰恰相反,这是程序员应该具备 ...
- 迎来OO的曙光,总结规格的意义——OO第四次博客总结
一切都要结束了,砥砺前行~ 一.测试与正确性论证的效果差异 测试,顾名思义就是我们暴力用大量数据轰炸编写的程序的过程.日常的OO过程中,我们经常互相寻求“测试集”,正是因为测试使用特定数据对我们的功能 ...
- 用java进行简单的万年历编写
import java.util.Scanner; public class PrintCalendarDemo1 { public static void main(String[] args) { ...
- 第三次作业— C++计算器项目的初始部分
作业题目: C++计算器项目的初始部分 仓库 代码: Scan.h #ifndef SCAN_H #define SCAN_H #include<string> #include<i ...
- Git的基本使用方法和安装&心得体会(使用git命令行)
这是补发的,使用命令行操作的. (1)选择本地repository的路径 找到后点鼠标右键,选择git bash here. (2) clone到本地 在命令行输入 git clone ADDRESS ...
- Unity3D游戏开发——收集当前关卡游戏中分散的物件
运用场景 许多游戏中会有一些供玩家拾起的物件,例如装备.血包.道具等.当玩家与物件进行碰撞后,则会进入仓库. 本篇介绍了简单的碰撞过程. 原理 基本的碰撞机制,用到OnTriggerEnter()碰撞 ...
- 读书笔记-《Linux/Unix设计思想》
本书主要讲的是Unix程序设计思想,具体涉及到linux的内容不多. 整本书的一个基本出发点是开源.其中主要强调的观点包括: 1.小即是美 作者持有的主要观点是程序应该以小为美.小程序实现小功能,每个 ...