POJ 2631 Roads in the North（求树的直径，两次遍历 or 树DP）

题目链接：http://poj.org/problem?id=2631

Description

Building and maintaining roads among communities in the far North is an expensive business. With this in mind, the roads are build such that there is only one route from a village to a village that does not pass through some other village twice. 
Given is an area in the far North comprising a number of villages and roads among them such that any village can be reached by road from any other village. Your job is to find the road distance between the two most remote villages in the area. 

The area has up to 10,000 villages connected by road segments. The villages are numbered from 1. 

Input

Input to the problem is a sequence of lines, each containing three positive integers: the number of a village, the number of a different village, and the length of the road segment connecting the villages in kilometers. All road segments are two-way.

Output

You are to output a single integer: the road distance between the two most remote villages in the area.

 

题目大意：给你一棵树，求树上最远的两点距离为多少。（题目描述比较奇葩，你可以认为点数等于边数+1，所有运算不会超过2^31-1）

思路：从任意一点开始，找到离这个点最远的点x。再从x开始做单源最短路，离x最远的点和x就是这棵树的最远点对。（当然基于树的分治也是可以做的）

证明（参考DISCUSS）：

设最长链是MN->已知[1]

设由A开始DFS得到最长路为AB->已知[2]

结论[1] MN与AB有公共点.否则MN<AM+AN<=AM+AB=BM 与已知[1]矛盾

结论[2] B是最长链的一个端点.否则由结论[1] 设K是AB上距B最近且在MN上的点 则MN=MK+KN=MK+AN-AK<=MK+AB-AK=MK+BK=BM 当取等号时MB与MN等长 符合结论[2] 否则与已知[1]矛盾  [这里假定了A不在NK上.若A在NK上 只须将上面式子中MN交换位置即可 不影响结论]

结论[3] 从B开始DFS得到的最长路径是一条最长链.由结论[2].B是最长链的一端

至此证毕

 

代码（0MS）：

 1 #include <iostream>

 2 #include <cstring>

 3 #include <cstdio>

 4 #include <algorithm>

 5 #include <queue>

 6 using namespace std;

 7 

 8 const int MAXN = 10010;

 9 const int MAXE = 20010;

10 

11 int dis[MAXN], head[MAXN];

12 int to[MAXE], next[MAXE], cost[MAXE];

13 int n, ecnt;

14 

15 void init() {

16     memset(head, -1, sizeof(head));

17     ecnt = 0;

18 }

19 

20 void add_edge(int u, int v, int c) {

21     to[ecnt] = v; cost[ecnt] = c; next[ecnt] = head[u]; head[u] = ecnt++;

22     to[ecnt] = u; cost[ecnt] = c; next[ecnt] = head[v]; head[v] = ecnt++;

23 }

24 

25 int bfs(int st) {

26     memset(dis, 0x3f, sizeof(dis));

27     queue<int> que; que.push(st);

28     dis[st] = 0;

29     while(!que.empty()) {

30         int u = que.front(); que.pop();

31         for(int p = head[u]; ~p; p = next[p]) {

32             int &v = to[p];

33             if(dis[v] > dis[u] + cost[p]) {

34                 dis[v] = dis[u] + cost[p];

35                 que.push(v);

36             }

37         }

38     }

39     int ed = st;

40     for(int i = 1; i <= n; ++i)

41         if(dis[i] > dis[ed]) ed = i;

42     return ed;

43 }

44 

45 int main() {

46     init();

47     n = 1;

48     int u, v, c;

49     while(scanf("%d%d%d", &u, &v, &c) != EOF) {

50         add_edge(u, v, c);

51         ++n;

52     }

53     u = bfs(1);

54     v = bfs(u);

55     printf("%d\n", dis[v]);

56 }

View Code

————————————————————————————————时间的分割线————————————————————————————————————————————————

后记（2014-7-26）：

当然这题还可以树DP，之前觉得麻烦用了上面的解法。现在有一个比较高大上的写法，边权是负数的时候也适用（据说上面的做法边权有负数就没用了，我懒得验证了有兴趣去验证一下吧……）。

思路很简单，树的直径必然是树上某一个点开始往下的最长链和次长链之和。而实现的时候只要保留最长链的大小就可以了，比较简单不讲了看代码就好。

 

代码（16MS）：

 1 #include <iostream>

 2 #include <cstring>

 3 #include <cstdio>

 4 #include <algorithm>

 5 using namespace std;

 6 

 7 const int MAXV = 10010;

 8 const int MAXE = 20010;

 9 

10 int head[MAXV];

11 int to[MAXE], next[MAXE], cost[MAXE];

12 int n, ecnt;

13 

14 void init() {

15     memset(head, -1, sizeof(head));

16     ecnt = 0;

17 }

18 

19 void add_edge(int u, int v, int c) {

20     to[ecnt] = v; cost[ecnt] = c; next[ecnt] = head[u]; head[u] = ecnt++;

21     to[ecnt] = u; cost[ecnt] = c; next[ecnt] = head[v]; head[v] = ecnt++;

22 }

23 

24 int dfs(int u, int f, int &ans) {

25     int maxdep = 0;

26     for(int p = head[u]; ~p; p = next[p]) {

27         int &v = to[p];

28         if(v == f) continue;

29         int tmp = dfs(v, u, ans) + cost[p];

30         ans = max(ans, maxdep + tmp);

31         maxdep = max(tmp, maxdep);

32     }

33     return maxdep;

34 }

35 

36 int main() {

37     init();

38     n = 1;

39     int u, v, c;

40     while(scanf("%d%d%d", &u, &v, &c) != EOF) {

41         add_edge(u, v, c);

42         ++n;

43     }

44     int ans = 0;

45     dfs(1, 0, ans);

46     printf("%d\n", ans);

47 }

View Code