[bzoj1832][bzoj1787][lca]聚会 & Meet 紧急集合

109 篇文章 4 订阅
10 篇文章 0 订阅

Description

Y岛风景美丽宜人,气候温和,物产丰富。Y岛上有N个城市,有N-1条城市间的道路连接着它们。每一条道路都连接某两个城市。幸运的是,小可可通过这些道路可以走遍Y岛的所有城市。神奇的是,乘车经过每条道路所需要的费用都是一样的。小可可,小卡卡和小YY经常想聚会,每次聚会,他们都会选择一个城市,使得3个人到达这个城市的总费用最小。
由于他们计划中还会有很多次聚会,每次都选择一个地点是很烦人的事情,所以他们决定把这件事情交给你来完成。他们会提供给你地图以及若干次聚会前他们所处的位置,希望你为他们的每一次聚会选择一个合适的地点。

Input

第一行两个正整数,N和M。分别表示城市个数和聚会次数。后面有N-1行,每行用两个正整数A和B表示编号为A和编号为B的城市之间有一条路。城市的编号是从1到N的。再后面有M行,每行用三个正整数表示一次聚会的情况:小可可所在的城市编号,小卡卡所在的城市编号以及小YY所在的城市编号。

Output

一共有M行,每行两个数Pos和Cost,用一个空格隔开。表示第i次聚会的地点选择在编号为Pos的城市,总共的费用是经过Cost条道路所花费的费用。

Sample Input

6 4
1 2
2 3
2 4
4 5
5 6
4 5 6
6 3 1
2 4 4
6 6 6

Sample Output

5 2
2 5
4 1
6 0

HINT

数据范围:
100%的数据中,N<=500000,M<=500000。
40%的数据中N<=2000,M<=2000。

题解

upd 2017/12/15
双倍经验 1832&1787一样的,样例都一样输出都没改= =

哇塞一道水题,bzoj居然会有这样的好东西嘿嘿嘿
这个。。很明显可以看出,答案一定在这三个数的lca之一嘛。口胡一下,如果聚会点不在lca上,譬如在lca的父亲节点。那么深度较深的两个节点走到这里的花费就要多2,而较浅的节点走到这里的花费只少了1。暴力三个lca然后算出余下的一个节点到这个lca的权,记录答案即可
因为这题所有边权相同哈哈哈哈哈就可以lca水过了

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    int x,y,next;
}a[1110000];int len,last[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[22],n,m;
int fa[510000][22],dep[510000];
void pre_tree_node(int x)
{
    for(int i=1;i<=20 && bin[i]<=dep[x];i++)fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=fa[x][0])
        {
            fa[y][0]=x;
            dep[y]=dep[x]+1;
            pre_tree_node(y);
        }
    }
}
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=fa[x][i];
    if(x==y)return x;
    for(int i=20;i>=0;i--)if(dep[x]>=bin[i] && fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}
int u[4];
bool cmp(int n1,int n2){return dep[n1]<dep[n2];}
int main()
{
    scanf("%d%d",&n,&m);
    bin[0]=1;
    for(int i=1;i<=20;i++)bin[i]=bin[i-1]*2;
    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);
    }
    dep[1]=1;fa[1][0]=0;pre_tree_node(1);
    while(m--)
    {
        int ans=999999999;
        scanf("%d%d%d",&u[1],&u[2],&u[3]);
        sort(u+1,u+1+3,cmp);
        int num,cnt,tmp,ret;
        cnt=0;;
        num=lca(u[1],u[2]);
        cnt+=dep[u[1]]+dep[u[2]]-dep[num]*2;
        tmp=lca(num,u[3]);
        cnt+=dep[u[3]]+dep[num]-dep[tmp]*2;
        if(ans>cnt)ans=cnt,ret=num;

        cnt=0;
        num=lca(u[1],u[3]);
        cnt+=dep[u[1]]+dep[u[3]]-dep[num]*2;
        tmp=lca(num,u[2]);
        cnt+=dep[u[2]]+dep[num]-dep[tmp]*2;
        if(ans>cnt)ans=cnt,ret=num;

        cnt=0;
        num=lca(u[2],u[3]);
        cnt+=dep[u[2]]+dep[u[3]]-dep[num]*2;
        tmp=lca(num,u[1]);
        cnt+=dep[u[1]]+dep[num]-dep[tmp]*2;
        if(ans>cnt)ans=cnt,ret=num;
        printf("%d %d\n",ret,ans);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值