BZOJ3514: Codechef MARCH14 GERALD07加强版【LCT】【主席树】【思维】
Description
N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数。
Input
第一行四个整数N、M、K、type,代表点数、边数、询问数以及询问是否加密。
接下来M行,代表图中的每条边。
接下来K行,每行两个整数L、R代表一组询问。对于type=0的测试点,读入的L和R即为询问的L、R;对于type=1的测试点,每组询问的L、R应为L xor lastans和R xor lastans。
Output
K行每行一个整数代表该组询问的联通块个数。
Sample Input
3 5 4 0
1 3
1 2
2 1
3 2
2 2
2 3
1 5
5 5
1 2
Sample Output
2
1
3
1
HINT
对于100%的数据,1≤N、M、K≤200,000。
2016.2.26提高时限至60s
思路
考虑怎么解决联通块个数
首先我们如果从左向右扫所有的边加进联通块中
可以用LCT维护一个动态的生成树
然后一旦出现环就删除环中编号最小的一个边
然后我们考虑一下左端点如果大于当前删除的边的编号,实际上还是存在这个联通块的
所以说我们只需要统计在\([l,r]\)这个区间中左端点小于\(l\)的个数就可以知道有多少个联通块了
然后当当前节点存在自环的时候就会发现他的删除的边的编号就是自己,不然会有多余的贡献
然后主席树维护就没了
#include<bits/stdc++.h>
using namespace std;
const int N = 4e5 + 10;
int n, m, k, typ;
int u[N], v[N];
int lastans = 0;
namespace Link_Cut_Tree {
int ch[N][2], fa[N], minval[N], val[N], rev[N], cnt = 0;
bool isroot(int t) {
return ch[fa[t]][0] != t && ch[fa[t]][1] != t;
}
void pushup(int t) {
minval[t] = val[t];
if (ch[t][0]) minval[t] = min(minval[t], minval[ch[t][0]]);
if (ch[t][1]) minval[t] = min(minval[t], minval[ch[t][1]]);
}
void pushnow(int t) {
swap(ch[t][0], ch[t][1]);
rev[t] ^= 1;
}
void pushdown(int t) {
if (!isroot(t)) pushdown(fa[t]);
if (rev[t]) {
pushnow(ch[t][0]);
pushnow(ch[t][1]);
rev[t] = 0;
}
}
void newnode(int vl) {
++cnt;
fa[cnt] = ch[cnt][0] = ch[cnt][1] = 0;
minval[cnt] = val[cnt] = vl;
rev[cnt] = 0;
}
bool son(int t) {
return t == ch[fa[t]][1];
}
void rotate(int t) {
int f = fa[t], g = fa[f];
bool a = son(t), b = a ^ 1;
if (!isroot(f)) ch[g][son(f)] = t;
fa[t] = g;
ch[f][a] = ch[t][b];
fa[ch[t][b]] = f;
ch[t][b] = f;
fa[f] = t;
pushup(f);
pushup(t);
}
void splay(int t) {
pushdown(t);
while (!isroot(t)) {
int f = fa[t];
if (!isroot(f)) {
if (son(f) ^ son(t)) rotate(t);
else rotate(f);
}
rotate(t);
}
}
void access(int t) {
int tmp = 0; // 需要设定初值
while (t) {
splay(t);
ch[t][1] = tmp;
pushup(t);
tmp = t;
t = fa[t];
}
}
void makeroot(int t) {
access(t);
splay(t);
pushnow(t);
}
void link(int x, int y) {
makeroot(x);
fa[x] = y;
}
void cut(int x, int y) {
makeroot(x);
access(y);
splay(y);
fa[x] = ch[y][0] = 0;
pushup(y);
}
};
using Link_Cut_Tree::minval;
using Link_Cut_Tree::newnode;
using Link_Cut_Tree::link;
using Link_Cut_Tree::cut;
using Link_Cut_Tree::makeroot;
using Link_Cut_Tree::access;
namespace Functional_Segment_Tree {
const int LOG = 30;
int cnt = 0;
int siz[N * LOG], ls[N * LOG], rs[N * LOG], rt[N];
void insert(int &t, int last, int l, int r, int pos) {
t = ++cnt;
ls[t] = ls[last];
rs[t] = rs[last];
siz[t] = siz[last] + 1;
if (l == r) return;
int mid = (l + r) >> 1;
if (pos <= mid) insert(ls[t], ls[last], l, mid, pos);
else insert(rs[t], rs[last], mid + 1, r, pos);
}
int query(int t, int last, int l, int r, int ql, int qr) {
if (ql <= l && r <= qr) return siz[t] - siz[last];
int mid = (l + r) >> 1;
if (qr <= mid) return query(ls[t], ls[last], l, mid, ql, qr);
else if (ql > mid) return query(rs[t], rs[last], mid + 1, r, ql, qr);
else return query(ls[t], ls[last], l, mid, ql, mid) + query(rs[t], rs[last], mid + 1, r, mid + 1, qr);
}
};
using Functional_Segment_Tree::rt;
using Functional_Segment_Tree::insert;
using Functional_Segment_Tree::query;
namespace Union_Find {
int fa[N << 1];
void init() {
for (int i = 1; i <= n; i++)
fa[i] = i;
}
int find(int x) {
return x == fa[x] ? x : fa[x] = find(fa[x]);
}
bool merge(int x, int y) {
int fax = find(x), fay = find(y);
if (fax == fay) return 0;
fa[fax] = fay;
return 1;
}
}
using Union_Find::init;
using Union_Find::merge;
int main() {
#ifdef dream_maker
freopen("input.txt", "r", stdin);
#endif
scanf("%d %d %d %d", &n, &m, &k, &typ);
for (int i = 1; i <= n; i++) newnode(m + 1);
for (int i = 1; i <= m; i++) newnode(i);
init();
for (int i = 1; i <= m; i++) {
scanf("%d %d", &u[i], &v[i]);
if (u[i] == v[i]) {
insert(rt[i], rt[i - 1], 0, m, i);
continue;
}
if (merge(u[i], v[i])) {
insert(rt[i], rt[i - 1], 0, m, 0);
link(u[i], n + i);
link(v[i], n + i);
} else {
makeroot(u[i]);
access(v[i]);
makeroot(v[i]);
int cur = minval[v[i]];
cut(v[cur], n + cur);
cut(u[cur], n + cur);
link(v[i], n + i);
link(u[i], n + i);
insert(rt[i], rt[i - 1], 0, m, cur);
}
}
for (int i = 1; i <= k; i++) {
int l, r; scanf("%d %d", &l, &r);
if (typ) l ^= lastans, r ^= lastans;
if (l > r) swap(l, r);
lastans = n - query(rt[r], rt[l - 1], 0, m, 0, l - 1);
printf("%d\n", lastans);
}
return 0;
}
BZOJ3514: Codechef MARCH14 GERALD07加强版【LCT】【主席树】【思维】的更多相关文章
- [BZOJ3514]CodeChef MARCH14 GERALD07加强版(LCT+主席树)
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 2177 Solved: 834 ...
- BZOJ 3514: Codechef MARCH14 GERALD07加强版 [LCT 主席树 kruskal]
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1312 Solved: 501 ...
- BZOJ 3514: Codechef MARCH14 GERALD07加强版( LCT + 主席树 )
从左到右加边, 假如+的边e形成环, 那么记下这个环上最早加入的边_e, 当且仅当询问区间的左端点> _e加入的时间, e对答案有贡献(脑补一下). 然后一开始是N个连通块, 假如有x条边有贡献 ...
- 【BZOJ3514】Codechef MARCH14 GERALD07加强版 LCT+主席树
题解: 还是比较简单的 首先我们的思路是 确定起点 然后之后贪心的选择边(也就是越靠前越希望选) 我们发现我们只需要将起点从后向前枚举 然后用lct维护连通性 因为强制在线,所以用主席树记录状态就可以 ...
- bzoj3514 Codechef MARCH14 GERALD07加强版 lct预处理+主席树
Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1951 Solved: 746[Submi ...
- 【BZOJ-3514】Codechef MARCH14 GERALD07加强版 LinkCutTree + 主席树
3514: Codechef MARCH14 GERALD07加强版 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1288 Solved: 490 ...
- BZOJ3514: Codechef MARCH14 GERALD07加强版(LCT,主席树)
Description N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数. Input 第一行四个整数N.M.K.type,代表点数.边数.询问数以及询问是否加密.接下来M ...
- BZOJ3514 Codechef MARCH14 GERALD07加强版 LCT维护最大生成树 主席树
题面 考虑没有询问,直接给你一个图问联通块怎么做. 并查集是吧. 现在想要动态地做,那么应该要用LCT. 考虑新加进来一条边,想要让它能够减少一个联通块的条件就是现在边的两个端点还没有联通. 如果联通 ...
- BZOJ3514 Codechef MARCH14 GERALD07加强版 LCT+可持久化线段树
自己独自想出来并切掉还是很开心的~ Code: #include <bits/stdc++.h> #define N 400005 #define inf 1000000000 #defi ...
- BZOJ3514 Codechef MARCH14 GERALD07加强版 LCT
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3514 题意概括 N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数. N ...
随机推荐
- AES SBox的构造(python)
几点需要注意的,求解逆元的时候使用的是拓展欧几里得,但是那些运算规则需要变一变,模2的加减乘除(或者可以理解为多项式的运算) 在进行字节的仿射变换不用进行矩阵的运算. 一个矩阵和一个列向量进行运算的时 ...
- java日期操作 大全
先来一个: 取得指定月份的第一天与取得指定月份的最后一天 http://iamin.blogdriver.com/iamin/847990.html )); } ...
- Java IO流-File类
2017-10-24 23:50:22 File类 File类:文件和目录路径名的抽象表示形式.该文件或者目录不一定真实存在. * 构造方法 File类有四种构造方法,主要采用的构造方法师第二种,也就 ...
- php json josn_decode()返回的是对像,如何把对像转成数组
php json josn_decode()返回的是对像,如何把对像转成数组 a.php传值页面,使用 json_encode($array)对数组进行加密码. b.php页面在接收a.php传过来的 ...
- 封装DLL并调用
c# DLL封装并调用 1.封装自己的dll: a.打开visual studio - 文件 - 新建 - 项目- 类库 - 名称MyTestDll: b.右键Class1.cs - 修改为 Te ...
- 20170503xlVBA房地产数据分类连接
Sub NextSeven_CodeFrame4() Application.ScreenUpdating = False Application.DisplayAlerts = False Appl ...
- 12月14日 bs-grid , destroy_all()
bootstrap Grid : The Bs grid system has four classes: xs (phones), sm (tablets), md (desktops), and ...
- H5基础知识(一)
一.概述 HTML5 是html4.0 升级版 结构 Html5 .样式 css3 .行为: API 都有所增强 HTML5并不仅仅只是做为HTML标记语言的一个最新版本,更重要的是它制定了Web ...
- spting Boot 创建一个springBoot项目
spting Boot 创建一个springBoot项目 1)学习springBoot使用软件:IDEA软件(前面的文章有安装idea的过程). 也可以使用另一种方法在https://start.sp ...
- ES批量索引写入时的ID自动生成算法
对bulk request的处理流程: 1.遍历所有的request,对其做一些加工,主要包括:获取routing(如果mapping里有的话).指定的timestamp(如果没有带timestamp ...