WC模拟（1.6） T2 送你一棵圣诞树

最新推荐文章于 2020-08-16 19:07:49 发布
最新推荐文章于 2020-08-16 19:07:49 发布
阅读量372 收藏
点赞数
分类专栏： NOIP解题报告 线段树 树状数组 LCA 文章标签： NOI 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/scar_lyw/article/details/79035054
版权

送你一棵圣诞树

题目背景：

1.6 WC模拟T2

分析：树状数组 + 线段树 + lca

 

考场上没有想到nlog2n的做法······只想到nlog3n，然后自己并不想写，然后写了应该也是T到飞起······显然因为是子树操作，很简单的想法就是直接求出dfs序，然后转化为序列查询，然后我们发现，我们有两种可行的方法，第一种用主席树，外层颜色，内层区间，然后区间上在加上权值线段树，对于每一个i记录上一个和它颜色相同的点pre[i]，那么对于一个询问就是求某一段颜色中，某一段区间中pre[i] < dfn[u]的个数。然后发现，它还需要进行修改，那么最外层还有一个树状数组，然后华丽丽的三个log，还有一种方法是分块，考虑首先按照区间分块，块内按照pre排序，然后对于每一个块针对颜色区间建主席树，然后对于一个块的询问，直接二分找到pre = l - 1的位置，在主席树里面直接查询颜色个数就可以了，（块中只记录每一个颜色的第一个出现位置），而对于修改，直接暴力重建所在块就可以了，这个复杂度是O(n * sqrt(n) * logn)，这两种方法都是基于你只想到了在序列上处理，而忽略了本身一棵树的性质，我们先把相同的颜色的位置的dfs序提出来排序，然后在相同颜色的相邻dfs序的lca处减一，这样就可以直接变成序列求和了，主席树就可以解决了，因为带修改，那么变成树状数组套线段树就可以了，复杂度O(nlog2n)

 

Source:

/*
	created by scarlyw
*/
#include <cstdio>
#include <string>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
#include <cctype>
#include <vector>
#include <set>
#include <queue>
#include <ctime>
#include <bitset>

inline char read() {
	static const int IN_LEN = 1024 * 1024;
	static char buf[IN_LEN], *s, *t;
	if (s == t) {
		t = (s = buf) + fread(buf, 1, IN_LEN, stdin);
		if (s == t) return -1;
	}
	return *s++;
}

