[2018.12.26]BZOJ1787 [Ahoi2008]Meet 紧急集合

一个结论:

树上任意三点两两之间的LCA必然有两个是相同的。

至于为什么。。。找个规律就好了干嘛那么烦,暂时没找到好的证明。

可以发现答案必然在这两个LCA之一。

发现只出现一次的那个更优。

大致就是因为如果走到重复的那个LCA,它和不重复的LCA之间的路径被多走了一次。而从后者出发的话路径是不会重叠的。画个图会更清楚。

由于之前没看出结论,试图树剖强上,所以懒得重新写,求LCA用了树剖。

当然树剖很慢(过这道题总用时1432ms,最大点用时632ms),当然也可能是蒟蒻选手自带大常数

code:

#include<bits/stdc++.h>
using namespace std;
int n,m,u,v,w,dep[500010],f[500010],sz[500010],dfn[500010],tp[500010],siz,l1,l2,l3;
vector<int>e[500010];
void dfs1(int x){
	sz[x]=1;
	for(int i=0;i<e[x].size();i++){
		if(f[x]!=e[x][i])f[e[x][i]]=x,dep[e[x][i]]=dep[x]+1,dfs1(e[x][i]),sz[x]+=sz[e[x][i]];
	}
}
void dfs2(int x){
	dfn[x]=++siz;
	int mx=0,d=0;
	for(int i=0;i<e[x].size();i++){
		if(f[x]!=e[x][i]&&sz[e[x][i]]>mx)mx=sz[e[x][i]],d=e[x][i];
	}
	if(d)tp[d]=tp[x],dfs2(d);
	for(int i=0;i<e[x].size();i++){
		if(f[x]!=e[x][i]&&e[x][i]!=d)tp[e[x][i]]=e[x][i],dfs2(e[x][i]);
	}
}
int FindLCA(int x,int y){//找两点LCA
	while(tp[x]!=tp[y]){
		if(dep[tp[x]]<dep[tp[y]])swap(x,y);
		x=f[tp[x]];
	}
	return dfn[x]<dfn[y]?x:y;
}
int Getdis(int x,int y){//求两点距离(当然可以和上面那个和在一起写)
	int dis=0;
	while(tp[x]!=tp[y]){
		if(dep[tp[x]]<dep[tp[y]])swap(x,y);
		dis+=dfn[x]-dfn[tp[x]]+1;
		x=f[tp[x]];
	}
	dis+=abs(dfn[x]-dfn[y]);
	return dis;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<n;i++)scanf("%d%d",&u,&v),e[u].push_back(v),e[v].push_back(u);
	dep[1]=1;
	dfs1(1);
	tp[1]=1;
	dfs2(1);
	while(m--){
		scanf("%d%d%d",&u,&v,&w);
		l1=FindLCA(u,v);
		l2=FindLCA(v,w);
		l3=FindLCA(w,u);
		if(l1==l2)printf("%d %d\n",l3,Getdis(l3,u)+Getdis(l3,v)+Getdis(l3,w));
		else if(l2==l3)	printf("%d %d\n",l1,Getdis(l1,u)+Getdis(l1,v)+Getdis(l1,w));
		else printf("%d %d\n",l2,Getdis(l2,u)+Getdis(l2,v)+Getdis(l2,w));
	}
	return 0;
}
posted @ 2019-03-15 15:37  xryjr233  阅读(103)  评论(0编辑  收藏  举报