【wikioi】1036商务旅行

题目描述  Description

某首都城市的商人要常常到各城镇去作生意，他们按本身的路线去作，目的是为了更好的节约时间。

假设有N个城镇，首都编号为1，商人从首都出发，其余各城镇之间都有道路链接，任意两个城镇之间若是有直连道路，在他们之间行驶须要花费单位时间。该国公路网络发达，从首都出发能到达任意一个城镇，而且公路网络不会存在环。

你的任务是帮助该商人计算一下他的最短旅行时间。

 

 

输入描述  Input Description

输入文件中的第一行有一个整数N，1<=n<=30 000，为城镇的数目。下面N-1行，每行由两个整数a 和b (1<=a, b<=n; a<>b)组成，表示城镇a和城镇b有公路链接。在第N+1行为一个整数M，下面的M行，每行有该商人须要顺次通过的各城镇编号。

输出描述  Output Description

    在输出文件中输出该商人旅行的最短期。

样例输入  Sample Input

5

1 2

1 5

3 5

4 5

4

1

3

2

5

样例输出  Sample Output

7

 

 

【思路】：

典型的lca，因为是顺次通过点，因此从1出发，而后依次求每两个点之间的lca，而后过程当中把deep加起来。。话说这道题被坑了很久。。从下午到晚上。。而后发现是无向边，因而要加双向边，果断推翻重写。。而后终于AC了。。

【代码】：

 1 #include<cstdio>
 2 #include<algorithm>
 3 
 4 using namespace std;
 5 
 6 struct edge{
 7     int u,v,next;
 8 }e[30000 * 2 + 5];
 9 int head[37000];
10 int deep[37000];
11 int p[37000][37];
12 int n,m,tot,k=1;
13 
14 void adde(int u,int v,int k)
15 {
16     e[k].next=head[u];
17     head[u]=k;
18     e[k].u=u;
19     e[k].v=v;
20 }
21 void dfs(int fa,int x)
22 {
23     if(deep[x]) return ;
24     deep[x]=deep[fa]+1;
25     p[x][0]=fa;
26     for(int i=1;p[p[x][i-1]][i-1];i++) {
27         p[x][i]=p[p[x][i-1]][i-1];
28     }
29     for(int i=head[x];i;i=e[i].next) {
30         dfs(x,e[i].v);
31     }
32 }
33 int lca(int u,int v)
34 {
35     if(deep[u]<deep[v]) swap(u,v);
36     if(deep[u]>deep[v]) {
37         int dva=deep[u]-deep[v];
38         for(int i=0;i<=24;i++) {
39             if(dva&(1<<i)) {
40                 u=p[u][i];
41             }
42         }
43     }
44     if(u!=v) {
45         for(int i=24;i>=0;i--) {
46             if(p[u][i]!=p[v][i]) {
47                 u=p[u][i];
48                 v=p[v][i];
49             }
50         }
51     }
52     if(u==v) return u;
53     else return p[u][0];
54 }
55 int main()
56 {
57     scanf("%d",&n);
58     for(int i=1;i<n;i++) {
59         int u,v;
60         scanf("%d%d",&u,&v);
61         adde(u,v,k++);
62         adde(v,u,k++);
63     }
64     deep[1]=1;
65     for(int i=head[1];i;i=e[i].next) {
66         dfs(1,e[i].v);
67     }
68     int a = 1;
69     scanf("%d",&m);
70     for(int i=1;i<=m;i++) {
71         int b;
72         scanf("%d",&b);
73         int Lca=lca(a,b);
74         tot+=deep[a]+deep[b]-2*deep[Lca];
75         a=b;
76     }
77     printf("%d",tot);
78     return 0;
79 }

View Code