[题解] 洛谷 P3603 雪辉

最新推荐文章于 2024-10-29 21:24:58 发布
最新推荐文章于 2024-10-29 21:24:58 发布
阅读量109 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/yydyz/p/10599491.html
版权

模拟赛中遇到了这个题，当时我这个沙雕因为把一个\(y\)打成了\(x\)而爆零。回来重新写这道题，莫名其妙的拿了rank1。。。
我的解法与其他几位的题解有些不同我太蒻了。并没有选取所谓的关键点，而是用树链剖分将树上问题转化为序列的问题。

序列问题的解决

  • 现在的问题是，维护一个序列，要求查询若干个区间的并集的权值种类和\(mex\)，无修改操作。
  • 如果只有一个区间，那么这个问题就可以很愉快的用主席树解决了。
  • 然而要查询的是区间并的权值种类和\(mex\)，而且这两个东西都不满足区间可加性，于是可以放弃线段树的思路了。
  • 可以想到用分块解决这个问题，设每一块大小为\(B\)
  • 考虑到点权\(\leq 30000\)，又要合并若干个区间，于是可以想到用\(bitset\)维护区间权值信息。
  • 设\(f[][]\)为一个\(bitset\)的二维数组，\(f[l][r]\)表示从第\(l\)块到第\(r\)块(包含\(l,r\))之间的权值集合，这个数组可以在\(O(n\frac{30000}{32})\)的时间内预处理出来。每次查询的时候，整块的信息可以在\(O(\frac{30000}{32})\)的时间内解决。
  • 剩余的部分直接暴力处理，时间复杂度为\(O(B)\)。这样单次查询可以做到\(O(\frac{30000}{32}+B)\)
  • 至于区间并，直接开一个全局变量记录答案(我的代码中用的是\(cur\))，依次处理每个区间。因为重叠的部分并不会对答案有影响。现在序列的问题已经解决了。

序列处理部分的代码

int b[N]; // 记录区间上每一个位置属于哪个块
int L[N], R[N]; // 每个块的左右断点
bitset<W> cur; // 查询时用到的全局变量
bitset<W> f[110][110]; // 预处理的 f 数组

void preWork() {
  // 预处理
  for (int i = 1; i <= n; ++i) {
    // 先计算单个块的权值情况
    b[i] = (i-1) / B + 1;
    f[b[i]][b[i]].set(a[i]);
  }
  // 处理每个块的左右端点
  for (int i = 1; i <= b[n]; ++i)
    L[i] = R[i-1] + 1, R[i] = i * B;
  R[b[n]] = n; // 最后一个块的右端点要特判
  for (int i = 1; i < b[n]; ++i)
    for (int j = i+1; j <= b[n]; ++j) // 计算 f 数组
      f[i][j] = f[i][j-1] | f[j][j];
}

void queryOnBlock(int l, int r) {
  if (b[l] == b[r]) {
    // 特判左右端点在同一个块内的情况
    for (int i = l; i <= r; ++i) cur.set(a[i]);
    return;
  }
  cur |= f[b[l]+1][b[r]-1]; // 两块之间的部分直接查询
  for (int i = l; i <= R[b[l]]; ++i) cur.set(a[i]); // 左边的剩余部分
  for (int i = L[b[r]]; i <= r; ++i) cur.set(a[i]); // 右边的剩余部分
}

int mex(bitset<W> &s) {
  // 暴力求 mex
  for (int i = 0; i < W; ++i)
    if (!s.test(i)) return i;
  return 1e9;
}

将树上问题转化为序列问题

树链剖分的板子(我这个蒟蒻写挂的部分)。。。
直接贴代码了

int G[N], ed = 1, w[N]; // 树的存储
struct Edge {
  int to, nxt;
  Edge() { to = nxt = 0; }
  Edge(int to, int nxt) : to(to), nxt(nxt) {}
} e[N<<1];

inline void addEdge(int x, int y) {
  e[++ed] = Edge(y, G[x]), G[x] = ed;
  e[++ed] = Edge(x, G[y]), G[y] = ed;
}

// 树链剖分相关
int dfn[N]; // dfs 序
int fa[N]; // 父结点
int son[N]; // 重儿子
int top[N]; // 重链顶端
int size[N]; // 子数大小
int dep[N]; // 深度
int a[N]; // 转化的序列

void dfs1(int x, int p) {
  size[x] = 1, fa[x] = p;
  for (int i = G[x]; i != 0; i = e[i].nxt) {
    int y = e[i].to;
    if (y == p) continue;
    dep[y] = dep[x] + 1;
    dfs1(y, x);
    size[x] += size[y];
    if (size[son[x]] < size[y])
      son[x] = y;
  }
}

