http://www.lydsy.com/JudgeOnline/problem.php?id=1602

題目描述

N頭牛（2<=n<=1000）別人被標記為1到n，在同樣被標記1到n的n塊土地上吃草，第i頭牛在第i塊牧場吃草。 這n塊土地被n-1條邊連接。 奶牛可以在邊上行走，第i條邊連接第Ai，Bi塊牧場，第i條邊的長度是Li（1<=Li<=10000）。 這些邊被安排成任意兩頭奶牛都可以通過這些邊到達的情況，所以說這是一棵樹。 這些奶牛是非常喜歡交際的，經常會去互相訪問,他們想讓你去幫助他們計算Q(1<=q<=1000)對奶牛之間的距離。

輸入

*第一行：兩個被空格隔開的整數：N和Q

*第二行到第n行：第i+1行有兩個被空格隔開的整數：AI，BI，LI

*第n+1行到n+Q行：每一行有兩個空格隔開的整數：P1，P2，表示兩頭奶牛的編號。

輸出

*第1行到第Q行：每行輸出一個數，表示那兩頭奶牛之間的距離。

樣例輸入

樣例輸出

提示

來源

資格賽

很慘的一下午：倍增死活沒過，然後改樹剖就A了~~~QAQ

先貼個沒過的請大佬找找錯

 1 #include <algorithm>
 2 #include <cstdio>
 3
 
  
 4 using namespace std;
 5  
 6 const int N(1000+15);
 7 int n,q,x,y,z;
 8 int dad[N<<1][23],deep[N<<1];
 9 int dis[N<<1];
10 int sumedge,head[N<<1];
11 struct Edge
12 {
13     int from,to,next,dis;
14     Edge(int from=0,int to=0,int next=0,int dis=0) :
15         from
 
(from),to(to),next(next),dis(dis) {}
16 }edge[N<<1];
17  
18 int add(int from,int to,int dis)
19 {
20     edge[++sumedge]=Edge(from,to,head[from],dis);
21     return head[from]=sumedge;
22 }
23  
24 void DFS(int x)
25 {
26     deep[x]=deep[dad[x][0]]+1;
27     for(int i=0;dad[x][i];i++)
28         dad[x][i+1]=dad[dad[x][i]][i];
29     for(int i=head[x];i;i=edge[i].next)
30         if(!deep[edge[i].to])
31         {
32             dad[edge[i].to][0]=x;
33             dis[edge[i].to]=dis[x]+edge[i].dis;
34             DFS(edge[i].to);
35         }
36 }
37  
38 int LCA(int x,int y)
39 {
40     if(deep[x]>deep[y]) swap(x,y);
41     for(int i=20;i>=0;i--)
42         if(deep[x]<=deep[dad[y][i]]) y=dad[y][i];
43     if(x==y) return x;
44     for(int i=20;i>=0;i--)
45         if(deep[dad[x][i]]!=deep[dad[y][i]])
46             x=dad[x][i],y=dad[y][i];
47     return dad[x][0];
48 }
49  
50 int main()
51 {
52     scanf("%d%d",&n,&q);
53     for(int i=1;i<n;i++)
54     {
55         scanf("%d%d%d",&x,&y,&z);
56         add(x,y,z); add(y,x,z);
57     }
58     DFS(1);
59     for(;q;q--)
60     {
61         scanf("%d%d",&x,&y);
62         printf("%d\n",dis[x]+dis[y]-(dis[LCA(x,y)]<<1));
63     }
64     return 0;
65 }

然後是樹剖的~

 1 #include <algorithm>
 2 #include <cstdio>
 3 #include <vector>
 4 
 5 using namespace std;
 6 
 7 const int N(1000+15);
 8 vector< pair<int,int> >vec[N];
 9 int n,q,x,y,z,dis[N];
10 int dad[N],deep[N],size[N],top[N];
11 
12 void DFS(int x)
13 {
14     size[x]=1; deep[x]=deep[dad[x]]+1;
15     for(int i=0;i<vec[x].size();i++)
16         if(dad[x]!=vec[x][i].first)
17         {
18             dad[vec[x][i].first]=x;
19             dis[vec[x][i].first]=dis[x]+vec[x][i].second;
20             DFS(vec[x][i].first);
21             size[x]+=size[vec[x][i].first];
22         }
23 }
24 
25 void DFS_(int x)
26 {
27     int t=0; if(!top[x]) top[x]=x;
28     for(int i=0;i<vec[x].size();i++)
29         if(dad[x]!=vec[x][i].first&&size[t]<size[vec[x][i].first]) t=vec[x][i].first;
30     if(t) top[t]=top[x],DFS_(t);
31     for(int i=0;i<vec[x].size();i++)
32         if(dad[x]!=vec[x][i].first&&t!=vec[x][i].first) DFS_(vec[x][i].first);
33 }
34 
35 int LCA(int x,int y)
36 {
37     while(top[x]!=top[y])
38     {
39         if(deep[top[x]]<deep[top[y]]) swap(x,y);
40         x=dad[top[x]];
41     }
42     if(deep[x]>deep[y]) swap(x,y);
43     return x;
44 }
45 
46 int main()
47 {
48     scanf("%d%d",&n,&q);
49     for(int i=1;i<n;i++)
50     {
51         scanf("%d%d%d",&x,&y,&z);
52         vec[x].push_back(make_pair(y,z));
53         vec[y].push_back(make_pair(x,z));
54     }
55     DFS(1); DFS_(1);
56     for(;q;q--)
57     {
58         scanf("%d%d",&x,&y);
59         printf("%d\n",dis[x]+dis[y]-dis[LCA(x,y)]*2);
60     }
61     return 0;
62 }

BZOJ——1602: [Usaco2008 Oct]牧場行走