luogu P4281 [AHOI2008]紧急集合 / 聚会

题面传送门
其实这道弱省省选还蛮简单的。
稍微思考一下我们就可以得出一个结论:集合点肯定在三个点中某两个点的\(lca\)上。
因为无论这个点朝哪个方向移都是靠近一个点远离两个点肯定不会更优。
所以我们可以先跑出三个点的\(lca\)然后再分别算答案输出。
但是这样要跑\(6\)遍\(lca\)，在\(qzezoj\)上过不去淦。
考虑优化。
然后我们又思考一下就可以优化一下结论:与其他两个不一样的\(lca\)为最优点。
因为其他两个点肯定高于原来最高的一个点，而这个点更靠近下面。又和更低的两个点接近。所以现在只要\(4\)次\(lca\)就好了。
别忘了加\(fread\)
代码实现:

#include<cstdio>
#include<cstring>
using namespace std;
int a[500039],fa[500039][25],n,m,k,lg[500039],d[500039],f[500039];
struct yyy{int to,z;};
struct ljb{
	yyy f[1000039];
	int head,h[500039];
	inline void add(int x,int y){
		f[++head]=(yyy){y,h[x]};
		h[x]=head;
	}
}s;
inline void dfs1(register int x,register int last){
	fa[x][0]=last;
	d[x]=d[last]+1;
	for(register int i=1;i<lg[d[x]];i++) fa[x][i]=fa[fa[x][i-1]][i-1];
	int cur=s.h[x];
	yyy tmp;
	while(cur!=-1){
		tmp=s.f[cur];
		if(tmp.to!=last) dfs1(tmp.to,x);
		cur=tmp.z;
	}
}
inline void swap(register int &x,register int &y){x^=y,y^=x,x^=y;}
inline int lca(register int x,register int y){
	if(d[x]<d[y]) swap(x,y);
	while(d[x]>d[y]) x=fa[x][lg[d[x]-d[y]]-1];
	if(x==y) return x;
	for(register int i=lg[d[x]]-1;i>=0;i--) if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
inline int dist(register int x,register int y){return d[x]+d[y]-d[lca(x,y)]*2;}
static char buf[1000000],*p1=buf,*p2=buf;
#define getchar() p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++
inline void read(register int &x){
	x=0;char s=getchar();
	while(s<'0'||s>'9') s=getchar();
	while(s>='0'&&s<='9') x=(x<<3)+(x<<1)+(s^48),s=getchar();
}
inline void print(int x){
	if(x>9) print(x/10);
	putchar(x%10+48);
}
int main(){
//	freopen("1.in","r",stdin);
//	freopen("1.out","w",stdout);
	register int i,j,x,y,z,ans1,ans2,ans3,lcax,lcay,lcaz;
	memset(s.h,-1,sizeof(s.h));
	read(n);read(m);
	for(i=1;i<=n;i++) lg[i]=lg[i-1]+(1<<lg[i-1]==i);
	for(i=1;i<n;i++) read(x),read(y),s.add(x,y),s.add(y,x);
	dfs1(1,0);
	for(i=1;i<=m;i++){
		read(x),read(y),read(z);
		lcax=lca(x,y);
		lcay=lca(y,z);
		lcaz=lca(x,z);
		if(lcay==lcaz) print(lcax),putchar(' '),print(d[z]+d[x]+d[y]-d[lcax]-d[lca(lcax,z)]*2),putchar('\n');
		else if(lcax==lcaz) print(lcay),putchar(' '),print(d[x]+d[y]+d[z]-d[lcay]-d[lca(x,lcay)]*2),putchar('\n');
		else print(lcaz),putchar(' '),print(d[y]+d[x]+d[z]-d[lcaz]-d[lca(y,lcaz)]*2),putchar('\n');
	}
}