[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;
}

转载于:https://www.cnblogs.com/xryjr233/p/BZOJ1787.html