[AHOI2008] BZOJ 1787 Meet紧急集合-LCA

题目链接：右转进入题目

题目大意：请自行参考题目；

做法：

其实注意到，三个点，两两搭配，会有三个LCA；可以证明，这三个LCA至少有两个是相同的；

不妨设相同的是LCA1和LCA2，可以证明，如果选择LCA3作为答案，由于边权为1，那么花费达到最小值。

更进一步的说，这个最小值就是虚树的边数。

证明方法十分简单。

先证至少有两个LCA相同。

不妨设三个点是a,b,c。

现在考虑LCA(b,c)和LCA(a,b)。显然LCA(b,c)和LCA(a,b)都在b到根的路径上。

如果LCA(b,c)==LCA(a,b)，结论成立。

如果LCA(b,c)是LCA(a,b)的祖先，则显然LCA(a,c)==LCA(b,c)，结论成立。

如果LCA(a,b)是LCA(b,c)的祖先，则显然LCA(a,c)==LCA(a,b)，结论成立。

综上至少有两个LCA相同。

再证如果LCA1=LCA2，则选LCA3作为答案花费最小。

不难发现无论如何，选择LCA3作为答案，花费都是虚树的边数。

从第一个证明不难发现，相同的那个两个LCA一定是第三个LCA的祖先（或者本身）。

无论哪种情况，显然如果不选LCA3，都会包含虚树的边数且会有重复计算。

因此选择LCA3是最优的。

一个可以简化代码的地方：

再求u,v的LCA的时候，可以顺便统计出u到v的路径长度，不放记为ai（1<=i<=3）。

那么观察可知ans=(a1+a2+a3)/2。（a1+a2+a3就是虚树边数*2）

代码如下：

//AHOI2008
//BZOJ 1787
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<vector>
#define MAXN 500010
#define MAX_LOG 20
using namespace std;
int up[MAXN][MAX_LOG],depth[MAXN];
int sumn[MAXN][MAX_LOG];
vector<int> g[MAXN];
int dfs(int rt,int fa)
{
	for(int i=1;i<MAX_LOG;i++)
		up[rt][i]=up[up[rt][i-1]][i-1],
		sumn[rt][i]=sumn[rt][i-1]+sumn[up[rt][i-1]][i-1];
	for(int i=g[rt].size()-1;i>=0;i--)
		if(g[rt][i]!=fa)
		{
			up[g[rt][i]][0]=rt;
			sumn[g[rt][i]][0]=1;
			depth[g[rt][i]]=depth[rt]+1;
			dfs(g[rt][i],rt);
		}
	return 0;
}
int LCA(int u,int v,int &lca)
{
	int ans=0;
	if(depth[u]<depth[v]) swap(u,v);
	for(int i=MAX_LOG-1;i>=0;i--)
		if(depth[up[u][i]]>=depth[v])
		{
			ans+=sumn[u][i];
			u=up[u][i];
		}
	if(u==v)
	{
		lca=u;return ans;
	}
	for(int i=MAX_LOG-1;i>=0;i--)
		if(up[u][i]!=up[v][i])
		{
			ans+=sumn[u][i]+sumn[v][i];
			u=up[u][i];v=up[v][i];
		}
	ans+=sumn[u][0]+sumn[v][0];
	lca=up[u][0];
	return ans;
}
int n,m;
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<n;i++)
	{
		int u,v;scanf("%d%d",&u,&v);
		g[u].push_back(v);
		g[v].push_back(u);
	}
	int rt;dfs(depth[rt=1]=1,0);
	while(m--)
	{
		int l1,l2,l3,a1,a2,a3,a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		a1=LCA(a,b,l1);a2=LCA(b,c,l2);a3=LCA(a,c,l3);
		if(l1==l2) printf("%d ",l3);
		else if(l1==l3) printf("%d ",l2);
		else printf("%d ",l1);
		printf("%d\n",(a1+a2+a3)/2);
	}
	return 0;
}