Wannafly挑战赛1 C MMSet2(lca)

本文介绍了一种基于树形结构的动态规划算法,用于解决特定类型的最短路径问题。通过预处理得到每个节点的深度及最近公共祖先(LCA),算法能够在给定节点集合中找到两点间距离的最大值对应的中点,进而计算出该集合的最小最大值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

MMSet2
时间限制:3秒 空间限制:131072K

题目描述

给定一棵n个节点的树,点编号为1…n。
Q次询问,每次询问给定一个点集S,令
你需要求出
其中dist(u,v)表示树上路径(u,v)的边数。

输入描述:

第一行一个整数n,接下来n−1行每行两个整数表示树上的一条边。
接下来一行一个整数Q,接着Q行,每行第一个数是|S|,剩下|S|个互不相同的数代表这个集合。

输出描述:

输出Q行,每行一个整数表示答案。
示例1

输入

3
1 2
1 3
1
2 2 3

输出

1

备注:

 n≤3×105,|S|≥1,∑|S|≤106




这题要求最大值得最小值,那么答案点一定就在集合中两个点的中心点

因为在这条路径外的点到这两个点的路径长度虽然大,但是求得是最小值所以 一定是最长的一条路径的中心点

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<vector>
#include<map>
#include <bits/stdc++.h>
using namespace std;
const int N = 300000+10;
typedef  long long LL;
const LL mod = 998244353;
int a[N];
vector<int>p[N];
int dep[N],fa[N][20];
void dfs(int u,int d,int f)
{
    dep[u]=d,fa[u][0]=f;
    for(int i=0; i<p[u].size(); i++)
    {
        int v=p[u][i];
        if(v==f) continue;
        dfs(v,d+1,u);
    }
    return ;
}
void init(int n)
{
    int k=(int)(1.0*log(1.0*n)/log(2.0));
    for(int i=1;i<=k;i++)
        for(int j=1;j<=n;j++)
            fa[j][i]=fa[fa[j][i-1]][i-1];
    return ;
}
int get(int x,int y,int n)
{
    if(dep[x]<dep[y]) swap(x,y);
    int d=dep[x]-dep[y];
    int k=(int)(log(1.0*n)/log(1.0*2));
    for(int i=0;i<=k;i++)
        if((1<<i)&d)x=fa[x][i];
    if(x==y) return x;
    if(x!=y)
    {
        for(int i=k;i>=0;i--)
            if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
    }
    return fa[x][0];
}
int lca(int x,int y,int n)
{
    int z=get(x,y,n);
    return dep[x]+dep[y]-2*dep[z];
}

int main()
{
    int n;
    scanf("%d", &n);
    for(int i=1; i<n; i++)
    {
        int x, y;
        scanf("%d %d", &x, &y);
        p[x].push_back(y),p[y].push_back(x);
    }
    dep[0]=0,fa[0][0]=0;
    dep[1]=1;
    dfs(1,1,0);
    init(n);
    int q;
    scanf("%d", &q);
    while(q--)
    {
        int m, ans=0;
        scanf("%d", &m);
        int root=0;
        for(int i=1; i<=m; i++)
        {
            scanf("%d", &a[i]);
            if(dep[a[i]]>dep[root]) root=a[i];
        }
        for(int i=1;i<=m;i++)
            if(a[i]!=root) ans=max(ans,lca(root,a[i],n));
        printf("%d\n",(ans+1)/2);
    }
    return 0;
}




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值