bzoj 1787 && bzoj 1832: [Ahoi2008]Meet 紧急集合（倍增LCA）

1832: [AHOI2008]聚会

Time Limit: 10 Sec   Memory Limit: 64 MB
Submit: 1539   Solved: 663
[ Submit][ Status][ Discuss]

Description

Y岛风景美丽宜人，气候温和，物产丰富。Y岛上有N个城市，有N-1条城市间的道路连接着它们。每一条道路都连接某两个城市。幸运的是，小可可通过这些道路可以走遍Y岛的所有城市。神奇的是，乘车经过每条道路所需要的费用都是一样的。小可可，小卡卡和小YY经常想聚会，每次聚会，他们都会选择一个城市，使得3个人到达这个城市的总费用最小。 由于他们计划中还会有很多次聚会，每次都选择一个地点是很烦人的事情，所以他们决定把这件事情交给你来完成。他们会提供给你地图以及若干次聚会前他们所处的位置，希望你为他们的每一次聚会选择一个合适的地点。

Input

第一行两个正整数，N和M。分别表示城市个数和聚会次数。后面有N-1行，每行用两个正整数A和B表示编号为A和编号为B的城市之间有一条路。城市的编号是从1到N的。再后面有M行，每行用三个正整数表示一次聚会的情况：小可可所在的城市编号，小卡卡所在的城市编号以及小YY所在的城市编号。

Output

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

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

答案一定是其中两个点的LCA，然后第三个点走到这个LCA上

暴力三种情况就好了

#include<stdio.h>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> G[500005];
int fa[500005][24], d[500005];
void Sech(int u)
{
	int i, j, v;
	for(i=0;i<G[u].size();i++)
	{
		v = G[u][i];
		if(v==fa[u][0])
			continue;
		fa[v][0] = u;
		for(j=0;fa[fa[v][j]][j]!=0;j++)
			fa[v][j+1] = fa[fa[v][j]][j];
		d[v] = d[u]+1;
		Sech(v);
	}
}
int LCA(int u, int v)
{
	int i;
	if(d[u]<d[v])
		swap(u, v);
	for(i=21;i>=0;i--)
	{
		if(d[fa[u][i]]>=d[v])
			u = fa[u][i];
	}
	if(u==v)
		return u;
	for(i=21;i>=0;i--)
	{
		if(fa[u][i]!=fa[v][i])
			u = fa[u][i], v = fa[v][i];
	}
	return fa[u][0];
}
int main(void)
{
	int n, m, i, x, y, z, ans, p, q, temp;
	scanf("%d%d", &n, &m);
	for(i=1;i<=n-1;i++)
	{
		scanf("%d%d", &x, &y);
		G[x].push_back(y);
		G[y].push_back(x);
	}
	d[1] = 1;
	Sech(1);
	while(m--)
	{
		scanf("%d%d%d", &x, &y, &z);
		p = LCA(x, y);
		q = LCA(z, p);
		ans = d[x]+d[y]-d[p]*2+d[z]-d[q]+d[p]-d[q], temp = p;
		p = LCA(x, z);
		q = LCA(y, p);
		if(d[x]+d[z]-d[p]*2+d[y]-d[q]+d[p]-d[q]<ans)
			ans = d[x]+d[z]-d[p]*2+d[y]-d[q]+d[p]-d[q], temp = p;
		p = LCA(y, z);
		q = LCA(x, p);
		if(d[z]+d[y]-d[p]*2+d[x]-d[q]+d[p]-d[q]<ans)
			ans = d[z]+d[y]-d[p]*2+d[x]-d[q]+d[p]-d[q], temp = p;
		printf("%d %d\n", temp, ans);
	}
	return 0;
}