void dfs2(int x, int t) {
  static int cur = 0;
  dfn[x] = ++cur, a[cur] = w[x], top[x] = t;
  if (!son[x]) return;
  dfs2(son[x], t);
  for (int i = G[x]; i != 0; i = e[i].nxt) {
    int y = e[i].to;
    if (y == son[x] || y == fa[x]) continue;
    dfs2(y, y);
  }
}

void queryOnTree(int x, int y) {
  // 树上查询
  while (top[x] != top[y]) {
    if (dep[top[x]] < dep[top[y]]) swap(x, y);
    queryOnBlock(dfn[top[x]], dfn[x]);
    x = fa[top[x]];
  }
  if (dfn[x] > dfn[y]) swap(x, y);
  queryOnBlock(dfn[x], dfn[y]);
}

总代码

// 2598ms 57.76MB 无O2
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <bitset>
#define LL long long
using namespace std;

inline int getint() {
  int x = 0, f = 1; char c = getchar();
  while (!isdigit(c)) { if (c == '-') f = 0; c = getchar(); }
  while (isdigit(c)) { x = (x*10) + (c-'0'); c = getchar(); }
  return f ? x : -x;
}

const int N = 1e5 + 10;
const int B = 1e3;
const int W = 30001;

int n, m, flag; // flag 判断是否需要异或
int ans1, ans2, lastans; // ans1 为权值种类，ans2 为权值 mex

int G[N], ed = 1, w[N]; // 树的存储
struct Edge {
  int to, nxt;
  Edge() { to = nxt = 0; }
  Edge(int to, int nxt) : to(to), nxt(nxt) {}
} e[N<<1];

inline void addEdge(int x, int y) {
  e[++ed] = Edge(y, G[x]), G[x] = ed;
  e[++ed] = Edge(x, G[y]), G[y] = ed;
}

// 树链剖分相关
int dfn[N]; // dfs 序
int fa[N]; // 父结点
int son[N]; // 重儿子
int top[N]; // 重链顶端
int size[N]; // 子数大小
int dep[N]; // 深度
int a[N]; // 转化的序列

void dfs1(int x, int p) {
  size[x] = 1, fa[x] = p;
  for (int i = G[x]; i != 0; i = e[i].nxt) {
    int y = e[i].to;
    if (y == p) continue;
    dep[y] = dep[x] + 1;
    dfs1(y, x);
    size[x] += size[y];
    if (size[son[x]] < size[y])
      son[x] = y;
  }
}

void dfs2(int x, int t) {
  static int cur = 0;
  dfn[x] = ++cur, a[cur] = w[x], top[x] = t;
  if (!son[x]) return;
  dfs2(son[x], t);
  for (int i = G[x]; i != 0; i = e[i].nxt) {
    int y = e[i].to;
    if (y == son[x] || y == fa[x]) continue;
    dfs2(y, y);
  }
}

int b[N]; // 记录区间上每一个位置属于哪个块
int L[N], R[N]; // 每个块的左右断点
bitset<W> cur; // 查询时用到的全局变量
bitset<W> f[110][110]; // 预处理的 f 数组

void preWork() {
  // 预处理
  for (int i = 1; i <= n; ++i) {
    // 先计算单个块的权值情况
    b[i] = (i-1) / B + 1;
    f[b[i]][b[i]].set(a[i]);
  }
  // 处理每个块的左右端点
  for (int i = 1; i <= b[n]; ++i)
    L[i] = R[i-1] + 1, R[i] = i * B;
  R[b[n]] = n; // 最后一个块的右端点要特判
  for (int i = 1; i < b[n]; ++i)
    for (int j = i+1; j <= b[n]; ++j) // 计算 f 数组
      f[i][j] = f[i][j-1] | f[j][j];
}

void queryOnBlock(int l, int r) {
  if (b[l] == b[r]) {
    // 特判左右端点在同一个块内的情况
    for (int i = l; i <= r; ++i) cur.set(a[i]);
    return;
  }
  cur |= f[b[l]+1][b[r]-1]; // 两块之间的部分直接查询
  for (int i = l; i <= R[b[l]]; ++i) cur.set(a[i]); // 左边的剩余部分
  for (int i = L[b[r]]; i <= r; ++i) cur.set(a[i]); // 右边的剩余部分
}

void queryOnTree(int x, int y) {
  // 树上查询
  while (top[x] != top[y]) {
    if (dep[top[x]] < dep[top[y]]) swap(x, y);
    queryOnBlock(dfn[top[x]], dfn[x]);
    x = fa[top[x]];
  }
  if (dfn[x] > dfn[y]) swap(x, y);
  queryOnBlock(dfn[x], dfn[y]);
}

int mex(bitset<W> &s) {
  // 暴力求 mex
  for (int i = 0; i < W; ++i)
    if (!s.test(i)) return i;
  return 1e9;
}

