### 代码

```#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#include <queue>

using namespace std;

const int MaxN = 100000 + 5, MaxM = 200000 + 5, MaxLog = 20, INF = 999999999;

int n, m, Index;
int Father[MaxN], Depth[MaxN], Top[MaxN], Size[MaxN], Son[MaxN], Pos[MaxN];
int d[MaxN], D[MaxN * 4], Jump[MaxN][MaxLog + 3];

struct Edge
{
int u, v, w;
bool Mark;
Edge *Next;
} E[MaxM * 2], *P = E, *Pre[MaxN], *Point[MaxN];

inline void AddEdge(int x, int y, int z) {
++P; P -> u = x; P -> v = y; P -> w = z; P -> Mark = false;
P -> Next = Point[x]; Point[x] = P;
}

struct ES
{
int x, y;
ES() {}
ES(int a, int b) {
x = a; y = b;
}
};

struct Cmp
{
bool operator () (ES a, ES b) {
return a.y > b.y;
}
};

priority_queue<ES, vector<ES>, Cmp> Q;

bool Visit[MaxN];

void Dijkstra() {
while (!Q.empty()) Q.pop();
for (int i = 1; i <= n; ++i) {
d[i] = INF; Visit[i] = false;
}
d[1] = 0;
for (int i = 1; i <= n; ++i) Q.push(ES(i, d[i]));
ES Now;
int x;
while (!Q.empty()) {
Now = Q.top(); Q.pop();
x = Now.x;
if (Visit[x]) continue;
Visit[x] = true;
for (Edge *j = Point[x]; j; j = j -> Next) {
if (d[x] + (j -> w) < d[j -> v]) {
d[j -> v] = d[x] + j -> w;
if (Pre[j -> v] != NULL) Pre[j -> v] -> Mark = false;
Pre[j -> v] = j;
j -> Mark = true;
Q.push(ES(j -> v, d[j -> v]));
}
}
}
}

int DFS_1(int x, int Dep, int Fa) {
Depth[x] = Dep; Father[x] = Fa;
Size[x] = 1;
int SonSize, MaxSonSize;
SonSize = MaxSonSize = 0;
for (Edge *j = Point[x]; j; j = j -> Next) {
if (j -> v == Fa || j -> Mark == false) continue;
SonSize = DFS_1(j -> v, Dep + 1, x);
if (SonSize > MaxSonSize) {
MaxSonSize = SonSize;
Son[x] = j -> v;
}
Size[x] += SonSize;
}
return Size[x];
}

void DFS_2(int x) {
if (x == 0) return;
if (x == Son[Father[x]]) Top[x] = Top[Father[x]];
else Top[x] = x;
Pos[x] = ++Index;
DFS_2(Son[x]);
for (Edge *j = Point[x]; j; j = j -> Next) {
if (j -> v == Father[x] || j -> v == Son[x] || j -> Mark == false) continue;
DFS_2(j -> v);
}
}

void Build_Tree(int x, int s, int t) {
D[x] = INF;
if (s == t) return;
int m = (s + t) >> 1;
Build_Tree(x << 1, s, m);
Build_Tree(x << 1 | 1, m + 1, t);
}

void Init_LCA() {
for (int i = 1; i <= n; ++i) Jump[i][0] = Father[i];
for (int j = 1; j <= MaxLog; ++j) {
for (int i = 1; i <= n; ++i) {
if (Depth[i] < (1 << j)) continue;
Jump[i][j] = Jump[Jump[i][j - 1]][j- 1];
}
}
}

int LCA(int x, int y) {
int Dif;
if (Depth[x] < Depth[y]) swap(x, y);
Dif = Depth[x] - Depth[y];
if (Dif) {
for (int i = 0; i <= MaxLog; ++i) {
if (Dif & (1 << i)) x = Jump[x][i];
}
}
if (x == y) return x;
for (int i = MaxLog; i >= 0; --i) {
if (Jump[x][i] != Jump[y][i]) {
x = Jump[x][i];
y = Jump[y][i];
}
}
return Father[x];
}

inline int gmin(int a, int b) {return a < b ? a : b;}

void Paint(int x, int Num) {
if (Num >= D[x]) return;
D[x] = Num;
}

void PushDown(int x) {
if (D[x] == INF) return;
Paint(x << 1, D[x]);
Paint(x << 1 | 1, D[x]);
D[x] = INF;
}

void Change(int x, int s, int t, int l, int r, int Num) {
if (l <= s && r >= t) {
Paint(x, Num);
return;
}
PushDown(x);
int m = (s + t) >> 1;
if (l <= m) Change(x << 1, s, m, l, r, Num);
if (r >= m + 1) Change(x << 1 | 1, m + 1, t, l, r, Num);
}

void EChange(int x, int y, int z) {
int fx, fy;
fx = Top[x]; fy = Top[y];
while (fx != fy) {
Change(1, 1, n, Pos[fx], Pos[x], z);
x = Father[fx];
fx = Top[x];
}
if (x != y) Change(1, 1, n, Pos[y] + 1, Pos[x], z);
}

int Get(int x, int s, int t, int p) {
if (s == t) return D[x];
PushDown(x);
int m = (s + t) >> 1;
int ret;
if (p <= m) ret = Get(x << 1, s, m, p);
else ret = Get(x << 1 | 1, m + 1, t, p);
return ret;
}

int main()
{
scanf("%d%d", &n, &m);
int a, b, c;
for (int i = 1; i <= m; ++i) {
scanf("%d%d%d", &a, &b, &c);
}
Dijkstra();
DFS_1(1, 0, 0);
Index = 0;
DFS_2(1);
Build_Tree(1, 1, n);
Init_LCA();
int t;
for (Edge *j = E + 1; ; ++j) {
if (j -> Mark) continue;
t = LCA(j -> u, j -> v);
EChange(j -> v, t, d[j -> u] + j -> w + d[j -> v]);
if (j == P) break;
}
int Temp;
for (int i = 2; i <= n; ++i) {
Temp = Get(1, 1, n, Pos[i]);
if (Temp < INF) printf("%d\n", Temp - d[i]);
else printf("-1\n");
}
return 0;
}
```

