目录
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+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('\n');}
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\n",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\n",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\n",lcaxy,Dis);
}
}
return 0;
}
转载于:https://www.cnblogs.com/Apocrypha/p/9453125.html
相关资源：JAVA上百实例源码以及开源项目
