传送门

显然树上第k大直接主席树

如果连边的话,我们重构小的那一棵,连到另一棵上。

说起来简单,调了我一晚上。

总的来说3个错误:

1.离散化写错位置写到了后面

2."="写成了"=="

3.加双向边时加成了单向边

3个错误3个小时。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 80011
#define M 10000000 using namespace std; int n, m, T, cnt, tot, test, last;
int head[N], to[N << 2], nex[N << 2], val[N], ntr[N], deep[N], f[N][21], root[N], sum[M], ls[M], rs[M], fa[N], size[N];
bool vis[N]; inline int read()
{
int x = 0, f = 1;
char ch = getchar();
for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -1;
for(; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + ch - '0';
return x * f;
} inline void add(int x, int y)
{
to[cnt] = y;
nex[cnt] = head[x];
head[x] = cnt++;
} inline int find(int x)
{
return x == fa[x] ? x : fa[x] = find(fa[x]);
} inline void Union(int x, int y)
{
int fx = find(x), fy = find(y);
if(fx != fy) fa[fx] = fy, size[fy] += size[fx];
} inline int query(int a, int b, int c, int d, int l, int r, int x)
{
if(l == r) return l;
int mid = (l + r) >> 1;
if(sum[ls[a]] + sum[ls[b]] - sum[ls[c]] - sum[ls[d]] >= x) return query(ls[a], ls[b], ls[c], ls[d], l, mid, x);
else return query(rs[a], rs[b], rs[c], rs[d], mid + 1, r, x - (sum[ls[a]] + sum[ls[b]] - sum[ls[c]] - sum[ls[d]]));
} inline void insert(int &now, int last, int l, int r, int x)
{
now = ++tot;
ls[now] = ls[last];
rs[now] = rs[last];
sum[now] = sum[last] + 1;
if(l == r) return;
int mid = (l + r) >> 1;
if(x <= mid) insert(ls[now], ls[last], l, mid, x);
else insert(rs[now], rs[last], mid + 1, r, x);
} inline void dfs(int u)
{
int i, v;
vis[u] = 1;
deep[u] = deep[f[u][0]] + 1;
insert(root[u], root[f[u][0]], 1, m, val[u]);
for(i = 0; f[u][i]; i++) f[u][i + 1] = f[f[u][i]][i];
for(; i <= 20; i++) f[u][i] = 0;
for(i = head[u]; ~i; i = nex[i])
{
v = to[i];
if(!vis[v])
{
f[v][0] = u;
dfs(v);
}
}
vis[u] = 0;
} inline int lca(int x, int y)
{
int i;
if(deep[x] < deep[y]) swap(x, y);
for(i = 20; i >= 0; i--)
if(deep[f[x][i]] >= deep[y]) x = f[x][i];
if(x == y) return x;
for(i = 20; i >= 0; i--)
if(f[x][i] != f[y][i]) x = f[x][i], y = f[y][i];
return f[x][0];
} int main()
{
char s[10];
int i, x, y, k, fx, fy;
test = read();
n = read();
m = read();
T = read();
memset(head, -1, sizeof(head));
for(i = 1; i <= n; i++)
{
fa[i] = i, size[i] = 1;
val[i] = ntr[i] = read();
}
for(i = 1; i <= m; i++)
{
x = read();
y = read();
add(x, y);
add(y, x);
Union(x, y);
}
sort(ntr + 1, ntr + n + 1);
m = unique(ntr + 1, ntr + n + 1) - ntr - 1;
for(i = 1; i <= n; i++) val[i] = lower_bound(ntr + 1, ntr + m + 1, val[i]) - ntr;
for(i = 1; i <= n; i++)
if(!deep[i]) dfs(i);
while(T--)
{
scanf("%s", s);
x = read() ^ last;
y = read() ^ last;
if(s[0] == 'Q')
{
k = read() ^ last;
printf("%d\n", last = ntr[query(root[x], root[y], root[lca(x, y)], root[f[lca(x, y)][0]], 1, m, k)]);
}
else
{
fx = find(x), fy = find(y);
if(size[fx] > size[fy]) swap(x, y);
Union(x, y);
f[x][0] = y;
dfs(x);
add(x, y);
add(y, x);
}
}
return 0;
}

  

[luoguP3302] [SDOI2013]森林(主席树 + 启发式合并 + lca)的更多相关文章

  1. luoguP3302 [SDOI2013]森林 主席树 启发式合并

    题目链接 luoguP3302 [SDOI2013]森林 题解 本来这题树上主席树暴力启发式合并就完了 结果把lca写错了... 以后再也不这么写了 复杂度\(O(nlog^2n)\) "f ...

  2. Bzoj 3123: [Sdoi2013]森林(主席树+启发式合并)

    3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当前 ...

  3. [bzoj3123] [SDOI2013]森林 主席树+启发式合并+LCT

    Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...

  4. BZOJ 3123: [Sdoi2013]森林 [主席树启发式合并]

    3123: [Sdoi2013]森林 题意:一个森林,加边,询问路径上k小值.保证任意时刻是森林 LCT没法搞,树上kth肯定要用树上主席树 加边?启发式合并就好了,小的树dfs重建一下 注意 测试点 ...

  5. [BZOJ3123][Sdoi2013]森林 主席树+启发式合并

    3123: [Sdoi2013]森林 Time Limit: 20 Sec  Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当 ...

  6. [SDOI2013]森林 主席树+启发式合并

    这题的想法真的很妙啊. 看到题的第一眼,我先想到树链剖分,并把\(DFS\)序当成一段区间上主席树.但是会发现在询问的时候,可能会非常复杂,因为你需要把路径拆成很多条轻链和重链,它们还不一定连续,很难 ...

  7. 【BZOJ 3123】 [Sdoi2013]森林 主席树启发式合并

    我们直接按父子关系建主席树,然后记录倍增方便以后求LCA,同时用并查集维护根节点,而且还要记录根节点对应的size,用来对其启发式合并,然后每当我们合并的时候我们都要暴力拆小的一部分重复以上部分,总时 ...

  8. 【BZOJ-3123】森林 主席树 + 启发式合并

    3123: [Sdoi2013]森林 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 2738  Solved: 806[Submit][Status] ...

  9. bzoj 3123 [Sdoi2013]森林(主席树+启发式合并+LCA)

    Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...

随机推荐

  1. [VC]在VC++中实现让程序只运行一个实例的方法且实现该实例

    方法一: 有时候在开发应用程序时,希望控制程序运行唯一的实例.例如,最常用的mp3播放软 件Winamp,由于它需要独占计算机中的音频设备,因此该程序只允许自身运行唯一的一个例程.在Visual C+ ...

  2. 撤销git pull命令

    比如:在master分支上执行了git pull命令,想回到pull之前分支所在的commit位置. 步骤一:用 git reflog master 查看master分支的历史变动记录,其中有一个就是 ...

  3. python 基础之while无限循环

    用户登录程序 username = "chenxi" passwed = "testki" counter = 0 while counter < 3: ...

  4. leetcode 105. Construct Binary Tree from Preorder and Inorder Traversal,剑指offer 6 重建二叉树

    不用迭代器的代码 class Solution { public: TreeNode* reConstructBinaryTree(vector<int> pre,vector<in ...

  5. python_82_标准库_random模块

    import random print(help(random.random)) #随机整数 print(random.randint(1,7))#生成一个[1, 7]的随机整数 print(rand ...

  6. @Param注解在dao层的使用

    有时在前台用ajax传过来许多参数,不知道在mybatis如何封装,就要用到@Param注解了,这时就不需要在映射文件写传入参数了,这种方法虽然比较取巧,但还是很实用的,如下图:

  7. HTML5<header>元素

    HTML5<header>元素 1.header元素描述了文档的头部区域,主要用于定义内容的介绍展示区域. 2.实例: <header> <h2>heder元素描述 ...

  8. iOS 后台传输服务

    后台传输服务 — 我们用水壶来比喻 (0:14) 后天传输服务是 iOS 7 引进的 API,它准许应用暂停或者中止之后,在后台继续执行网络服务(比如下载或者上传).举个例子,这正是 Dropbox ...

  9. CPL学习笔记(二)

    数组 数组(array) 声明通用格式为: typeName arrayName[arraySize]. arrayName指定元素数目,必须为整型常量或const值.(不能是变量) int a[3] ...

  10. 基于Centos7.2搭建Cobbler自动化批量部署操作系统服务

    1       Cobbler服务器端系统环境配置 1.1     系统基本环境准备 [root@cobbler-server ~]# cat /etc/redhat-release CentOS L ...