## [BZOJ 1576] [Usaco2009 Jan] 安全路经Travel 【树链剖分】的更多相关文章

1. bzoj 1576: [Usaco2009 Jan]安全路经Travel 树链剖分

1576: [Usaco2009 Jan]安全路经Travel Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 665  Solved: 227[Sub ...

2. BZOJ1576: [Usaco2009 Jan]安全路经Travel(树链剖分)

Description Input * 第一行: 两个空格分开的数, N和M * 第2..M+1行: 三个空格分开的数a_i, b_i,和t_i Output * 第1..N-1行: 第i行包含一个数 ...

3. bzoj 1576 [Usaco2009 Jan]安全路经Travel（树链剖分，线段树）

[题意] 给定一个无向图,找到1-i所有的次短路经,要求与最短路径的最后一条边不重叠. [思路] 首先用dijkstra算法构造以1为根的最短路树. 将一条无向边看作两条有向边,考察一条不在最短路树上 ...

4. bzoj 1576: [Usaco2009 Jan]安全路经Travel【spfa+树链剖分+线段树】

这几天写USACO水题脑子锈住了--上来就贪心,一交就WA 事实上这个是一个叫最短路树的东西,因为能保证只有一条最短路,所以所有最短路合起来是一棵以1为根的树,并且在这棵树上,每个点被精灵占据的路是它 ...

5. BZOJ.1576.[Usaco2009 Jan]安全路经Travel(树形DP 并查集)

题目链接 BZOJ 洛谷 先求最短路树.考虑每一条非树边(u,v,len),设w=LCA(u,v),这条边会对w->v上的点x(x!=w)有dis[u]+dis[v]-dis[x]+len的距离 ...

6. bzoj 1576: [Usaco2009 Jan]安全路经Travel——并查集+dijkstra

Description Input * 第一行: 两个空格分开的数, N和M * 第2..M+1行: 三个空格分开的数a_i, b_i,和t_i Output * 第1..N-1行: 第i行包含一个数 ...

7. BZOJ 1576: [Usaco2009 Jan]安全路经Travel

日常自闭半小时后看题解,太弱了qwq. 感觉这道题还是比较难的,解法十分巧妙,不容易想到. 首先题目说了起点到每个点的最短路都是唯一的,那么对这个图求最短路图必定是一棵树,而且这棵树是唯一的. 那么我 ...

8. 【BZOJ】1576 [Usaco2009 Jan]安全路经Travel

[算法]最短路树+(树链剖分+线段树)||最短路树+并查集 [题解] 两种方法的思想是一样的,首先题目限制了最短路树唯一. 那么建出最短路树后,就是询问对于每个点断掉父边后重新找路径的最小值,其它路径 ...

9. 【BZOJ1576】[Usaco2009 Jan]安全路经Travel 最短路+并查集

[BZOJ1576][Usaco2009 Jan]安全路经Travel Description Input * 第一行: 两个空格分开的数, N和M * 第2..M+1行: 三个空格分开的数a_i, ...

## 随机推荐

1. 获取文件CRC和MD5

unit untCRCMD5; interface { 获取文件CRC校验码 } function GetFileCRC(const iFileName: string): String; { 获取字 ...

2. hdoj Last non-zero Digit in N! 【数论】

找规律! 求N!最后非0位的值.比方2是120的最后一个不是0的值. 输入N比較大,要大数保存. 注意到最后0的个数是与5的因数的个数相等.设f(n)为n!的最后非0位. 那么f(n)=((n%5)! ...

3. hdu 2222 Keywords Search ac自己主动机

点击打开链接题目链接 Keywords Search Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Ja ...

4. [转] EF Configuring a DbContext

本文转自:https://docs.microsoft.com/en-us/ef/core/miscellaneous/configuring-dbcontext Note This document ...

5. 关于NetworkInfo对象的isConnected()与isAvailable()

public class MainActivity extends Activity{    /** Called when the activity is first created. */   ...

6. Active Desktop--桌面字体背景被修改

怎么修改回来 步骤如下 方法一.在桌面上点击右键 -- 排列图标 -- 去掉“在桌面上锁定Web项目”上的勾. 方法二.右键点击我的电脑 -- 属性 -- 高级 -- 点击“性能”下面的“设置”按钮, ...

7. Axiom3D学习日记 1.程序配置

1.需要引用的库 Axiom Axiom.Framework Axiom.Platforms.Win32 Axiom.Plugins.FreeImageCodecs Axiom.Plugins.Par ...

8. .NET中的三种Timer的区别和用法(转)

最近正好做一个WEB中定期执行的程序,而.NET中有3个不同的定时器.所以正好研究研究.这3个定时器分别是: //1.实现按用户定义的时间间隔引发事件的计时器.此计时器最宜用于 Windows 窗 ...

9. ashx页面 “检测到有潜在危险的 Request.Form 值”的解决方法（控制单个处理程序不检测html标签）

如题: 使用web.config的configuration/location节点. 在configuration节点内新建一个location节点,注意这个节点和system.webserver那些 ...

10. mysqldump备份、还原数据库路径名含有空格的处理方法(如：Program Files)

虽然以下的方法也可以解决,不过最简单直接的,还是直接在路径前后加双引号-" ",这个方法简单有效. 首先要说明的是mysqldump.exe在哪里不重要,重要的是要处理好路径中的非 ...