首页 > 其他 > 详细

BZOJ1832 [AHOI2008]聚会

时间:2019-06-13 13:26:13      阅读:109      评论:0      收藏:0      [点我收藏+]

1832: [AHOI2008]聚会

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 2132  Solved: 911
[Submit][Status][Discuss]

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

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

HINT

Source

[Submit][Status][Discuss]
?
HOME Back

题解

画个图就很清楚要找的是什么点了。以这个要找的节点为根后,三个询问节点都在不同的子树内(或者与他重合)。

两两求lca后,有两个相同且深度较浅的,有一个深度较深的,要找的就是后者。

那么枚举判断一下就好了。时间复杂度\(O((n+m)\log n)\)

#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<class T>il T read(){
    rg T data=0,w=1;rg char ch=getchar();
    for(;!isdigit(ch);ch=getchar())if(ch=='-') w=-w;
    for(;isdigit(ch);ch=getchar()) data=data*10+ch-'0';
    return data*w;
}
template<class T>il T read(rg T&x) {return x=read<T>();}
typedef long long ll;
using namespace std;

co int N=5e5+1;
int n,m,t,f[N][20],d[N];
bool v[N];
queue<int> q;
vector<int> e[N];

void bfs(){
    q.push(1);
    v[1]=1,d[1]=1;
    while(q.size()){
        int x=q.front();q.pop();
        for(unsigned i=0;i<e[x].size();++i){
            int y=e[x][i];
            if(v[y]) continue;
            q.push(y);
            v[y]=1,d[y]=d[x]+1;
            f[y][0]=x;
            for(int j=1;j<=t;++j)
                f[y][j]=f[f[y][j-1]][j-1];
        }
    }
}
int lca(int x,int y){
    if(d[x]>d[y]) swap(x,y);
    for(int i=t;i>=0;--i)
        if(d[f[y][i]]>=d[x]) y=f[y][i];
    if(x==y) return x;
    for(int i=t;i>=0;--i)
        if(f[x][i]!=f[y][i])  x=f[x][i],y=f[y][i];
    return f[x][0];
}
int main(){
    t=log(read(n))/log(2),read(m);
    for(int i=1,x,y;i<n;++i){
        read(x),read(y);
        e[x].push_back(y),e[y].push_back(x);
    }
    bfs();
    while(m--){
        int x=read<int>(),y=read<int>(),z=read<int>();
        int xy=lca(x,y),yz=lca(y,z),xz=lca(x,z);
        int s=d[x]+d[y]-2*d[xy];
        if(d[xy]<d[yz]){
            xy=yz;
            s=d[y]+d[z]-2*d[yz];
            z=x;
        }
        if(d[xy]<d[xz]){
            xy=xz;
            s=d[x]+d[z]-2*d[xz];
            z=y;
        }
        s+=d[z]+d[xy]-2*d[lca(z,xy)];
        printf("%d %d\n",xy,s);
    }
    return 0;
}

BZOJ1832 [AHOI2008]聚会

原文:https://www.cnblogs.com/autoint/p/11015803.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!