bzoj 1787 [Ahoi2008]Meet 紧急集合（1832 [AHOI2008]聚会）

1787: [Ahoi2008]Meet 紧急集合

Time Limit: 20 Sec  Memory Limit: 162 MB
Submit: 1841  Solved: 857
[Submit][Status][Discuss]

Description

Input

Output

Sample Input

6 4
1 2
2 3
2 4
4 5
5 6
4 5 6
6 3 1
2 4 4
6 6 6

Sample Output

5 2
2 5
4 1
6 0

HINT

Source

Day1

 

【思路】

       Lca。

       求出三点之间的lca，会出现两个点相同的情况，集合点为另外一个点。

 

【代码】

/**************************************************************
    Problem: 1787
    User: hahalidaxin2
    Language: C++
    Result: Accepted
    Time:3716 ms
    Memory:65216 kb
****************************************************************/
 
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<iostream>
using namespace std;
  
const int maxn = 500000+10;
const int maxd = 19;
  
struct Edge{ int u,v;
};
vector<int> G[maxn];
vector<Edge> es;
  
int d[maxn];
int p[maxn][maxd];
  
void addedge(int u,int v) {
    es.push_back((Edge){u,v});
    int m=es.size(); G[u].push_back(m-1);
}
void dfs(int u,int fa) { 
    for(int i=1;i<=maxd;i++) {       //构造倍增数组 
        if(d[u]<(1<<i)) break;
        p[u][i]=p[p[u][i-1]][i-1];
    }
    for(int i=0;i<G[u].size();i++) {
        Edge e=es[G[u][i]]; int v=e.v;
        if(v!=fa)
            d[v]=d[u]+1 , p[v][0]=u , dfs(v,u);
    }
}
int lca(int u,int v) {
    if(d[v]>d[u]) swap(u,v);
    int dep=d[u]-d[v];
    for(int i=0;i<maxd;i++)
        if((1<<i)&dep) u=p[u][i]; 
    if(u==v) return u;
    for(int i=maxd-1;i>=0;i--)
        if(p[u][i]!=p[v][i])
            u=p[u][i] , v=p[v][i];
    return p[u][0];
}
int dist(int x,int y) { return d[x]+d[y]-(d[lca(x,y)]<<1); }
int n,m;
  
void read(int& x) {
    char c=getchar(); 
    while(!isdigit(c)) c=getchar();
    x=0;
    while(isdigit(c))
        x=x*10+c-'0' , c=getchar();
}
int main() {
    read(n),read(m);
    int u,v;
    for(int i=0;i<n-1;i++) {
        read(u),read(v);
        addedge(u,v); addedge(v,u);
    }
    dfs(n>>1,-1);
    int a,b,c,lab,lac,lbc,s;
    while(m--) {
        read(a),read(b),read(c);
        lab=lca(a,b),lac=lca(a,c),lbc=lca(b,c);
        if(lab==lac) s=lbc;
        else if(lab==lbc) s=lac;
        else s=lab;
        printf("%d %d\n",s,dist(a,s)+dist(b,s)+dist(c,s));
    }
    return 0;
}