商务旅行

Description

　　某首都城市的商人要经常到各城镇去做生意，他们按自己的路线去做，目的是为了更好的节约时间。
　　假设有N个城镇，首都编号为1，商人从首都出发，其他各城镇之间都有道路连接，任意两个城镇之间如果有直连道路，在他们之间行驶需要花费单位时间。该国公路网络发达，从首都出发能到达任意一个城镇，并且公路网络不会存在环。
　　你的任务是帮助该商人计算一下他的最短旅行时间。

Input

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

Output

　　在输出文件kom.out中输出该商人旅行的最短时间。

Sample Input

5
1 2
1 5
3 5
4 5
4
1
3
2
5
Sample Output

7

题解

lca的模板题，先跑一次记录没给点的深度，（s,t）的最短距离是d[s]+d[t]-2*d[fa(s,t)]~~o(>_<)o ~~，第一次做，多打几遍就会了。

代码

type  ab=^node;    node=record    a,b:longint;    c:ab;end;var  n,m,tot,sum:longint;  p:array[0..30000]of ab;  v:array[0..30000]of boolean;  f,f1:array[0..30000,0..200]of longint;  d:array[0..30000]of longint;procedure put(x,y:longint);var  i:ab;begin  new(i);  i^.a:=y;  i^.c:=p[x];  p[x]:=i;end;procedure fq(x,y,z:longint);var  i,j:longint;begin  f[y,0]:=x;  f1[y,0]:=z;  i:=x; j:=1;  while f[i,j-1]>0 do    begin      f[y,j]:=f[i,j-1];      f1[y,j]:=f1[y,j-1]+f1[i,j-1];      i:=f[i,j-1];      inc(j);    end;end;procedure find(x,y:longint);var  i:ab;begin  v[x]:=true;  d[x]:=y;  i:=p[x];  while i<>nil do    begin      if v[i^.a]=false then        begin          fq(x,i^.a,1);          find(i^.a,y+1);        end;      i:=i^.c;    end;end;function fa(x,y:longint):longint;var  i,j:longint;begin  if d[x]<d[y] then    begin      i:=x;x:=y;y:=i;    end;  i:=d[x]-d[y];  j:=0;  while i>0 do    begin      if i mod 2=1 then        begin          sum:=f1[x,j]+sum;          x:=f[x,j];        end;       i:=i div 2;       j:=j+1;    end;  if x=y then exit(sum);  j:=0;  while f[x,0]<>f[y,0] do    begin       while f[x,j]=f[y,j] do dec(j);       sum:=sum+f1[x,j];       x:=f[x,j];       sum:=sum+f1[y,j];       y:=f[y,j];       j:=j+1;    end;  sum:=sum+f1[x,0];  sum:=sum+f1[y,0];  exit(sum);end;var  i,s,t:longint;begin  readln(n);  for i:=1 to n-1 do    begin     readln(s,t);     put(s,t);     put(t,s);    end;  find(1,1);  readln(m);  readln(s);  for i:=2 to m do    begin      sum:=0;      readln(t);      tot:=tot+fa(s,t);      s:=t;    end;  writeln(tot);end.