值得一做》一道类似于货车运输的题目（BZOJ3732）（easy+）

 1 #include<stdio.h>
 2 #include<algorithm>
 3 using namespace std;
 4 struct shit{
 5     int aim,from,lon,next;
 6 }e[20000],s[61000];
 7 struct ass{
 8     int fat,mx;
 9 }f[30][20000];
10 int point,dep[20000],fa[20000],fff[20000],cnt,head[20000],n,m,k;
11 bool vis[20000];
12 int LCA_(int x,int y)
13 {
14     int mather_fucker=0;
15     if(dep[x]>dep[y])swap(x,y);
16     for(int i=14;~i;--i)
17     if(dep[f[i][y].fat]>=dep[x])
18     {
19         mather_fucker=max(mather_fucker,f[i][y].mx);
20         y=f[i][y].fat;
21     }
22     if(y==x)return mather_fucker;
23     for(int i=15;~i;--i)
24     if(f[i][y].fat!=f[i][x].fat){
25         mather_fucker=max(mather_fucker,max(f[i][y].mx,f[i][x].mx));
26         y=f[i][y].fat;
27         x=f[i][x].fat;
28     }
29     mather_fucker=max(mather_fucker,max(f[0][y].mx,f[0][x].mx));
30     return mather_fucker;
31 }
32 void fuck(int x,int y,int z)
33 {
34     e[++point].aim=y;e[point].from=x;e[point].lon=z;e[point].next=head[x];head[x]=point;
35     e[++point].aim=x;e[point].from=y;e[point].lon=z;e[point].next=head[y];head[y]=point;
36 }
37 void get(int w,int x,int y,int z)
38 {
39     s[w].aim=y,s[w].from=x,s[w].lon=z;
40 }
41 int find(int x)
42 {
43     return fa[x]==x?x:fa[x]=find(fa[x]);
44 }
45 bool cmp(shit x,shit y)
46 {
47     if(x.lon<y.lon)return true;
48     else return false;
49 }
50 void work(int x,int d)
51 {
52     dep[x]=d;
53     vis[x]=true;
54     for(int i=head[x];i;i=e[i].next)
55     {
56         int u=e[i].aim;
57         if(vis[u])continue;
58         f[0][u].fat=x;
59         f[0][u].mx=e[i].lon;
60         work(u,d+1);
61     }
62 }
63 int main()
64 {
65     int a,b,c;
66     scanf("%d%d%d",&n,&m,&k);
67     for(int i=1;i<=m;i++)
68     {
69         scanf("%d%d%d",&a,&b,&c);
70         get(i,a,b,c);
71     }
72     for(int i=1;i<=n;i++)fa[i]=i;
73     sort(s+1,s+m+1,cmp);
74     for(int i=1;i<=m;++i)
75     {
76         a=find(s[i].aim);
77         b=find(s[i].from);
78         if(a==b)continue;
79         fa[a]=b;
80         fuck(a,b,s[i].lon);
81         ++cnt;
82         if(cnt==n-1)break;
83     }
84     f[0][n/2].fat=n/2,f[0][n/2].mx=0,dep[n/2]=1,vis[n/2]=true;
85     work(n/2,1);
86     for(int i=1;i<=15;i++)
87         for(int j=1;j<=n;j++)
88         {
89             f[i][j].fat=f[i-1][f[i-1][j].fat].fat;
90             f[i][j].mx=max(f[i-1][j].mx,f[i-1][f[i-1][j].fat].mx);
91         }
92     for(int i=1;i<=k;++i)
93     {
94         scanf("%d%d",&a,&b);
95         printf("%d\n",LCA_(a,b));
96     }
97     return 0;
98 }