\(\color{#0066ff}{题目描述}\)

本题的故事发生在魔力之都,在这里我们将为你介绍一些必要的设定。 魔力之都可以抽象成一个 n 个节点、m 条边的无向连通图(节点的编号从 1 至 n)。我们依次用 l,a 描述一条边的长度、海拔。 作为季风气候的代表城市,魔力之都时常有雨水相伴,因此道路积水总是不可避免 的。由于整个城市的排水系统连通,因此有积水的边一定是海拔相对最低的一些边。我们用水位线来描述降雨的程度,它的意义是:所有海拔不超过水位线的边都是有积水的。

Yazid 是一名来自魔力之都的OIer,刚参加完ION2018 的他将踏上归程,回到他 温暖的家。 Yazid 的家恰好在魔力之都的 1 号节点。对于接下来 Q 天,每一天Yazid 都会告诉你他的出发点 v ,以及当天的水位线p。 每一天,Yazid 在出发点都拥有一辆车。这辆车由于一些故障不能经过有积水的边。 Yazid 可以在任意节点下车,这样接下来他就可以步行经过有积水的边。但车会被留在他下车的节点并不会再被使用。 需要特殊说明的是,第二天车会被重置,这意味着:

车会在新的出发点被准备好。

Yazid 不能利用之前在某处停放的车。

Yazid 非常讨厌在雨天步行,因此他希望在完成回家这一目标的同时,最小化他步行经过的边的总长度。请你帮助 Yazid 进行计算。 本题的部分测试点将强制在线,具体细节请见【输入格式】和【子任务】。

\(\color{#0066ff}{输入格式}\)

单个测试点中包含多组数据。输入的第一行为一个非负整数T,表示数据的组数。

接下来依次描述每组数据,对于每组数据:

第一行 2 个非负整数 n,m,分别表示节点数、边数。

接下来 m 行,每行 4 个正整数 u,v,l,a,描述一条连接节点 u,v 的、长度为 l、海拔为 a 的边。 在这里,我们保证\(1 \leq u,v \leq n\)。

接下来一行 3 个非负数 Q,K,S ,其中 Q 表示总天数,\(K \in {0,1}\) 是一个会在下面被用到的系数,S 表示的是可能的最高水位线。

接下来 Q 行依次描述每天的状况。每行 2 个整数 \(v_0; p_0\) 描述一天:

这一天的出发节点为\(v = (v_0 + K \times \mathrm{lastans} - 1) \bmod n + 1\)

这一天的水位线为\(p = (p_0 + K \times \mathrm{lastans}) \bmod (S + 1)\)

其中 lastans 表示上一天的答案(最小步行总路程)。特别地,我们规定第 1 天时 lastans = 0。 在这里,我们保证\(1 \leq v_0 \leq n,0 \leq p_0 \leq S\)

对于输入中的每一行,如果该行包含多个数,则用单个空格将它们隔开。

\(\color{#0066ff}{输出格式}\)

依次输出各组数据的答案。对于每组数据:

输出 Q 行每行一个整数,依次表示每天的最小步行总路程。

\(\color{#0066ff}{输入样例}\)

1
4 3
1 2 50 1
2 3 100 2
3 4 50 1
5 0 2
3 0
2 1
4 1
3 1
3 2

\(\color{#0066ff}{输出样例}\)

0
50
200
50
150

\(\color{#0066ff}{数据范围与提示}\)

\(\color{#0066ff}{题解}\)

Kruskal重构树

在跑Kruskal的时候,对于加进去的边

建立一个虚点

虚点的点权等于当前边的边权

将虚点与两个并查集的根相连

这样最后会构成一棵树

叶子节点为原来的节点

其余节点为虚点

可以发现这是一个堆

对于本题来说

我们按海拔排序,弄成一个小根堆

对于询问的一个点x

倍增在树上跳到一个最浅的点,使得海拔恰好高于水平面

那么这个点的子树内的所有叶子节点都可以互相到达

我们输出这些叶子节点到1的最短距离即可

这个可以在开始时跑dij预处理

#include<bits/stdc++.h>
#define LL long long
LL in() {
char ch; int x = 0, f = 1;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
return x * f;
}
const int maxn = 4e5 + 100;
const int inf = 0x7fffffff;
struct E {
int x, y, z, h;
const bool operator < (const E &b) const {
return h > b.h;
}
}e[maxn];
int r[maxn];
int dep[maxn], f[maxn][26], n, m, min[maxn], val[maxn], ans;
int cnt, dis[maxn];
using std::vector;
using std::pair;
using std::make_pair;
std::priority_queue<pair<int, int>, vector<pair<int, int> >, std::greater<pair<int, int> > > q;
struct node {
int to, dis;
node *nxt;
node(int to = 0, int dis = 0, node *nxt = NULL): to(to), dis(dis), nxt(nxt) {}
void *operator new (size_t) {
static node *S = NULL, *T = NULL;
return (S == T) && (T = (S = new node[1024]) + 1024), S++;
}
};
node *head[maxn], *h[maxn];
bool vis[maxn];
void add(int from, int to, int dis) {
h[from] = new node(to, dis, h[from]);
}
void add(int from, int to) {
head[from] = new node(to, 0, head[from]);
}
void dfs(int x, int fa) {
f[x][0] = fa;
dep[x] = dep[fa] + 1;
for(node *i = head[x]; i; i = i->nxt)
if(i->to != fa) dfs(i->to, x), min[x] = std::min(min[x], min[i->to]);
}
void dij() {
for(int i = 1; i <= n; i++) dis[i] = inf, vis[i] = 0;
q.push(make_pair(dis[1] = 0, 1));
while(!q.empty()) {
int tp = q.top().second;
q.pop();
if(vis[tp]) continue;
vis[tp] = true;
for(node *i = h[tp]; i; i = i->nxt)
if(dis[i->to] > dis[tp] + i->dis)
q.push(make_pair(dis[i->to] = dis[tp] + i->dis, i->to));
}
} void beizeng() {
dfs(cnt, 0);
for(int j = 1; j <= 25; j++)
for(int i = 1; i <= cnt; i++)
f[i][j] = f[f[i][j - 1]][j - 1];
}
int findset(int x) { return x == r[x]? r[x] : r[x] = findset(r[x]); }
int query(int x, int y) {
for(int i = 25; i >= 0; i--) if(val[f[x][i]] > y) x = f[x][i];
return min[x];
}
void Kruskal() {
for(int i = 1; i <= n; i++) r[i] = i;
std::sort(e + 1, e + m + 1);
int tot = 0;
cnt = n;
for(int i = 1; i <= m; i++) {
int x = findset(e[i].x);
int y = findset(e[i].y);
if(x != y) {
r[x] = r[y] = ++cnt;
r[cnt] = cnt;
add(cnt, x);
add(cnt, y);
val[cnt] = e[i].h;
tot++;
}
if(tot == n - 1) break;
}
for(int i = n + 1; i <= cnt; i++) min[i] = inf;
beizeng();
int T = in(), k = in(), s = in();
int x, y;
ans = 0;
while(T --> 0) {
x = (in() + k * ans - 1) % n + 1;
y = (in() + k * ans) % (s + 1);
printf("%d\n", ans = query(x, y));
}
} int main() {
for(int T = in(); T --> 0;) {
n = in(), m = in();
memset(head, 0, sizeof head);
memset(h, 0, sizeof h);
for(int i = 1; i <= m; i++) {
E &o = e[i];
o.x = in(), o.y = in(), o.z = in(), o.h = in();
add(o.x, o.y, o.z);
add(o.y, o.x, o.z);
}
dij();
for(int i = 1; i <= n; i++) min[i] = dis[i];
Kruskal();
}
return 0;
}

P4768 [NOI2018]归程的更多相关文章

  1. Luogu P4768 [NOI2018]归程(Dijkstra+Kruskal重构树)

    P4768 [NOI2018]归程 题面 题目描述 本题的故事发生在魔力之都,在这里我们将为你介绍一些必要的设定. 魔力之都可以抽象成一个 \(n\) 个节点. \(m\) 条边的无向连通图(节点的编 ...

  2. P4768 [NOI2018]归程(kruskal 重构树)

    洛谷P4768 [NOI2018]归程 LOJ#2718.「NOI2018」归程 用到 kruskal 重构树,所以先说这是个啥 显然,这和 kruskal 算法有关系 (废话 这个重构树是一个有点权 ...

  3. [洛谷P4768] [NOI2018]归程 (kruskal重构树模板讲解)

    洛谷题目链接:[NOI2018]归程 因为题面复制过来有点炸格式,所以要看题目就点一下链接吧\(qwq\) 题意: 在一张无向图上,每一条边都有一个长度和海拔高度,小\(Y\)的家在\(1\)节点,并 ...

  4. 洛谷P4768 [NOI2018]归程 [可持久化并查集,Dijkstra]

    题目传送门 归程 格式难调,题面就不放了. 分析: 之前同步赛的时候反正就一脸懵逼,然后场场暴力大战,现在呢,还是不会$Kruskal$重构树,于是就拿可持久化并查集做. 但是之前做可持久化并查集的时 ...

  5. 洛谷P4768 [NOI2018]归程(Kruskal重构树)

    题意 直接看题目吧,不好描述 Sol 考虑暴力做法 首先预处理出从$1$到每个节点的最短路, 对于每次询问,暴力的从这个点BFS,从能走到的点里面取$min$ 考虑如何优化,这里要用到Kruskal重 ...

  6. Luogu P4768 [NOI2018]归程

    题目链接 \(Click\) \(Here\) \(Kruskal\)重构树的好题.想到的话就很好写,想不到乱搞的难度反而相当高. 按照点的水位,建出来满足小根队性质的\(Kruskal\)重构树,这 ...

  7. 洛谷P4768 [NOI2018]归程(可持久化并查集,最短路)

    闲话 一个蒟蒻,在网络同步赛上进行了这样的表演-- T2组合计数不会,T3字符串数据结构不会,于是爆肝T1 一开始以为整个地图都有车,然后写了2h+的树套树,终于发现样例过不去 然后写可持久化并查集D ...

  8. 洛谷 P4768 [NOI2018]归程

    洛谷 361行代码的由来 数据分治大发好啊- NOI的签到题,可怜我在家打了一下午才搞了80分. 正解应该是kruskal重构树或排序+可持久化并查集. 我就分点来讲暴力80分做法吧(毕竟正解我也没太 ...

  9. 洛谷P4768 [NOI2018]归程(克鲁斯卡尔重构树+最短路)

    传送门 前置技能,克鲁斯卡尔重构树 我们按道路的高度建一个最大生成树,然后建好克鲁斯卡尔重构树 那么我们需要知道一颗子树内到1点距离最近是多少(除此之外到子树内任何一个点都不需要代价) 可以一开始直接 ...

随机推荐

  1. nginx upstream的几种配置方式

    nginx 的upstream目前支持4种方式的分配 1.轮询(默认) 每个请求按时间顺序逐一分配到不同的后端服务器 ,如果后端服务器down掉,能自动剔除. 2.weight指定轮询几率,weigh ...

  2. 问题:oracle decode;结果:oracle中的decode的使用

    oracle中的decode的使用 Oracle 中 decode 函数用法 含义解释:decode(条件,值1,返回值1,值2,返回值2,...值n,返回值n,缺省值) 该函数的含义如下:IF 条件 ...

  3. sort排序应用

      private static int SortTestObj2Compare(TestSortClass obj1, TestSortClass obj2)        {            ...

  4. day70-oracle PLSQL_02光标

    涨工资之前员工的工资. 如果PLSQL程序没有commit的话,命令行这边的客户端是无法读到的.这是oracle数据库的隔离级别. 为什么在PLSQL程序中commit之后还是不行呢? PLSQL程序 ...

  5. js面试题知识点全解(一原型和原型链)

    1.如何准确判断一个变量是数组类型2.写一个原型链继承的例子3.描述new一个对象的过程4.zepto(或其他框架)源码中如何使用原型链知识点:1.构造函数2.构造函数-扩展3.原型规则和示例4.原型 ...

  6. HDU 6395(2018多校第7场1010)Sequence

    不久前做过POJ3070,所以知道这题要用矩阵快速幂优化,但是这个题的递推公式中有一项⌊p/n⌋,场上就不会了... 下来才知道要用分块矩阵快速幂,因为⌊p/n⌋最多有2√p块,可以对每一块使用快速幂 ...

  7. opencv3.2加载图片

    #include <opencv2/opencv.hpp>#include <iostream> using namespace cv;using namespace std; ...

  8. Python Matplotlib.plot Update image Questions

    1. 最近在测试一款设备,采集了一些设备后需要一帧一帧显示图像,经常使用Python,所以选用了Matplotlib进行图像操作 数据结构: timesatamp polar_distance hor ...

  9. C#调用C++类库的几种方式

    1.  直接调用C++类库中的公共方法 使用DllImport特性对方法进行调用,比如一个C++类库SampleCppWrapper.dll中的公共方法: extern "C" _ ...

  10. OpenGL编程

    一.简介 OpenGL™ 是行业领域中最为广泛接纳的 2D/3D 图形 API, 其自诞生至今已催生了各种计算机平台及设备上的数千优秀应用程序.OpenGL™ 是独立于视窗操作系统或其它操作系统的,亦 ...