洛谷传送门
BZOJ传送门
题目描述
本题的故事发生在魔力之都,在这里我们将为你介绍一些必要的设定。 魔力之都可以抽象成一个 n n 个节点、 条边的无向连通图(节点的编号从 1 1 至 )。我们依次用 l,a l , a 描述一条边的长度、海拔。 作为季风气候的代表城市,魔力之都时常有雨水相伴,因此道路积水总是不可避免 的。由于整个城市的排水系统连通,因此有积水的边一定是海拔相对最低的一些边。我们用水位线来描述降雨的程度,它的意义是:所有海拔不超过水位线的边都是有积水的。
Yazid 是一名来自魔力之都的OIer,刚参加完ION2018 的他将踏上归程,回到他 温暖的家。 Yazid 的家恰好在魔力之都的 11 号节点。对于接下来 QQ 天,每一天Yazid 都会告诉你他的出发点 vv ,以及当天的水位线pp。 每一天,Yazid 在出发点都拥有一辆车。这辆车由于一些故障不能经过有积水的边。 Yazid 可以在任意节点下车,这样接下来他就可以步行经过有积水的边。但车会被留在他下车的节点并不会再被使用。 需要特殊说明的是,第二天车会被重置,这意味着:
- 车会在新的出发点被准备好。
- Yazid 不能利用之前在某处停放的车。 Yazid 非常讨厌在雨天步行,因此他希望在完成回家这一目标的同时,最小化他步行经过的边的总长度。请你帮助 Yazid 进行计算。 本题的部分测试点将强制在线,具体细节请见【输入格式】和【子任务】。
输入输出格式
输入格式:
单个测试点中包含多组数据。输入的第一行为一个非负整数 T T ,表示数据的组数。
接下来依次描述每组数据,对于每组数据:
第一行 个非负整数 n,m n , m ,分别表示节点数、边数。
接下来 m m 行,每行 个正整数 u,v,l,a u , v , l , a ,描述一条连接节点 u,v u , v 的、长度为 l l 、海拔为 的边。 在这里,我们保证 1≤u,v≤n 1 ≤ u , v ≤ n 。
接下来一行 3 3 个非负数 ,其中 Q Q 表示总天数, 是一个会在下面被用到的系数, S S 表示的是可能的最高水位线。
接下来 行依次描述每天的状况。每行
2
2
个整数 描述一天:
这一天的出发节点为
v=(v0+K×lastans−1)modn+1
v
=
(
v
0
+
K
×
l
a
s
t
a
n
s
−
1
)
mod
n
+
1
。
这一天的水位线为
p=(p0+K×lastans)mod(S+1)
p
=
(
p
0
+
K
×
l
a
s
t
a
n
s
)
mod
(
S
+
1
)
。
其中 lastans
表示上一天的答案(最小步行总路程)。特别地,我们规定第
1
1
天时 lastans = 0
。 在这里,我们保证 。
对于输入中的每一行,如果该行包含多个数,则用单个空格将它们隔开。
输出格式:
依次输出各组数据的答案。对于每组数据:
- 输出 Q Q 行每行一个整数,依次表示每天的最小步行总路程。
输入输出样例
输入样例#1:
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
输出样例#1:
0
50
200
50
150
输入样例#2:
1
5 5
1 2 1 2
2 3 1 2
4 3 1 2
5 3 1 2
1 5 2 1
4 1 3
5 1
5 2
2 0
4 0
输出样例#2:
0
2
3
1
说明
【样例1 解释】 第一天没有降水,Yazid 可以坐车直接回到家中。
第二天、第三天、第四天的积水情况相同,均为连接 号节点的边、连接 3,4 3 , 4 号点的边有积水。
对于第二天,Yazid 从 2 2 号点出发坐车只能去往 号节点,对回家没有帮助。因此 Yazid 只能纯靠徒步回家。
对于第三天,从 4 4 号节点出发的唯一一条边是有积水的,车也就变得无用了。Yazid只能纯靠徒步回家。
对于第四天,Yazid 可以坐车先到达 号节点,再步行回家。
第五天所有的边都积水了,因此Yazid 只能纯靠徒步回家。
所有测试点均保证 T≤3 T ≤ 3 ,所有测试点中的所有数据均满足如下限制:
- n≤2×105 n ≤ 2 × 10 5 , m≤4×105 m ≤ 4 × 10 5 , Q≤4×105 Q ≤ 4 × 10 5 , K∈{0,1} K ∈ { 0 , 1 } , 1≤S≤109 1 ≤ S ≤ 10 9 。
- 对于所有边: l≤104 l ≤ 10 4 , a≤109 a ≤ 10 9 。
- 任意两点之间都直接或间接通过边相连。
为了方便你快速理解,我们在表格中使用了一些简单易懂的表述。在此,我们对这些内容作形式化的说明:
- 图形态:对于表格中该项为“一棵树”或“一条链”的测试点,保证 m=n−1 m = n − 1 。 除此之外,这两类测试点分别满足如下限制:
- 一棵树:保证输入的图是一棵树,即保证边不会构成回路。
- 一条链:保证所有边满足 u+1=v u + 1 = v 。
- 海拔:对于表格中该项为“一种”的测试点,保证对于所有边有 a=1 a = 1 。
- 强制在线:对于表格中该项为“是”的测试点,保证 K=1 K = 1 ;如果该项为“否”, 则有 K=0 K = 0 。
- 对于所有测试点,如果上述对应项为“不保证”,则对该项内容不作任何保证。
n n | Q= Q = | 测试点 | 形态 | 海拔 | 强制在线 | |
---|---|---|---|---|---|---|
≤1 ≤ 1 | ≤0 ≤ 0 | 0 0 | 1 | 不保证 | 一种 | 否 |
≤10 ≤ 10 | 10 10 | 2 | 不保证 | 一种 | 否 | |
≤50 ≤ 50 | ≤150 ≤ 150 | 100 100 | 3 | 不保证 | 一种 | 否 |
≤100 ≤ 100 | ≤300 ≤ 300 | 200 200 | 4 | 不保证 | 一种 | 否 |
≤1500 ≤ 1500 | ≤4000 ≤ 4000 | 2000 2000 | 5 | 不保证 | 一种 | 否 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 6 | 不保证 | 一种 | 否 |
≤1500 ≤ 1500 | =n−1 = n − 1 | 2000 2000 | 7 | 一条链 | 不保证 | 否 |
≤1500 ≤ 1500 | =n−1 = n − 1 | 2000 2000 | 8 | 一条链 | 不保证 | 否 |
≤1500 ≤ 1500 | =n−1 = n − 1 | 2000 2000 | 9 | 一条链 | 不保证 | 否 |
≤200000 ≤ 200000 | =n−1 = n − 1 | 100000 100000 | 10 | 一棵树 | 不保证 | 否 |
≤200000 ≤ 200000 | =n−1 = n − 1 | 100000 100000 | 11 | 一棵树 | 不保证 | 是 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 12 | 不保证 | 不保证 | 否 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 13 | 不保证 | 不保证 | 否 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 14 | 不保证 | 不保证 | 否 |
≤1500 ≤ 1500 | ≤4000 ≤ 4000 | 2000 2000 | 15 | 不保证 | 不保证 | 是 |
≤1500 ≤ 1500 | ≤4000 ≤ 4000 | 2000 2000 | 16 | 不保证 | 不保证 | 是 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 17 | 不保证 | 不保证 | 是 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 100000 100000 | 18 | 不保证 | 不保证 | 是 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 400000 400000 | 19 | 不保证 | 不保证 | 是 |
≤200000 ≤ 200000 | ≤400000 ≤ 400000 | 400000 400000 | 20 | 不保证 | 不保证 | 是 |
解题分析
先 dj d j 跑出每个点到 1 1 号点的最短路, 再按边的高度建立最大生成树。 生成树的一棵子树可以通过不低于子树根节点的路径互相到达, 也就是坐车可以到达的区域。 我们可以维护子树内点到号点的最小值, 并通过倍增找到这样的一个最高点。
总复杂度 O(TN(log(N)) O ( T N ( l o g ( N ) ) 。
代码如下:
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdlib>
#include <ext/pb_ds/priority_queue.hpp>
#include <limits.h>
#define R register
#define IN inline
#define W while
#define MX 800050
#define gc getchar()
#define ll long long
template <class T>
IN void in(T &x)
{
x = 0; R char c = gc;
for (; !isdigit(c); c = gc);
for (; isdigit(c); c = gc)
x = (x << 1) + (x << 3) + c - 48;
}
int dot, line, Q, cnt1, cnt2, bd, arr, K, S;
int head1[MX], head2[MX], h[MX], to[MX][22], bel[MX];
ll dis[MX], dmin[MX];
bool vis[MX];
struct EDGE {int to, nex, len, height;} ed[MX << 1], edge[MX << 1];
struct Edge {int to, nex;} e[MX << 1];
struct INFO {int now, dis;};
IN bool operator < (const EDGE &x, const EDGE &y)
{return x.height == y.height ? x.len < y.len : x.height > y.height;}
IN bool operator < (const INFO &x, const INFO &y)
{return x.dis > y.dis;}
IN void add1(R int from, R int to, R int len, R int height)
{edge[++cnt1] = {to, head1[from], len, height}; head1[from] = cnt1;}
IN void add2(R int from, R int to)
{e[++cnt2] = {to, head2[from]}, head2[from] = cnt2;}
__gnu_pbds :: priority_queue <INFO> :: point_iterator hdle[MX];
__gnu_pbds :: priority_queue <INFO> que;
int find (R int now) {return now == bel[now] ? now : bel[now] = find(bel[now]);}
void DJ()
{
std::memset(dis, 63, sizeof(dis));
dis[1] = 0; que.push({1, 0});
INFO cur; R int now, i;
W (!que.empty())
{
cur = que.top(); que.pop();
now = cur.now;
for (i = head1[now]; i; i = edge[i].nex)
{
if(dis[edge[i].to] > dis[now] + edge[i].len)
{
dis[edge[i].to] = dis[now] + edge[i].len;
if(!vis[edge[i].to]) hdle[edge[i].to] = que.push({edge[i].to, dis[edge[i].to]}), vis[edge[i].to] = true;
else que.modify(hdle[edge[i].to], {edge[i].to, dis[edge[i].to]});
}
}
}
}
void build()
{
int b1, b2;
std::sort(ed + 1, ed + 1 + line);
for (R int i = 1; i <= bd; ++i) bel[i] = i;
for (R int i = 1; i <= line; ++i)
{
b1 = find(ed[i].to), b2 = find(ed[i].nex);
if(b1 ^ b2)
{
++arr; bel[b1] = bel[b2] = arr + dot;
h[arr + dot] = ed[i].height;
add2(arr + dot, b1), add2(arr + dot, b2);
if(arr == dot - 1) break;
}
}
}
void DFS(R int now, R int fa)
{
if(now <= dot) dmin[now] = dis[now];
else dmin[now] = LONG_LONG_MAX;
to[now][0] = fa;
for (R int i = 1; i <= 21; ++i)
{
if(to[to[now][i - 1]][i - 1])
to[now][i] = to[to[now][i - 1]][i - 1];
else break;
}
for (R int i = head2[now]; i; i = e[i].nex)
DFS(e[i].to, now), dmin[now] = std::min(dmin[e[i].to], dmin[now]);
}
IN ll query(R int now, R int lim)
{
for (R int i = 21; ~i; --i)
if(h[to[now][i]] > lim) now = to[now][i];
return dmin[now];
}
int main(void)
{
int a, b, c, d;
ll lastans;
int T;
in(T);
W (T--)
{
std::memset(head1, arr = cnt1 = 0, sizeof(head1));
std::memset(head2, cnt2 = 0, sizeof(head2));
std::memset(vis, lastans = 0, sizeof(vis));
std::memset(to, 0, sizeof(to));
std::memset(h, 0, sizeof(h));
in(dot), in(line); bd = dot * 2 - 1;
for (R int i = 1; i <= line; ++i)
{
in(a), in(b), in(c), in(d);
add1(a, b, c, d), add1(b, a, c, d);
ed[i] = {a, b, c, d};
}
DJ();
build(); DFS(bd, 0);
in(Q), in(K), in(S);
W (Q--)
{
in(a), in(b);
a = (a + lastans * K - 1) % dot + 1;
b = (b + lastans * K) % (S + 1);
printf("%lld\n", lastans = query(a, b));
}
}
}