码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ1787: [Ahoi2008]Meet 紧急集合

时间:2017-12-28 23:24:24      阅读:174      评论:0      收藏:0      [点我收藏+]

标签:class   ahoi2008   gpo   bzoj   printf   大数据   post   blog   做了   

【传送门:BZOJ1787


简要题意:

  给出有n个点的图,n-1条无向边,保证任意两点之间能互相到达,每条边的权值为1,给出m个询问,每个询问输入x,y,z,求出一个点使得三个点到这个点的距离和最短,输出这个点的编号和最短距离和


题解:

  终于算是做了一道较大数据结构的题了

  n个点,n-1条边其实就是一棵树,又因为每条边的权值为1,相当于深度差

  那么就转化为求三个点的LCA

  一开始看到,三个点!!LCA???身为蒟蒻的我懵逼了

  结果画了一下样例的图,发现,哎,好像三个点的LCA必定是其中两个点的LCA

  又手动搞了一下其他数据,发现三个点两两之间的LCA中必定有两个相同,而且其中不同的那个就是三个点的LCA!!!!

  惊呆了,马上码,1A!!


参考代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
struct node
{
    int x,y,next;
}a[1100000];int len,last[510000];
int f[510000][25],dep[510000];
void ins(int x,int y)
{
    len++;
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
int bin[25];
void dfs(int x,int fa)
{
    dep[x]=dep[fa]+1;
    f[x][0]=fa;
    for(int i=1;bin[i]<=dep[x];i++) f[x][i]=f[f[x][i-1]][i-1];
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y==fa) continue;
        dfs(y,x);
    }
}
int LCA(int x, int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=20;i>=0;i--) if(dep[x]-dep[y]>=bin[i]) x=f[x][i];
    if(x==y) return x;
    for(int i=20;i>=0;i--) if(dep[x]>=(1<<i)&&f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}
int n,m;
int dis(int x,int y)
{
    int lca=LCA(x,y);
    return dep[x]+dep[y]-2*dep[lca];
}
int main()
{
    scanf("%d%d",&n,&m);
    bin[0]=1;
    for(int i=1;i<=20;i++) bin[i]=bin[i-1]<<1;
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<n;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        ins(x,y);ins(y,x);
    }
    dfs(1,0);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        int l1=LCA(x,y),l2=LCA(x,z),l3=LCA(y,z);
        if(l1==l2) printf("%d %d\n",l3,dis(x,l3)+dis(y,l3)+dis(z,l3));
        else if(l1==l3) printf("%d %d\n",l2,dis(x,l2)+dis(y,l2)+dis(z,l2));
        else if(l2==l3) printf("%d %d\n",l1,dis(x,l1)+dis(y,l1)+dis(z,l1));
    }
    return 0;
}

 

BZOJ1787: [Ahoi2008]Meet 紧急集合

标签:class   ahoi2008   gpo   bzoj   printf   大数据   post   blog   做了   

原文地址:https://www.cnblogs.com/Never-mind/p/8137787.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!