int main() {
  cin >> n >> m >> flag;
  for (int i = 1; i <= n; ++i) w[i] = getint();
  for (int i = 1; i < n; ++i) {
    int x = getint(), y = getint();
    addEdge(x, y);
  }
  dfs1(1, 0), dfs2(1, 1);
  preWork();
  for (int i = 1; i <= m; ++i) {
    cur.reset();
    int num = getint();
    for (int j = 1; j <= num; ++j) {
      int x = getint() ^ (flag*lastans);
      int y = getint() ^ (flag*lastans);
      queryOnTree(x, y);
    }
    ans1 = cur.count(), ans2 = mex(cur);
    lastans = ans1 + ans2;
    printf("%d %d\n", ans1, ans2);
  }
  return 0;
}

转载于:https://www.cnblogs.com/yydyz/p/10599491.html

weixin_30586257
关注
题解 洛谷 P1763 埃及分数
qq_14886903的博客
02-14 547
这么经典的题目居然这么少人写题解，那么我这个蒟蒻就来写一发 Java 题解罢！ 由于这道题的搜索树深度无限，所以排除掉深搜；每层能扩展的广度也无限，所以排除广搜。那么，我们就只能使用一种深搜广搜结合的算法——迭代加深搜索。迭代加深搜索就是不断放宽搜索树最大深度的深搜。 我们还需要一些剪枝优化。 优化1：为了避免重复搜索，我们强制让分母递增。 优化2：由于分母递增，那么如果当前搜索的分数过小，在...
[BZOJ4763][P3603]雪辉[手写bitset+静态分块]
weixin_33943347的博客
12-28 122
题意：给一个n个点的树，点有点权，有m次询问，每次询问多条链的并有多少种不同的点权以及它的mex mex就是一个集合中最小的没有出现的非负整数，注意0要算 rand出 \(\sqrt n\)个点，把每次查询拆成 x->fx0->fx->lca->fy->fy0->y #include <bits/stdc++.h> using namespace s...
P6177 Count on a tree II
orz
02-04 1442
给你一棵 nnn 个节点且带点权的树，mmm 个询问，每次查询链颜色数。 强制在线。 1≤n≤4×104,1≤m≤1051\leq n \leq 4\times 10^4,1\leq m\leq 10^51≤n≤4×104,1≤m≤105。 sol 首先如果不强制在线，用树上莫队即可。 但多了个强制在线，容易想到是预处理题。 查询链颜色数，比较好的一种方法是使用 bitset，对值域建 bitset，答案就是 bitset 中 111 的数量。 那么现在的问题就是怎么把一条路径上的 bitset 并起来。
[洛谷]CON1466 洛谷2017春节联欢赛 Hello Dingyou题解 & Bzoj4763雪辉
weixin_30815427的博客
03-14 278
题目来源：https://www.luogu.org/contest/show?tid=1466 创建时间：2017/3/13 18:33镇楼图： 猜猜她是谁~ 解题思路： 春节居然也有比赛，一脸萌比完全不知道啊。。其实我本来逛着Bzoj，发现未来日记出新题了——[bzoj4763雪辉]，要不shy也来乱搞一下>_<结果不会哦，那就来洛谷一道一道切题。做...
洛谷P1443 马的遍历 题解 bfs
weixin_51189382的博客
02-05 175
方法：bfs 先放AC代码： #include<iostream> #include<queue> #include<algorithm> using namespace std; typedef pair<int, int> PII; #define x first #define y second int dx[] = { -2,-1,1,2,2,1,-1,-2 }, dy[] = { 1,2,2,1,-1,-2,-2,-1 }; const int N
洛谷P3613 睡觉困难综合征
蒟蒻的博客
03-10 371
题目背景刚立完Flag我就挂了WC和THUWC。。。时间限制0.5s，空间限制128MB因为Claris大佬帮助一周目由乃通过了Deus的题，所以一周目的由乃前往二周目世界找雪辉去了由于二周目世界被破坏殆尽，所以由乃和雪辉天天都忙着重建世界（其实和MC差不多吧），Deus看到了题问她，总是被告知无可奉告Deus没办法只能去三周目世界问三周目的由乃OI题。。。三周目的世界中，因为没有未来日记，所以一...
洛谷P3352 [ZJOI2016]线段树
东方幻想乡真棒（^W^）
06-05 250
P3352 [ZJOI2016]线段树 (^ w ^) 题目描述 小Yuuka遇到了一个题目：有一个序列a_1,a_2,?,a_n，q次操作，每次把一个区间内的数改成区间内的最大值，问最后每个数是多少。小Yuuka很快地就使用了线段树解决了这个问题。 于是充满智慧的小Yuuka想，如果操作是随机的，即在这q次操作中每次等概率随机地选择一个区间l,r，然后将这个区间内的数改成区间内最大值（注意这样的...
（差分）洛谷P4231 三步必杀
东方幻想乡真棒（^W^）
03-31 338
洛谷P4231 三步必杀 题目背景 （三）旧都 离开狭窄的洞穴，眼前豁然开朗。 天空飘着不寻常的雪花。 一反之前的幽闭，现在面对的，是繁华的街市，可以听见酒碗碰撞的声音。 这是由被人们厌恶的鬼族和其他妖怪们组成的小社会，一片其乐融融的景象。 诶，不远处突然出现了一些密密麻麻的小点，好像大颗粒扬尘一样。 离得近了点，终于看清楚了。 长着角的鬼们聚在一起，围观着另一只鬼的表演。 那”扬尘”，其实都是弹...
（树状数组+逆元）洛谷P5142 区间方差
东方幻想乡真棒（^W^）
04-01 654
洛谷P5142 区间方差 (^ w ^) 题目背景 出题人并没有能力写有趣的题面…… 题目描述 对于一个长度为n的序列a1,a2,a3⋯ana_1,a_2,a_3\cdots a_na1​,a2​,a3​⋯an​，我们定义它的平均数a为: a=1n∑i=1naia=\frac{1}{n}\sum_{i=1}^{n}a_ia=n 1​∑i=1n​ai​ 并定义它的方差d为: d=1n∑i=1n(ai...
洛谷P1724题解（请不要抄袭！）
06-17
洛谷P1427题解，若要复制or下载，请与我联系，如要改编，也请联系，若发现未经我同意的复制、下载、改编，直接举报。
洛谷题解：P1421-小玉买文具
04-16
《洛谷题解：P1421-小玉买文具》是一道典型的编程入门题，主要考察编程者对于C++顺序结构编程的理解和应用。题目描述的是一个简单的数学问题，即小玉要用不超过某金额去购买文具，每件文具的价格是19分，问最多可以...
洛谷题解：P1001-A+B Problem
04-16
洛谷题解：P1001-A+B Problem主要涉及到知识点包括C++程序设计基础、程序输入输出操作、C++标准库的使用、位运算和C++编译预处理指令。 1. C++程序设计基础： - C++是编译型语言，需要编译成机器语言后才能运行。 ...
题解 洛谷P1093 奖学金
mabaoguochuanren的博客
12-08 261
题解 P1093 奖学金 本蒟蒻的题解 毕竟是2007年NOIP的题，所以看似很难 仔细一想，其实却很简单。 蒟蒻方法：结构体； 贴代码 1.创建结构体 #include<bits/stdc++.h>//只会这一个 using namespace std; struct zky{ int xh,y,s,w,zf; }a[305];//创建结构体 2.确定排序方法 bool cmp(const zky&nd1,const zky& nd2) { if(nd1.zf
HJ51输出单向链表中倒数第k个结点
最新发布
晴山の博客
10-29 300
HJ51输出单向链表中倒数第k个结点。
【LeetCode】每日一题 2024_10_29 生成不含相邻零的二进制字符串（DFS）
Locky136的博客
10-29 304
【LeetCode】每日一题 2024_10_29 生成不含相邻零的二进制字符串（DFS）
算法练习：四数之和
2302_80652761的博客
10-28 326
四数之和
【动态规划】力扣198.打家劫舍
J_pluto的博客
10-27 392
题目中指出不可以选相邻的房间，说明如果选了第 1 间，那么第 2 间一定不可以选，第 3 间房间可以选，也可以不选……假设是按照从第 1 间房间开始依次往后选，每一间房间的选择都会影响到后面的房间可不可以被选择。在递归的基础上，创建一个数组 dp[] 来存储之前计算过的结果，避免重复计算。但这样提交会发现超时。
【LeetCode】每日一题 2024_10_27 冗余连接（并查集）
Locky136的博客
10-27 337
【LeetCode】每日一题 2024_10_27 冗余连接（并查集）
Leetcode 合并 K 个升序链表
985小菜鸡
10-29 366
题目要求合并 ( K ) 个已经升序排列的链表，输出一个合并后的有序链表。代码使用**优先队列（PriorityQueue）**来实现，以确保每次取出当前最小的节点，这样可以高效地合并多个链表。这段代码通过优先队列实现了多链表的高效合并，每次从优先队列中取出最小节点并连接到结果链表上，从而保证了结果链表的有序性。这样就可以在相对较低的时间复杂度下完成 ( K ) 个有序链表的合并。
洛谷P1699C++题解
07-05
洛谷(P1699)是一道C++题目，通常这类编号表示的是该题目的难度和ID。然而，由于我无法直接查看实时的在线题库，我需要简述一下解题思路的一般框架。 P1699可能是指的一个算法题目，它可能涉及到数据结构、数学建模...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值