LOJ1036

AHOI 2008 聚会

Y 岛风景美丽宜人，气候温和，物产丰富。Y 岛上有 N 个城市，有 N-1 条城市间的道路连接着它们。每一条道路都连接某两个城市。幸运的是，小可可通过这些道路可以走遍 Y 岛的所有城市。神奇的是，乘车经过每条道路所需要的费用都是一样的。

小可可，小卡卡和小 YY 经常想聚会，每次聚会，他们都会选择一个城市，使得三个人到达这个城市的总费用最小。

由于他们计划中还会有很多次聚会，每次都选择一个地点是很烦人的事情，所以他们决定把这件事情交给你来完成。他们会提供给你地图以及若干次聚会前他们所处的位置，希望你为他们的每一次聚会选择一个合适的地点。

输入格式

第一行两个正整数，N 和 M。分别表示城市个数和聚会次数；

后面有 N−1 行，每行用两个正整数 A 和 B 表示编号为 A 和编号为 B 的城市之间有一条路。城市的编号是从 1 到 N 的；

再后面有 M 行，每行用三个正整数表示一次聚会的情况：小可可所在的城市编号，小卡卡所在的城市编号以及小 YY 所在的城市编号。

输出格式

一共有 M 行，每行两个数 P 和 C，用一个空格隔开。表示第 ii 次聚会的地点选择在编号为 P 的城市，总共的费用是经过 C 条道路所花费的费用。

样例

样例输入

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

样例输出

5 2
2 5
4 1
6 0

数据范围与提示

40% 的数据中，1≤N,M≤2×10^3；

100% 的数据中，1≤N,M≤5×10^5。

________________________________________________________________________________________________

题目还是很有意思的，但是调试时最后两个点超时，比较难调。主要是变成习惯不好！

题目大意：树上有3个点，找一点，使得到3点的距离和最小。输出所找到的点和距离和。

从图上看，3个点的关系有两种，一，第三点在前两点的路径上；二、第三点在前两点路径的某一点分出的岔路上；

情况一，中间点就是聚会点，距离就是前两点的和；

情况二， 那个“某一点”就是聚会点，距离为三点到“某一点”的距离。

（证明可以贪心一下）

那么最终结果就是，三点中两两求LCA，3个LCA求异或就是聚会点，而三个点到根的距离减去三个LCA到根的距离就是距离和。

恶心的调试，最后两个点超时，都是1010ms

改成优化读入，作用为0；

把LCA中的循环从20开始改成log(dep[u])开始，第9个点变成了978ms，10号点没变。

于是第一个循环改成

while(dep[u]>dep[v])u=f[u][lg[dep[u]-dep[v]]];
优化不大！！！

最后把所有的LL变成的int，当然最终的距离除外！就OK了！！！
而且时间大大优化！！！
所以说，int还是比LL要快不少的。

________________________________________________________________________________________________

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn=5e5+10;
ll n,m;
void readint(ll &x)
{
    ll f=1;
    char c=getchar();
    for(;c>'9' || c<'0';c=getchar())if(c=='-')f=-f;
    x=0;
    for(;c>='0' && c<='9';c=getchar())x=x*10+c-'0';
    x*=f;
}
struct edge
{
    ll u,v,nxt;
}e[maxn<<1];
ll head[maxn],js;
ll lg[maxn];
void init()
{
    lg[0]=-1;
    for(int i=1;i<n;++i)lg[i]=lg[i>>1]+1;
}
void addage(ll u,ll v)
{
    e[++js].u=u;e[js].v=v;
    e[js].nxt=head[u];head[u]=js;
}
ll f[maxn][20],dep[maxn];
void dfs(ll u,ll fa)
{
    f[u][0]=fa;dep[u]=dep[fa]+1;
    for(ll i=1;(1<<i)<dep[u];++i)f[u][i]=f[f[u][i-1]][i-1];
    for(ll i=head[u];i;i=e[i].nxt)
    {
        ll v=e[i].v;
        if(v!=fa)dfs(v,u);
    }
}
ll lca(ll u,ll v)
{
    if(dep[u]<dep[v])swap(u,v);
    while(dep[u]>dep[v])u=f[u][lg[dep[u]-dep[v]]];
    if(u==v)return u;
    for(ll i=lg[dep[u]];i>=0;--i)if(f[u][i]!=f[v][i])u=f[u][i],v=f[v][i];
    return f[u][0];
}
int main()
{
    readint(n);readint(m);
    for(ll u,v,i=1;i<n;++i)
    {
        readint(u);readint(v);
        addage(u,v);addage(v,u);
    }
    dfs(1,0);
    init();
    while(m--)
    {
        ll a,b,c,l1,l2,l3;
        readint(a);readint(b);readint(c);
        l1=lca(a,b);l2=lca(b,c);l3=lca(c,a);
        printf("%lld %lld
",l1^l2^l3,dep[a]+dep[b]+dep[c]-dep[l1]-dep[l2]-dep[l3]);
    }
    return 0;
}