BZOJ1787 Meet紧急集合

目录

• BZOJ1787 Meet紧急集合
• 题解
• ｃｏｄｅ

BZOJ1787 Meet紧急集合

题目传送门

题解

写起来很简单，但是想想又有点玄的一道题。刚开始想的是用树剖维护他们的路径什么的，后来发现自己就想错了。。我们如果要达到最优的答案，那么就是要最小化重复的路径，对于树上的任意三个点，我们总共只会有1个或者2个的LCA。如果这三个点的LCA只有一个，那么直接把这个目标点设为这个LCA即可，因为这样是不会有路径重复的，如果这三个点的LCA有两个，那么必定有两个LCA是重复的，那么取另一个就行了，因为如果选择LCA相同的那个点，则说明必定是有两条路径是重复的，这样就不是最优的方案了。至于距离，随便搞搞就行了。
双倍经验题：BZOJ1832，不过空间要稍微开少一点，不然会MLE(雾

ｃｏｄｅ

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
bool Finish_read;
template<class T>inline void read(T &x){Finish_read=0;x=0;int f=1;char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;if(ch==EOF)return;ch=getchar();}while(isdigit(ch))x=x*10+ch-'0',ch=getchar();x*=f;Finish_read=1;}
template<class T>inline void print(T x){if(x/10!=0)print(x/10);putchar(x%10+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('
');}
template<class T>inline void write(T x){if(x<0)putchar('-');x=abs(x);print(x);}
/*================Header Template==============*/
const int maxn=5e5+500;
int n,q,tot;
struct edge {
	int to,nxt;
}E[maxn<<2];
int head[maxn];
int dep[maxn],p[maxn][22],dis[maxn];
/*==================Define Area================*/
void addedge(int u,int v) {
	E[++tot].to=v;E[tot].nxt=head[u];head[u]=tot;
	E[++tot].to=u;E[tot].nxt=head[v];head[v]=tot;
}

namespace LCA {
	void Dfs(int o,int ff,int deep) {
		dep[o]=deep;p[o][0]=ff;dis[o]=dis[ff]+1;
		for(int i=head[o];~i;i=E[i].nxt) {
			int to=E[i].to;
			if(to==ff) continue;
			Dfs(to,o,deep+1);
		}
	}
	void init() {
		Dfs(1,0,1);
		for(int k=1;k<=20;k++) {
			for(int i=1;i<=n;i++) {
				p[i][k]=p[p[i][k-1]][k-1];
			}
		}
	}
	int Lca(int x,int y) {
		if(dep[x]<dep[y]) swap(x,y);
		for(int i=20;~i;i--) {
			if(dep[p[x][i]]>=dep[y]) x=p[x][i];
		}
		if(x==y) return x;
		for(int i=20;~i;i--) {
			if(p[x][i]!=p[y][i]) x=p[x][i],y=p[y][i];
		}
		return p[x][0];
	}
}
using namespace LCA;
int main() {
	memset(head,-1,sizeof head);
	read(n);read(q);
	for(int i=1,u,v;i<n;i++) {
		read(u);read(v);
		addedge(u,v);
	}
	init();
	while(q--) {
		int x,y,z;
		read(x);read(y);read(z);
		int lcaxy=Lca(x,y),lcayz=Lca(y,z),lcaxz=Lca(x,z);
		if(lcaxy==lcayz) {
			int Dis=dis[x]+dis[z]-2*dis[lcaxz];
			Dis+=dis[y]+dis[lcaxz]-2*dis[lcaxy];
			printf("%d %d
",lcaxz,Dis);
		}
		else if(lcaxy==lcaxz) {
			int Dis=dis[y]+dis[z]-2*dis[lcayz];
			Dis+=dis[x]+dis[lcayz]-2*dis[lcaxy];
			printf("%d %d
",lcayz,Dis);
		}
		else if(lcayz==lcaxz) {
			int Dis=dis[x]+dis[y]-2*dis[lcaxy];
			Dis+=dis[z]+dis[lcaxy]-2*dis[lcayz];
			printf("%d %d
",lcaxy,Dis);
		}
	}
	return 0;
}