///*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = read(), iosig = false; !isdigit(c); c = read()) {
		if (c == -1) return ;
		if (c == '-') iosig = true;	
	}
	for (x = 0; isdigit(c); c = read()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int OUT_LEN = 1024 * 1024;
char obuf[OUT_LEN], *oh = obuf;
inline void write_char(char c) {
	if (oh == obuf + OUT_LEN) fwrite(obuf, 1, OUT_LEN, stdout), oh = obuf;
	*oh++ = c;
}

template<class T>
inline void W(T x) {
	static int buf[30], cnt;
	if (x == 0) write_char('0');
	else {
		if (x < 0) write_char('-'), x = -x;
		for (cnt = 0; x; x /= 10) buf[++cnt] = x % 10 + 48;
		while (cnt) write_char(buf[cnt--]);
	}
}

inline void flush() {
	fwrite(obuf, 1, oh - obuf, stdout);
}

/*
template<class T>
inline void R(T &x) {
	static char c;
	static bool iosig;
	for (c = getchar(), iosig = false; !isdigit(c); c = getchar())
		if (c == '-') iosig = true;	
	for (x = 0; isdigit(c); c = getchar()) 
		x = ((x << 2) + x << 1) + (c ^ '0');
	if (iosig) x = -x;
}
//*/

const int MAXN = 100000 + 10;

int cnt, n, q, t, x, y, ind;
std::set<int> col[MAXN];
std::vector<int> edge[MAXN];
int c[MAXN], pos[MAXN], root[MAXN], in[MAXN], out[MAXN], size[MAXN], dep[MAXN];
int father[MAXN], son[MAXN], top[MAXN], last[MAXN];
std::set<int>::iterator it, it1, it2;

struct node {
	int left, right, sum;
} tree[MAXN * 500];

inline void insert(int &cur, int l, int r, int pos, int x) {
	if (cur == 0) cur = ++cnt;
	tree[cur].sum += x;
	if (l == r) return ;
	int mid = l + r >> 1;
	if (pos <= mid) insert(tree[cur].left, l, mid, pos, x);
	else insert(tree[cur].right, mid + 1, r, pos, x);
}

inline int query(int cur, int l, int r, int ql, int qr) {
	if (ql <= l && r <= qr) return tree[cur].sum;
	int mid = l + r >> 1, ret = 0;
	if (ql <= mid) ret += query(tree[cur].left, l, mid, ql, qr);
	if (qr > mid) ret += query(tree[cur].right, mid + 1, r, ql, qr);
	return ret;
}

inline void insert(int i, int c, int x) {
	for (; i <= n; i += (i & -i)) insert(root[i], 1, n, c, x);
}

inline void add_edge(int x, int y) {
	edge[x].push_back(y), edge[y].push_back(x);
}

inline void dfs1(int cur, int fa) {
	in[cur] = ++ind, size[cur] = 1, dep[cur] = dep[fa] + 1, father[cur] = fa;
	for (int p = 0; p < edge[cur].size(); ++p) {
		int v = edge[cur][p];
		if (v != fa) {
			dfs1(v, cur), size[cur] += size[v];
			if (size[v] > size[son[cur]]) son[cur] = v;
		}
	}
	out[cur] = ind, pos[in[cur]] = cur;
}

inline void dfs2(int cur, int tp) {
	top[cur] = tp;
	if (son[cur]) dfs2(son[cur], tp);
	for (int p = 0; p < edge[cur].size(); ++p) {
		int v = edge[cur][p];
		if (top[v] == 0) dfs2(v, v);
	}
}

inline int query_lca(int u, int v) {
	while (top[u] != top[v]) 
		(dep[top[u]] > dep[top[v]]) ? u = father[top[u]] : v = father[top[v]];
	return (dep[u] > dep[v]) ? v : u;
}

inline void solve_tree() {
	R(n), R(q), R(t);
	for (int i = 1; i <= n; ++i) R(c[i]);
	for (int i = 1; i < n; ++i) R(x), R(y), add_edge(x, y);
	dfs1(1, 0), dfs2(1, 1);
}

inline solve_struct() {
	for (int i = 1; i <= n; ++i) col[i].insert(0);
	for (int i = 1; i <= n; ++i) col[c[pos[i]]].insert(i);
	for (int i = 1; i <= n; ++i) {
		x = c[pos[i]], insert(i, x, 1);
		if (last[x] != 0) insert(in[query_lca(last[x], pos[i])], x, -1);
		last[x] = pos[i];
	}
}

inline void solve_modify(int cur, int c) {
	int pos = in[cur], x = ::c[cur];
	if (c == x) return ;
	it = it1 = it2 = col[x].lower_bound(pos), it1--, it2++;
	if (*it1 != 0) insert(in[query_lca(::pos[*it1], cur)], x, 1);
	if (it2 != col[x].end()) insert(in[query_lca(::pos[*it2], cur)], x, 1);
	if (*it1 != 0 && it2 != col[x].end()) 
		insert(in[query_lca(::pos[*it1], ::pos[*it2])], x, -1);
	insert(pos, x, -1), col[x].erase(it), x = ::c[cur] = c, insert(pos, x, 1);
	col[x].insert(pos), it1 = it2 = col[x].lower_bound(pos), it1--, it2++;
	if (*it1 != 0) insert(in[query_lca(::pos[*it1], cur)], x, -1);
	if (it2 != col[x].end()) insert(in[query_lca(::pos[*it2], cur)], x, -1);
	if (*it1 != 0 && it2 != col[x].end()) 
		insert(in[query_lca(::pos[*it1], ::pos[*it2])], x, 1);
}

inline int solve_query(int l, int r, int ql, int qr) {
	int ans = 0;
	for (int i = r; i; i -= (i & -i)) ans += query(root[i], 1, n, ql, qr);
	for (int i = l - 1; i; i -= (i & -i)) ans -= query(root[i], 1, n, ql, qr);
	return ans;
}

inline void solve() {
	static int type, u, l, r, c, ans = 0;
	solve_tree(), solve_struct();
	while (q--) {
		R(type);
		if (type == 1) {
			R(u), R(l), R(r), (t ? (u ^= ans, l ^= ans, r ^= ans) : 0);
			W(ans = solve_query(in[u], out[u], l, r)), write_char('\n');
		} else {
			R(u), R(c), (t ? (u ^= ans, c ^= ans) : 0);
			solve_modify(u, c);
		}
	}
}

int main() {
	freopen("xmastree1.in", "r", stdin);
	freopen("xmastree1.out", "w", stdout);
	solve();
	flush();
	return 0;
}


scarlyw
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
送你一颗圣诞树
欢迎您的光顾awa
08-17 549
1484.送你一颗圣诞树， dp + map记忆化。
[JZOJ5512] 送你一棵圣诞树
Never give in.
01-09 2026
Description 一棵 n 个点的树, 树根为 1. 一开始每个点上有一个 1…n 的颜色 ci, 不同点颜色可以相同. 现在有 q 次操作, 分为两种类型: • 1 u l r: 询问子树 u 中有多少种在 l 到 r 之间的颜色至少出现了一次 • 2 u c: 将 u 的颜色修改为 c 要求强制在线 n,q Solution 这题先不要转到DFS序列上做
WC模拟（1.6） T3 送你一棵圣诞树
scarlyw的博客
01-11 337
状压DP + 子集枚举
【清华冬令营2018模拟】送你一棵圣诞树
Jacky35
01-08 1237
Description 送你一棵n 个点的树, 树根为1. 一开始每个点上有一个1:::n 的颜色ci, 不同点颜色可以相同. 现在有q 次操作, 分为两种类型: • 1 u l r: 询问子树u 中有多少种在l 到r 之间的颜色至少出现了一次 • 2 u c: 将u 的颜色修改为c 部分测试点要求强制在线. Input 第一行三个整数n; q; t, 分别表示树的点数, 操作的个数
jzoj5512 送你一棵圣诞树 树状数组套线段树+set
olahiuj的博客
01-10 390
Description 一棵 n 个点的树, 树根为 1. 一开始每个点上有一个 1…n 的颜色 ci, 不同点颜色可以相同. 现在有 q 次操作, 分为两种类型: • 1 u l r: 询问子树 u 中有多少种在 l 到 r 之间的颜色至少出现了一次 • 2 u c: 将 u 的颜色修改为 c 对于前 20% 的数据, n, q ≤ 5000; 对于前 40% 的数据,
一天一个shell命令 linux文本操作系列-wc命令详解
01-11
wc是一个用于统计的工具，全拼为：Word Count(单词统计) 实例： 1. 统计行数 $wc –l file 2. 统计单词数 $wc –w file 3. 统计字符数 $wc –c file 4. 统计流中的字符数 $echo –n 1234 | wc –c 输出：4 解读–...
论文研究 - HFCVD方法在WC-Co工具上金刚石膜生长中衬底温度分布的模拟
05-26
基材温度是通过热丝化学气相沉积（HFCVD）制造的金刚石膜... 结果表明，钻头的表面覆盖有一层连续的金刚石膜，这证明了模拟沉积参数是有益的，并提供了一种新方法来调节用于沉积金刚石膜的CVD反应器中的基板温度分布。
wc.rar_压缩解压_小波 压缩_小波零树
09-23
基于小波零树量化的压缩解压缩程序，算法比较的多。
wc.rar_wc
09-22
在linux9中实现统计文件的功能，并且统计屏幕输入。在redhat9中已经调试通过
洛谷P3939填颜色
Liuzhou_lyt的博客
01-20 530
题目链接: 点击打开链接 看到题第一眼，线段树?树状数组?主席树(蒟蒻不会)全冒出来? 可是仔细看题发现这是道水题 我们先把每种颜色读入 预处理出每种颜色所在位置的排序 拿样例举例: 1 2 3 2 3 3 我们第一次初始化为b数组： 第一种颜色下标有:1 第二种颜色下标有:2,4 第三种颜色下标有:3,5,6 故b数组为1,2,4,3,5,6 用另
【JZOJ5512】送你一棵圣诞树
sadnohappy
01-11 645
Description一棵 n 个点的树, 树根为 1. 一开始每个点上有一个 1…n 的颜色 ci, 不同点颜色可以相同. 现在有 q 次操作, 分为两种类型: - 1 u l r: 询问子树 u 中有多少种在 l 到 r 之间的颜色至少出现了一次 - 2 u c: 将 u 的颜色修改为 c 强制在线
第二类斯特林数【清华冬令营2018模拟】送你一个DAG
Marcus
01-08 1067
题目 1.1 Description 送你一个n 个点m 条边的DAG 和参数k, 定义一条经过l 条边的路径的权值为l^k. 对于i=1…n, 求出所有1 到i 的路径的权值之和, 对998244353 取模. 1.2 Input Format 第一行三个整数n; m; k, 分别表示DAG 的点数, 边数和参数. 接下来m 行, 每行两个整数ui; vi, 表示一条从ui 到vi ...
【WC2014】时空穿梭（莫比乌斯反演）（组合数学）
FSYo 的博客
06-05 268
传送门 考虑枚举各维最大最小坐标的差量 Δi\Delta_iΔi​，可以写出式子：Ans=∑Δ(gcd(Δ1...n)−1c−2)∏i=1n(mi−Δi)=∑d(d−1c−2)∑Δ∏i=1n(mi−dΔi)[gcd(Δ1...n)=1]=∑d(d−1c−2)∑lμ(l)∑Δ(∏i=1nmi−dlΔi)=∑T∏i=1n(∑jmi/Tmi−Tj)∑l∣Tμ(l)(d−1c−2)Ans=\sum_{\Delta}\binom{gcd(\Delta_{1...n})-1}{c-2}\prod_{i=1}^n (m
送你一颗圣诞树【DP】【DFS】
//(^ o ^)//
08-16 454
>Description 再过三个多月就是圣诞节了，小R 想送小Y 一棵圣诞树作为节日礼物。因为他想让这棵圣诞树越大越好，所以当然是买不到能够让他满意的树的，因此他打算自己把这棵树拼出来。 现在，小R 开始画这棵树的设计图纸了。因为这棵树实在太大，所以他采用了一种比较方便的方法。首先他定义了m+ 1 棵树T0 到Tm。最开始他只画好了T0 的图纸：就只有一个点，编号为0。 接着，对于每一棵树Ti，他在第Tai 棵树的第ci 个点和第Tbi 棵树的第di 个点之间连上了一条长度为li 的边。在Ti 中，
[JZOJ5513] 送你一棵圣诞树2
Never give in.
01-09 488
Description 给出一棵N个点的树 在树上进行随机点分治，每个节点需要你确定一个不超过K的权值，点分树上的父亲的权值严格小于儿子的。 问有多少种不同的结果，两种结果不同当且仅当某一个节点在点分树上的父亲不同或者它们的权值不同 1 Solution 直接按照题意模拟，N基本上都要存在于指数中 考虑转化题意。 点分治是有着严格的大小关系的，不妨倒过来想，确定
NOIP2017 D2T3 列队
scarlyw的博客
12-19 1194
平衡树 or 线段树 or 树状数组
NOI模拟（5.11） BJOID2T2 链上二次求和 （bzoj5291）
scarlyw的博客
05-12 530
线段树
2017.10.6 BJOI2015 bzoj4336 骑士的旅行
scarlyw的博客
10-17 480
树链剖分 + 树状数组 + 权值线段树 + 二分
NOIP模拟（11.07）T3 小店购物
scarlyw的博客
11-07 405
线段树 + multiset
你的方法是不可行的，换一个
最新发布
06-03
如果您想切换到另一个应用程序的窗口，也可以使用 `dumpsys` 命令来获取当前所有的应用程序窗口，并根据窗口顺序来切换到下一个或上一个窗口。以下是一个示例脚本： ```bash #!/bin/bash # 获取所有应用程序窗口的信息 windows=$(adb shell dumpsys window windows | grep -E 'mCurrentFocus|mFocusedApp') # 获取当前窗口的包名和类名 current_package=$(echo "$windows" | grep mCurrentFocus | cut -d'/' -f1 | cut -d' ' -f7) current_class=$(echo "$windows" | grep mCurrentFocus | cut -d'/' -f2 | cut -d'}' -f1) # 获取当前应用程序的包名 current_app=$(echo "$windows" | grep mFocusedApp | cut -d'=' -f2 | cut -d' ' -f1 | cut -d'/' -f1) # 获取所有应用程序的包名和类名 packages=$(echo "$windows" | grep 'mBasePackageName' | cut -d'=' -f2 | cut -d' ' -f1) classes=$(echo "$windows" | grep 'mBasePackageName' | cut -d'=' -f2 | cut -d' ' -f2 | cut -d'}' -f1) # 将包名和类名合并为完整的组件名 components=$(paste -d/ <(echo "$packages") <(echo "$classes")) # 获取当前应用程序在所有应用程序中的位置 index=$(echo "$components" | grep -n "^$current_package/$current_class$" | cut -d':' -f1) # 计算下一个或上一个应用程序的位置 next_index=$((($index + 1) % $(echo "$components" | wc -l))) prev_index=$((($index - 1 + $(echo "$components" | wc -l)) % $(echo "$components" | wc -l))) # 获取下一个或上一个应用程序的组件名 next_component=$(echo "$components" | sed -n "${next_index}p") prev_component=$(echo "$components" | sed -n "${prev_index}p") # 切换到下一个或上一个应用程序 if [ "$1" == "next" ]; then adb shell am start -n "$next_component" elif [ "$1" == "prev" ]; then adb shell am start -n "$prev_component" fi ``` 以上脚本可以获取当前所有的应用程序窗口，并根据窗口顺序来计算下一个或上一个应用程序的位置和组件名，最后使用 `am start` 命令来启动下一个或上一个应用程序的窗口。例如，您可以使用以下命令来切换到下一个应用程序： ``` ./switch_app.sh next ``` 您也可以使用以下命令来切换到上一个应用程序： ``` ./switch_app.sh prev ``` 请注意，此脚本的执行需要您的设备已经连接到 adb，并且您需要有足够的权限来启动应用程序。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值