zoukankan      html  css  js  c++  java
  • bzoj 1787: [Ahoi2008]Meet 紧急集合

    2016-06-22

    这个题就是找树中三个点到哪个点距离和最短,很明显是LCA,找这三个点两两之间的LCA,画个图看一下就行,很简单。

     1 #include<cstdio>
     2 #include<iostream>
     3 #include<cstdlib>
     4 #include<cstring>
     5 #include<algorithm>
     6 #include<queue>
     7 #define ll long long
     8 #define M 1000009
     9 using namespace std;
    10 ll read()
    11 {
    12     char ch=getchar();
    13     ll x=0,f=1;
    14     for(;ch<'0'||ch>'9';ch=getchar())
    15       if(ch=='-')
    16         f=-1;
    17     for(;ch>='0'&&ch<='9';ch=getchar())
    18       x=x*10+ch-'0';
    19     return x*f;
    20 }
    21 int n,m,cnt,h[M],head[M],next[M],u[M],lc[M][22];
    22 void jia(int a1,int a2)
    23 {
    24     next[++cnt]=head[a1];
    25     head[a1]=cnt;
    26     u[cnt]=a2;
    27 }
    28 void dfs(int x)
    29 {
    30     for(int i=1;(1<<i)<=h[x];i++)
    31       lc[x][i]=lc[lc[x][i-1]][i-1];
    32     for(int i=head[x];i;i=next[i])
    33       if(u[i]!=lc[x][0])
    34       {
    35           h[u[i]]=h[x]+1;
    36           lc[u[i]][0]=x;
    37           dfs(u[i]);
    38       }
    39 }
    40 int lca(int a1,int a2)
    41 {
    42     if(h[a1]<h[a2])
    43       swap(a1,a2);
    44     int a3=h[a1]-h[a2];
    45     for(int i=0;i<=20;i++)
    46       if(a3&(1<<i))
    47         a1=lc[a1][i];
    48     for(int i=20;i>=0;i--)
    49       if(lc[a1][i]!=lc[a2][i])
    50         {
    51             a1=lc[a1][i];
    52             a2=lc[a2][i];
    53         }
    54     if(a1==a2)
    55        return a1;
    56     return lc[a1][0];
    57 }
    58 int main()
    59 {
    60     n=read();
    61     m=read();
    62     for(int i=1;i<n;i++)
    63       {
    64         int a1=read(),a2=read();
    65         jia(a1,a2);
    66         jia(a2,a1);
    67       }
    68     dfs(1);
    69     for(int i=1;i<=m;i++)
    70       {
    71           int a1=read(),a2=read(),a3=read(),a4;
    72           int t1=lca(a1,a2),t2=lca(a1,a3),t3=lca(a2,a3);
    73         if(h[t2]<h[t1])
    74           printf("%d %d
    ",t1,h[a1]+h[a2]-h[t1]+h[a3]-2*h[t2]);
    75         else
    76           {
    77             if(h[t2]>h[t3])
    78               a4=t2;
    79             else
    80               a4=t3;
    81             printf("%d %d
    ",a4,h[a1]+h[a2]-2*h[t1]+h[a3]-h[a4]);    
    82           }
    83       }
    84     return 0;
    85 }
  • 相关阅读:
    如何在一台主机上管理自己的多个git repository
    Linux deepin 中Jetbrain Idea等软件中文显示异常
    使用openmp进行并行编程
    使用pthread进行编程
    使用openmp进行共享内存编程
    并行与分布式计算之并行硬件与并行软件
    算法笔记刷题5(PAT A1025)
    算法笔记刷题4(PAT B1009)
    算法笔记刷题3(codeup 5901)
    算法笔记刷题2(codeup 1928)
  • 原文地址:https://www.cnblogs.com/xiw5/p/5605830.html
Copyright © 2011-2022 走看看