P4899 [IOI2018] werewolf 狼人(kruskal 重构树 + 主席树)

P4899 [IOI2018] werewolf 狼人

给定一个有nnn个点mmm条边的无向图,有QQQ个询问

每次输入S,E,L,RS, E, L, RS,E,L,R,表示你在SSS点出发,要到EEE点,且初始时你是人形态,你只能走[L,n][L, n][L,n]的点,

但是我们要以狼的形态走到终点,这个时候只能走[1,R][1, R][1,R]之间的点,

也就是说我们要在[L,R][L, R][L,R]的某个点变身,问能否从S−>ES->ES>E

一个比较容易想到的做法:

显然,初始状态我们要满足S∈[L,n],E∈[1,R]S \in [L, n], E \in[1, R]S[L,n],E[1,R],否则就没有讨论的意义了。

假设满足S∈[L,n],E∈[1,R]S \in[L, n], E \in[1,R]S[L,n],E[1,R]

我们对每条边按照{u,v,w=min(u,v)}\{u, v, w = min(u, v)\}{u,v,w=min(u,v)}这样的方式,跑一个降序kruskalkruskalkruskal树,然后从SSS点出发,向上跳,找到使点权不小于LLL,能到达的点uuu,

同样的我们对每条边按照{u,v,w=max(u,v}\{u, v, w = max(u, v\}{u,v,w=max(u,v}的方式,跑一个升序kruskalkruskalkruskal树,然后从EEE出发,向上跳,找到使点权不大于RRR,能到达的点vvv

接下来我们只需判断uuu所代表的点集SuS_uSuvvv所代表的点集SvS_vSv,判断这两个点集有没有交集即可。

设图一进栈dfsdfsdfs序为L1L1L1,出栈dfsdfsdfs序为R1R1R1,同理图二为L2,R2L2, R2L2,R2

x[i]x[i]x[i]为在图一dfsdfsdfs序为iii的点编号,y[i]y[i]y[i]为,点编号为iii的点在图二的dfsdfsdfs序,接下来我们只要在区间[L1u,R1u][L1_u, R1_u][L1u,R1u]的主席树上查找区间[L2v,R2v][L2_v, R2_v][L2v,R2v]是否有值即可。

#include <bits/stdc++.h>using namespace std;const int N = 4e5 + 10;struct Res {int u, v, w;bool operator < (const Res &t) const {return w < t.w;}
}edge[N];int n, m, Q, x[N], y[N];int root[N], ls[N << 5], rs[N<< 5], sum[N << 5], num;void update(int &rt, int pre, int l, int r, int x, int v) {rt = ++num;ls[rt] = ls[pre], rs[rt] = rs[pre], sum[rt] = sum[pre] + v;if (l == r) {return ;}int mid = l + r >> 1;if (x <= mid) {update(ls[rt], ls[pre], l, mid, x, v);}else {update(rs[rt], rs[pre], mid + 1, r, x, v);}
}int query(int rt1, int rt2, int l, int r, int L, int R) {// cout << l << " " << r << "\n";if (l >= L && r <= R) {return sum[rt2] - sum[rt1];}int ans = 0, mid = l + r >> 1;if (L <= mid) {ans += query(ls[rt1], ls[rt2], l, mid, L, R);}if (R > mid) {ans += query(rs[rt1], rs[rt2], mid + 1, r, L, R);}return ans;
}struct Kruskal {int head[N], to[N], nex[N], cnt = 1;int ff[N], value[N], fa[N][21], l[N], r[N], rk[N], tot, nn;void add(int x, int y) {to[cnt] = y;nex[cnt] = head[x];head[x] = cnt++;}int find(int rt) {return ff[rt] == rt ? rt : ff[rt] = find(ff[rt]);}void dfs(int rt, int f) {fa[rt][0] = f, l[rt] = ++tot, rk[tot] = rt;for (int i = 1; i <= 20; i++) {fa[rt][i] = fa[fa[rt][i - 1]][i - 1];}for (int i = head[rt]; i; i = nex[i]) {if (to[i] == f) {continue;}dfs(to[i], rt);}r[rt] = tot;}void kruskal(int rev) {sort(edge + 1, edge + 1 + m);if (rev) {reverse(edge + 1, edge + 1 + m);}for (int i = 1; i < N; i++) {ff[i] = i;}nn = n;for (int i = 1, cur = 1; i <= m && cur < n; i++) {int u = find(edge[i].u), v = find(edge[i].v);if (u != v) {cur++, nn++;ff[u] = ff[v] = nn;value[nn] = edge[i].w;add(nn, u), add(nn, v);if (u <= n) {value[u] = edge[i].w;}if (v <= n) {value[v] = edge[i].w;}}}dfs(nn, 0);}bool judge(int u, int x, int flag) {return flag ? value[u] >= x : value[u] <= x;}int work(int u, int x, int flag) {for (int i = 20; i >= 0; i--) {if (fa[u][i] && judge(fa[u][i], x, flag)) {u = fa[u][i];}}return u;}
}a, b;int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);scanf("%d %d %d", &n, &m, &Q);for (int i = 1; i <= m; i++) {scanf("%d %d", &edge[i].u, &edge[i].v);edge[i].u++, edge[i].v++;}for (int i = 1; i <= m; i++) {edge[i].w = min(edge[i].v, edge[i].u);}a.kruskal(1);for (int i = 1; i <= m; i++) {edge[i].w = max(edge[i].v, edge[i].u);}b.kruskal(0);for (int i = 1; i < 2 * n; i++) {x[i] = a.rk[i];y[i] = b.l[i];}for (int i = 1; i < 2 * n; i++) {root[i] = root[i - 1];if (x[i] <= n) {update(root[i], root[i - 1], 1, 2 * n, y[x[i]], 1);}}while (Q--) {int u, v, L, R;scanf("%d %d %d %d", &u, &v, &L, &R);u++, v++, L++, R++;u = a.work(u, L, 1), v = b.work(v, R, 0);int l1 = a.l[u], r1 = a.r[u], l2 = b.l[v], r2 = b.r[v];printf("%d\n", query(root[l1 - 1], root[r1], 1, 2 * n, l2, r2) != 0);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/313844.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

微软+开源,那些亲爱的以及热爱的

微软 Reactor 社区空间开幕式暨 Azure Meetup 社区活动已于9月7日在上海圆满结束&#xff01;但是…如何构建一个可持续发展的社区未来的路&#xff0c;仍然很长...你应该知道的微软 Reactor微软 Reactor 是微软为构建开发者社区而提供的一个社区空间&#xff0c;以“予力多元…

Java线程问题问答

1、多线程的作用&#xff1f; 发挥多核CPU的优势&#xff0c;提高效率防止阻塞便于建模&#xff1a;将一个任务拆分成多个子任务&#xff0c;分别建立程序模型 2&#xff1a;平时项目中使用锁和synchronized比较多&#xff0c;而很少使用volatile&#xff0c;难道就没有保证可见…

ASP.NET Core 2.2 项目升级至 3.0 备忘录

.NET Core 3.0及ASP.NET Core 3.0 前瞻ASP.NET Core 3.0 迁移避坑指南将 ASP.NET Core 2.2 迁移至 ASP.NET Core 3.0 需要注意的地方记录在这篇随笔中。TargetFramework 改为 netcoreapp3.0 <TargetFramework>netcoreapp3.0</TargetFramework>从 Web 项目&#xff…

Java 死锁

目录&#xff1a; 什么是死锁&#xff1f;死锁是怎么产生的&#xff1f;怎么排查死锁&#xff1f;死锁的预防拓展&#xff1a;Java CPU 100%排查 一 什么是死锁&#xff1f; 注&#xff1a;线程和进程都可能会产生死锁&#xff0c;以下以线程为例 死锁是指两个或两个以上的…

CF196E Opening Portals

CF196E Opening Portals 给定一个有nnn个节点&#xff0c;mmm条边的无向联通图&#xff0c;有kkk个点有portalsportalsportals&#xff0c;当经过了某个点&#xff0c;如果这个点有portalportalportal&#xff0c;它就会永久开启&#xff0c; 对于任意两个开启的portalportal…

.NET Core 微信小程序退款——(统一退款)

点击上方“dotNET名人堂”&#xff0c;选择“设为星标”用学习的姿态&#xff0c;步入工作的状态继上一篇".NET Core 微信小程序支付——&#xff08;统一下单&#xff09;后"&#xff0c;本文将实现统一退款功能&#xff0c;能支付就应该能退款嘛&#xff0c;一般涉…

Java 锁之 CAS

什么是CAS(compare and swap)&#xff1f; CAS&#xff08;Compare & Set&#xff0c;或是 Compare & Swap&#xff09;&#xff0c;即比较并交换&#xff0c;也是实现我们平时所说的自旋锁或乐观锁的核心操作。 它的实现很简单&#xff0c;就是用一个预期的值和内存…

牛客练习赛71 F 红蓝图(kruskal重构树)

红蓝图 给定两个参数x,tx, tx,t&#xff0c;删除边权大于ttt的红边&#xff0c;和边权小于ttt的蓝边&#xff0c;问对于所有的点yyy&#xff0c;既能通过红边走向xxx&#xff0c;又能通过蓝边走向xxx&#xff0c;的点有多少个。 考虑对红边按照边权升序建立一颗kruskalkruska…

.NET 分布式自增Id组件(解决自动分配机器Id、时间回拨问题)

IdHelper是一个.NET&#xff08;支持.NET45或.NET Standard2&#xff09;生成分布式趋势自增Id组件&#xff0c;有两个版本&#xff1a;原始版为基于雪花Id&#xff08;不了解请自行百度&#xff09;方案&#xff0c;需要手动管理设置WorkerId&#xff1b;完美版在原始版的基础…

悲观锁与乐观锁

悲观锁 总是假设最坏的情况&#xff0c;每次取数据的时候都认为别人会来修改&#xff0c;所以每次取数据的时候都会上锁。其它线程想要取这份数据就必须拿到相应的锁&#xff08;共享资源每次只供一个线程使用&#xff0c;其它线程阻塞&#xff0c;用完之后转让给其他线程&…

1190 最小公倍数之和 V2

1190 最小公倍数之和 V2 ∑iablcm(i,b)∑iabibgcd⁡(i,b)b∑d∣b∑i⌈ad⌉bdi[gcd(i,bd)1]b∑d∣b∑k∣bdμ(k)k∑i⌈⌈ad⌉k⌉abkib∑T∣n∑i⌈aT⌉bTi∑k∣Tμ(k)kb∑T∣n(bT⌈aT⌉)(bT−⌈aT⌉1)2∑k∣Tμ(k)k设f(n)∑d∣nμ(d)d,f(1)1,f(p)1−p,f(pk)1−p,且为积性函数\sum…

有关 VS Code 的五大谣言,背后的真相到底是如何的?

2015 年 4 月 29 日&#xff0c;在微软 Build 2015 大会上&#xff0c;微软发布了 Visual Studio Code 第一个预览版本。随着 Visual Studio Code 有了越来越多的使用者&#xff0c;随之而来的各类谣言也层出不穷。让我们就来看看有哪些与 VS Code 相关的谣言&#xff0c;背后的…

Mysql数据库锁机制

一&#xff1a;概念介绍 MySQL数据库锁管理机制&#xff1a; SQL层实现的锁机制    Meta-data元数据锁&#xff1a;在table cache缓存里实现的&#xff0c;为DDL&#xff08;Data Definition Language&#xff09;提供隔离操作。一种特别的meta-data元数据类型&#xff0c;…

HDU 6340 Problem I. Delightful Formulas(伯努利数 + 积性函数反演)

Problem I. Delightful Formulas 大概就是照着题解抄了一遍吧&#xff0c;这道题太神仙了…… aiik,si∑j1iajcalc∑i1nsi[gcd⁡(i,n)1]∑d∣nμ(d)∑i1ndsida_i i ^ k, s_i \sum_{j 1} ^{i} a_j\\ calc\ \sum_{i 1} ^{n} s_i[\gcd(i, n) 1]\\ \sum_{d \mid n} \mu(d) \s…

干货|亲测有效的N倍学习效果笔记法

这里是Z哥的个人公众号每周五11&#xff1a;45 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「108」篇原创敬上大家好&#xff0c;我是Z哥。先祝大家中秋快乐。我猜你现在心情不错&#xff0c;毕竟小长假的第一天才开始&#xff0c;后面还有60个小时的假期&a…

Java偏向锁、轻量级锁、重量级锁

先Mark&#xff0c;后补充 参照&#xff1a; https://www.infoq.cn/article/java-se-16-synchronized http://www.cnblogs.com/paddix/p/5405678.html http://www.cnblogs.com/lzh-blogs/p/7477157.html

.NET Core 3.0 可卸载程序集原理简析

文章转载授权级别&#xff1a;A 预计阅读时间&#xff1a;8分钟 损失发量&#xff1a;不好统计因为最近在群里被问到如何理解 .NET Core 3.0 可卸载程序集&#xff0c;所以就写了这篇简单的分析。因为时间实在很少&#xff0c;这篇文章只简单的罗列了相关的代码&…

P4331 [BalticOI 2004]Sequence 数字序列(左偏树)

P4331 [BalticOI 2004]Sequence 数字序列 给定一个序列整数a1,a2,a3,…,an−1,ana_1, a_2, a_3, \dots, a_{n - 1}, a_na1​,a2​,a3​,…,an−1​,an​&#xff0c;要找一个整数序列bbb&#xff0c;满足b1<b2<b3<⋯<bn−1<bnb_1 < b_2 < b_3< \dots&…

.NetCore技术研究-ConfigurationManager在单元测试下的坑

最近在将原有代码迁移.NET Core, 代码的迁移基本很快&#xff0c;当然也遇到了不少坑&#xff0c;重构了不少&#xff0c;后续逐步总结分享给大家。今天总结分享一下ConfigurationManager遇到的一个问题。先说一下场景&#xff1a;迁移.NET Core后&#xff0c;已有的配置文件&a…

E 速度即转发(牛客挑战赛48)(树套树)

速度即转发 给定一个长度为nnn的数组aaa&#xff0c;里面元素为a1,a2,a3,…,an−1,ana_1, a_2, a_3, \dots, a_{n - 1}, a_na1​,a2​,a3​,…,an−1​,an​。 有两种操作&#xff1a; 修改apka_p kap​k。给定l,r,kl, r, kl,r,k&#xff0c;设S(x)∑ilrmax(ai−x,0)S(x) …