Submit: 1574 Solved: 651
[Submit][Status][Discuss]
Description
给一个N个点M条边的连通无向图，满足每条边最多属于一个环，有Q组询问，每次询问两点之间的最短路径。
Input
输入的第一行包含三个整数，分别表示N和M和Q 下接M行，每行三个整数v，u，w表示一条无向边v-u，长度为w 最后Q行，每行两个整数v，u表示一组询问
Output
输出Q行，每行一个整数表示询问的答案
Sample Input
9 10 2
1 2 1
1 4 1
3 4 1
2 3 1
3 7 1
7 8 2
7 9 2
1 5 3
1 6 4
5 6 1
1 9
5 7
1 2 1
1 4 1
3 4 1
2 3 1
3 7 1
7 8 2
7 9 2
1 5 3
1 6 4
5 6 1
1 9
5 7
Sample Output
5
6
6
HINT
对于100%的数据，N<=10000，Q<=10000
Source
圆方树好毒瘤神奇啊qwq。
对于这题来说的话，首先把圆方树弄出来
考虑询问的两个点，如果他们的lca是圆点，那么直接树上前缀和相加
如果不是圆点，那么他们的lca一定是在一个环内，
我们可以维护出环内任意一点到环的根节点（也就是$dfn$最小的节点）的最短路径
然后分类讨论一下即可
对于一个点如何跳的它lca所在的环上，这里有个骚操作。
我们考虑树链剖分的过程
如果这个节点是重儿子，因为重儿子的dfs序是与父节点连续的，所以要求的点就是lca的dfs中对应的下一个节点
如果不是重儿子，那么我们一直跳，直到跳到lca处位置，那么要求的点就是来时离lca最近的那个节点
#include<cstdio> #include<vector> #include<cstring> using namespace std; const int MAXN = 1e6 + 10; inline int read() { char c = getchar(); int x = 0, f = 1; while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();} while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar(); return x * f; } int N, M, Q; struct Edge { int u, v, w, nxt; }E[MAXN]; int head[MAXN], num = 1; inline void AddEdge(int x, int y, int z) { E[num] = (Edge){x, y, z, head[x]}; head[x] = num++; } vector<Edge> vec[MAXN]; inline void add_edge(int x, int y, int z) { vec[x].push_back((Edge){x, y, z}); vec[y].push_back((Edge){y, x, z}); } int low[MAXN], dfn[MAXN], times, fa[MAXN], cdis[MAXN], dis[MAXN], tot; inline void Build(int x, int y, int w) { for(int i = y; i != x; i = fa[i]) cdis[i] = w, w += dis[i]; cdis[++tot] = w; add_edge(tot, x, 0);//建方点 for(int i = y; i != x; i = fa[i]) add_edge(tot, i, min(cdis[i], w - cdis[i])); } void Tarjan(int x, int _fa) { dfn[x] = low[x] = ++times; fa[x] = _fa; for(int v, i = head[x]; i != -1; i = E[i].nxt) { if((v = E[i].v) == _fa) continue; if(!dfn[v]) dis[v] = E[i].w, Tarjan(v, x), low[x] = min(low[x], low[v]); else low[x] = min(low[x], dfn[v]); if(low[v] > dfn[x]) add_edge(x, v, E[i].w); } // can I take this into for ? for(int v, i = head[x]; i != -1; i = E[i].nxt) if(fa[(v = E[i].v)] != x && dfn[v] > dfn[x])//如果不是负边，那么一定出现了环 Build(x, v, E[i].w); } int siz[MAXN], son[MAXN], top[MAXN], deep[MAXN], Index = 0, point[MAXN]; void dfs1(int x, int _fa) { fa[x] = _fa; siz[x] = 1; for(int i = 0, v; i < vec[x].size(); i++) { if((v = vec[x][i].v) == _fa) continue; dis[v] = dis[x] + vec[x][i].w; deep[v] = deep[x] + 1; dfs1(v, x); siz[x] += siz[v]; if(siz[v] > siz[son[x]]) son[x] = v; } } void dfs2(int x, int topf) { top[x] = topf; point[dfn[x] = ++Index] = x; if(!son[x]) return ; dfs2(son[x], topf); for(int i = 0, v; i < vec[x].size(); i++) if(!top[v = vec[x][i].v]) dfs2(v, v); } int LCA(int x, int y) { while(top[x] != top[y]) { if(deep[top[x]] < deep[top[y]]) swap(x, y); x = fa[top[x]]; } if(deep[x] < deep[y]) swap(x, y); return y; } int Jump(int x, int lca) { int las; while(top[x] != top[lca]) las = top[x], x = fa[top[x]];//las = top[x] not x return x == lca ? las : point[dfn[lca] + 1];//这里要写lca } int abs(int x) { return x < 0 ? -x : x; } int main() { #ifdef WIN32 freopen("a.in", "r", stdin); freopen("a.out", "w", stdout); #endif tot = N = read(); M = read(); Q = read(); memset(head, -1, sizeof(head)); for(int i = 1; i <= M; i++) { int x = read(), y = read(), z = read(); AddEdge(x, y, z); AddEdge(y, x, z); } Tarjan(1, 0); dis[1] = 0; deep[1] = 1; //for(int i = 1; i <= N; i++) printf("%d ", dfn[i]); puts(""); //for(int i = 1; i <= N; i++) printf("%d ", vec[i].size()); puts(""); //printf("%d ", tot); dfs1(1, 0); dfs2(1, 1); //for(int i = 1; i <= tot; i++) printf("%d ", top[i]); puts(""); while(Q--) { int x = read(), y = read(); int lca = LCA(x, y); if(lca <= N) {printf("%d ", dis[x] + dis[y] - (dis[lca] << 1)); continue;} int p1 = Jump(x, lca); int p2 = Jump(y, lca); int ans = dis[x] - dis[p1] + dis[y] - dis[p2] + min(abs(cdis[p2] - cdis[p1]), cdis[lca] - abs(cdis[p2] - cdis[p1])); printf("%d ", ans); } return 0; }