【线段树分治】【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 < ...
随机推荐
- [Hanani]高数相关知识记录
分部积分 \(\int uv'{\rm d}x=uv-\int u'v{\rm d}x\)
- Redux 入门教程(一):基本用法
转自http://www.ruanyifeng.com/blog/2016/09/redux_tutorial_part_one_basic_usages.html(仅供个人学习使用) 首先明确一点, ...
- 图片人脸检测(OpenCV版)
图片人脸检测 人脸检测使用到的技术是OpenCV,上一节已经介绍了OpenCV的环境安装,点击查看. 功能展示 识别一种图上的所有人的脸,并且标出人脸的位置,画出人眼以及嘴的位置,展示效果图如下: 多 ...
- 基础业务:滚动到指定位置导航固定(CSS实现)
最近公司做的业务都是使用Vue.Element写的,涉及到的相应的基础业务像轮播.预加载.懒加载,都是使用 NPM上的工具来实现,原理和基础还是要有的,就来实现几个项目中常用到的业务. 经常见到这样的 ...
- SQL 查一年内的数据
--查询今年的 select * from 表 where datediff(yy,时间字段,GETDATE())=0 --查询去年的 select * from 表 where datediff(y ...
- Beta版发布 - 感谢有你们
在本次Beta版开发的过程中,很感谢组长王航对我的信任,让我统筹大家的工作任务和进度,使我对项目管理有了深刻的理解. 我也要感谢邹双黛,在beta版开发中因为邹双黛的帮助,我对于文字类工作已经越来越得 ...
- 软工1816 · Beta冲刺(6/7)
团队信息 队名:爸爸饿了 组长博客:here 作业博客:here 组员情况 组员1(组长):王彬 过去两天完成了哪些任务 推进Web端完成开发 推进修改一些后端接口的逻辑 着手制作视频 接下来的计划 ...
- “吃神么,买神么”的第一个Sprint计划(第七天)
“吃神么,买神么”项目Sprint计划 ——5.25 星期一(第五天)立会内容与进度 摘要: 所有的部件都完成,在贴每个人负责的部件时发现很多问题,很多网页布局的运用不熟练,一部分的div会跑位置~ ...
- PHP面试题一
http://www.viphper.com/?p=28 1.用PHP打印出前一天的时间格式是2006-5-10 22:21:21(2分) $a = date("Y-m-d H:i:s&qu ...
- Swift-KVC构造函数中数据类型和私有属性