bzoj3124: [Sdoi2013]直径

该博客介绍了如何解决[Sdoi2013]直径问题。通过两个关键结论,确定直径必须在同一条线上且连续,利用这些信息遍历并确定直径上的路径。通过正向和反向询问,排除不可能构成直径的边,最终确定直径的范围并实现求解算法。

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

题目传送门

解法:
两个傻逼结论:
首先肯定在同一条直径上。因为每条直径都要经过这一段。
肯定是连续的一段。
因为如果你有两段的话。那么敢问这两段之间的路径你是怎么走过去的(树上两点路径唯一!)

这样的话就很好做了。
随便找条直径。
在这条直径上确定一条路径。使得这条路径上:
任意一点都不可能从别的分支扩展出直径。
懒得画图直接文字描述。
假设求出了一条直径是1->2->5->7
那么从2开始问(因为1肯定没有儿子了,要不然它不会是直径的端点)
那么我们可以记录1到2的距离,假设为t。
然后看看2能不能从别的分支走出一个距离t。如果可以那么2之前的边都不在所有直径的交集上

问到5就记录下1到5的距离,假设为T。
然后看看5能不能从别的分支走出一个T。可以的话那么5之前的边也都不行。
以此类推。

正着问一次倒着问一次。
就确定出了范围。
然后就可以求解了呀。

代码实现:

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
struct node {int x,y,next;ll d;}a[510000];int len,last[210000];
void ins(int x,int y,ll d) {len++;a[len].x=x;a[len].y=y;a[len].d=d;a[len].next=last[x];last[x]=len;}
int G,fa[210000],bian[210000];ll D;
void dfs(int x,ll d,int f) {
    if(D<d) {D=d;G=x;}
    for(int k=last[x];k;k=a[k].next) {
        int y=a[k].y;
        if(y!=f){fa[y]=x;bian[y]=k;dfs(y,d+a[k].d,x);}
    }
}
int s[210000];
ll T;bool v[210000];
void dfs1(int x,ll d,int f) {
    if(T<d)T=d;
    for(int k=last[x];k;k=a[k].next) {
        int y=a[k].y;
        if(v[y]==true&&y!=f) dfs1(y,d+a[k].d,x);
    }
}
int ans;
void dfs2(int x,int d,int f,int ed) {
    if(x==ed){ans=d;return ;}
    for(int k=last[x];k;k=a[k].next) {
        int y=a[k].y;
        if(y!=f) dfs2(y,d+1,x,ed);
    }
}
int main() {
    int n;scanf("%d",&n);len=0;memset(last,0,sizeof(last));
    for(int i=1;i<n;i++) {
        int x,y;ll d;scanf("%d%d%lld",&x,&y,&d);
        ins(x,y,d);ins(y,x,d);
    }
    int X,Y;
    memset(bian,0,sizeof(bian));dfs(1,0,0);X=G;D=0;dfs(X,0,0);Y=G;
    int x=Y,slen=0;memset(v,true,sizeof(v));
    while(x!=X) {s[++slen]=x;v[x]=false;x=fa[x];}
    s[++slen]=X;v[X]=false;
    int ansx=Y;ll K=a[bian[s[1]]].d;
    for(int i=2;i<=slen;i++) {
        T=0;dfs1(s[i],0,0);
        if(T==K)ansx=s[i];K+=a[bian[s[i]]].d;
    }
    for(int i=slen;i>=2;i--) bian[s[i]]=bian[s[i-1]];bian[s[1]]=0;
    int ansy=X;K=a[bian[s[slen]]].d;
    for(int i=slen-1;i>=1;i--) {
        T=0;dfs1(s[i],0,0);
        if(T==K)ansy=s[i];K+=a[bian[s[i]]].d;
    }dfs2(ansx,0,0,ansy);printf("%lld\n%d\n",D,ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值