• [AHOI2008]紧急集合 / 聚会


    题目描述

    欢乐岛上有个非常好玩的游戏,叫做“紧急集合”。在岛上分散有N个等待点,有N-1条道路连接着它们,每一条道路都连接某两个等待点,且通过这些道路可以走遍所有的等待点,通过道路从一个点到另一个点要花费一个游戏币。

    参加游戏的人三人一组,开始的时候,所有人员均任意分散在各个等待点上(每个点同时允许多个人等待),每个人均带有足够多的游戏币(用于支付使用道路的花费)、地图(标明等待点之间道路连接的情况)以及对话机(用于和同组的成员联系)。当集合号吹响后,每组成员之间迅速联系,了解到自己组所有成员所在的等待点后,迅速在N个等待点中确定一个集结点,组内所有成员将在该集合点集合,集合所用花费最少的组将是游戏的赢家。

    小可可和他的朋友邀请你一起参加这个游戏,由你来选择集合点,聪明的你能够完成这个任务,帮助小可可赢得游戏吗?

    输入输出格式

    输入格式:

    第一行两个正整数N和M(N<=500000,M<=500000),之间用一个空格隔开。分别表示等待点的个数(等待点也从1到N进行编号)和获奖所需要完成集合的次数。 随后有N-1行,每行用两个正整数A和B,之间用一个空格隔开,表示编号为A和编号为B的等待点之间有一条路。 接着还有M行,每行用三个正整数表示某次集合前小可可、小可可的朋友以及你所在等待点的编号。

    输出格式:

    一共有M行,每行两个数P,C,用一个空格隔开。其中第i行表示第i次集合点选择在编号为P的等待点,集合总共的花费是C个游戏币。

    输入输出样例

    输入样例#1: 复制
    6 4  
    1 2  
    2 3  
    2 4 
    4 5
    5 6
    4 5 6
    6 3 1
    2 4 4 
    6 6 6
    输出样例#1: 复制
    5 2
    2 5
    4 1
    6 0
    
    
    

    说明

    提示:

    40%的数据中N<=2000,M<=2000
    100%的数据中,N<=500000,M<=500000

    分析:

    首先考虑只有两个点的情况,很明显树上两点A,B之间距离为$ dep[a]+dep[b]-dep[LCA(a,b)] $ (dep表示该点到根节点的距离(即深度))

    那么这给我们的启示是:对于三个点能不能用也用LCA求解

    接下来我们就不难发现,对于树上三点a,b,c,它们两两之间都有一个LCA,一共就有三个,这三个中有两个是同一个且深度较小

    证明:我们假设$ x=LCA(a,b),y=LCA(a,c),且x>y $

    则对于$ z=LCA(b,c) $,首先存在路径$b->x,a->x->y$,则有路径$b->x->y$,又因为$c->y$,所以$ LCA(b,c)=y=LCA(a,c) $

    接下来我们先找集合点,假设集合点存在于那个较深的LCA(即上面的x)

    那么我们不难发现,不管是将x上移还是下移,都会使一个点的路径减少$w$,但另两个点的路径都会增加$w$

    因此将点放在那个较深的LCA最优

    关于计算路径,我们用上面所设的点表示出来是

    $$ ans=dep[a]+dep[b]-2*dep[x]+dep[x]+dep[c]-2*dep[y]$$

    $$=dep[a]+dep[b]+dep[c]-dep[x]+2*dep[y]$$

    $$=dep[a]+dep[b]+dep[c]-dep[x]-dep[y]-dep[z](y==z)$$

    我们发现不管a,b,c如何变化,由于x,y,z的值不会变化,所以ans的值也不会变化,直接代入求值即可

    #include<iostream>
    #include<string>
    #include<string.h>
    #include<stdio.h>
    #include<algorithm>
    #include<vector>
    #include<queue>
    #include<map>
    using namespace std;
    struct node{
        int to,nxt;
    }sq[1001000];
    int fa[500500][25],n,m,dep[500500],head[500500],all=0,low[500500];
    
    void add(int u,int v)
    {
        all++;sq[all].to=v;sq[all].nxt=head[u];head[u]=all;
        all++;sq[all].to=u;sq[all].nxt=head[v];head[v]=all;
    }
    
    void dfs(int u,int fu)
    {
        dep[u]=dep[fu]+1;
        fa[u][0]=fu;
        int i;
        for (i=1;i<=low[dep[u]];i++) fa[u][i]=fa[fa[u][i-1]][i-1];
        for (i=head[u];i;i=sq[i].nxt)
        {
            int v=sq[i].to;
            if (v!=fu) dfs(v,u);
        }
    }
    
    int LCA(int u,int v)
    {
        if (dep[u]<dep[v]) swap(u,v);
        while (dep[u]>dep[v]) u=fa[u][low[dep[u]-dep[v]]];
        if (u==v) return u;
        int i;
        for (i=low[dep[u]];i>=0;i--)
            if (fa[u][i]!=fa[v][i]) {u=fa[u][i];v=fa[v][i];}
        return fa[u][0];
    }
    
    int main()
    {
        scanf("%d%d",&n,&m);
        int i;
        memset(head,0,sizeof(head));
        for (i=1;i<n;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            add(u,v);
        }
        low[1]=0;
        for (i=2;i<=n;i++) low[i]=low[i>>1]+1;
        dep[0]=0;
        dfs(1,0);
        for (i=1;i<=m;i++)
        {
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            int x=LCA(a,b),y=LCA(a,c),z=LCA(b,c);
            int ansp,ansd;
            if (x==y) ansp=z;
            else if (x==z) ansp=y;
            else if (y==z) ansp=x;
            ansd=dep[a]+dep[b]+dep[c]-dep[x]-dep[y]-dep[z];
            printf("%d %d
    ",ansp,ansd);
        }
        return 0;
    }
  • 相关阅读:
    git 问题集
    es7集群安装配置及常用命令
    idea maven项目操作kafka--生产者和消费者
    linux中kafka集群搭建及常用命令
    xampp 支持php版本,支持php5.2的最后一版xampp——xampp-1.7.1
    TCP和UDP的区别及各自优缺点区别
    Qt容器(QHash/QMap等)基本学习记录
    linux系统剪切
    简单网络IP探索
    C++(Qt)线程与锁
  • 原文地址:https://www.cnblogs.com/encodetalker/p/9741671.html
Copyright © 2020-2023  润新知