LOJ #10136. 「一本通 4.4 练习 3」聚会

题目描述

原题来自：AHOI 2008

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

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

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

输入格式

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

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

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

输出格式

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

样例

样例输入

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%40\%40% 的数据中，1≤N,M≤2×1031\le N,M\le 2\times 10^31≤N,M≤2×103；

100%100\%100% 的数据中，1≤N,M≤5×1051\le N,M\le 5\times 10^51≤N,M≤5×105。

 

很容易看出来，求3个LCA，选深度最低的一个作为结点，计算

#include<cstdio>
#include<iostream>
using namespace std;
int read()
{
	int ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n,m,cnt,id,l1,l2,l3,ans;
const int N=1e6+5;
int lg[N],dep[N],f[N][21],he[N],to[N],nxt[N];

inline void add(int u,int v)
{
	to[++cnt]=v;
	nxt[cnt]=he[u];
	he[u]=cnt;
}

void dfs(int fa,int u)
{
	dep[u]=!fa?0:dep[fa]+1;
	f[u][0]=fa; 
	for(int i=1;i<=lg[dep[u]];i++)
		f[u][i]=f[f[u][i-1]][i-1];
	for(int e=he[u];e;e=nxt[e])
	{
		int v=to[e];
		if(v!=fa) dfs(u,v);
	}
}

int LCA(int u,int v)
{
	if(dep[u]>dep[v]) swap(u,v);
	ans+=dep[v]-dep[u];
	while(dep[u]<dep[v])
		v=f[v][lg[dep[v]-dep[u]]];
	for(int i=lg[dep[u]];i>=0;i--)
		if(f[u][i]!=f[v][i])
			ans+=1<<(i+1),
			u=f[u][i],v=f[v][i];
	if(u!=v) ans+=2,u=f[u][0];
	return u;
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<n;i++)
	{
		int u=read(),v=read();
		add(u,v),add(v,u);
	}
	lg[0]=-1;
	for(int i=1;i<=n;i++)
		lg[i]=lg[i>>1]+1;
		
	dfs(0,1);
	while(m--)
	{
		int u1=read(),u2=read(),u3=read();
		ans=0;
		l1=LCA(u1,u2);
		l2=LCA(u1,u3);
		l3=LCA(u2,u3);
		if(dep[l1]>dep[l2]) id=l1;
			else id=l2;
		if(dep[l3]>dep[id]) id=l3;
		printf("%d %d\n",id,ans>>1);
	}
	return 0;
}