BZOJ3242 快餐店【基环树直径】

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接: https://blog.csdn.net/C20181220_xiang_m_y/article/details/102596355

题目描述:

一棵正权边基环树,在树上选一个点(可以在边上),使得所有节点到它的距离的最大值最小,输出最小值。
3<=n<=105

题目描述:

如果是一棵树,我们比较容易证明直径的中点就是最优点。

考虑在基环树中有什么变化,多出了一条边。如果我们断掉一条边求直径,得出一个距离 d 2 \frac d2 ,显然这个距离是大于等于真实的所有点到这个直径中点的最大值的。

考虑答案的情形,假设最优点是 x x ,将所有点到 x x 的最短路径经过的边标记,那么基环树的环上至少有一条边不会被标记,所以删去这条边不会影响答案,当我们删去这条边时,最优点就是树的直径中点,所以答案>=直径二分之一,直径二分之一>=真实距离,而真实距离>=答案,所以答案就是此时直径的二分之一。

总的来说,我们枚举断掉一条边,求出直径,直径的二分之一一定大于等于答案,而答案一定会在枚举某条边时被计算到。

计算直径可以以1号点左右划分,用前后缀预处理然后端点相加,与每个子树中的直径取个max。摆一张图自行体会一下:
在这里插入图片描述

Code:

#include<bits/stdc++.h>
#define maxn 100005
#define LL long long
using namespace std;
int n,fa[maxn],dep[maxn],cir[maxn],D[maxn],m;
LL ans,f[maxn],g1[maxn],g2[maxn],h1[maxn],h2[maxn];
//f:maxdep  g:double chain  h:single chain
bool onc[maxn];
int fir[maxn],nxt[maxn<<1],to[maxn<<1],w[maxn<<1],tot;
inline void bothline(int x,int y,int z){
	nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,w[tot]=z;
	nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,w[tot]=z;
}
void findcircle(int u,int ff){
	dep[u]=dep[fa[u]=ff]+1;
	for(int i=fir[u],v;i;i=nxt[i]) if((v=to[i])!=ff){
		if(!dep[v]) findcircle(v,u);
		else if(dep[v]<dep[u])
			for(int j=u;;j=fa[j]) {onc[cir[++m]=j]=1;if(j==v) break;}
	}
}
void dfs(int u,int ff){
	for(int i=fir[u],v;i;i=nxt[i]) if((v=to[i])!=ff&&!onc[v]){
		dfs(v,u),f[v]+=w[i];
		ans=max(ans,f[u]+f[v]),f[u]=max(f[u],f[v]);
	}
}
int main()
{
	scanf("%d",&n);
	for(int i=1,x,y,z;i<=n;i++) scanf("%d%d%d",&x,&y,&z),bothline(x,y,z);
	findcircle(1,0);
	for(int i=1;i<=m;i++) dfs(cir[i],0);
	cir[m+1]=cir[1];
	for(int k=1;k<=m;k++)
		for(int i=fir[cir[k]];i;i=nxt[i]) if(to[i]==cir[k+1]) {D[k]=w[i];break;}
	LL sum=0,mx=f[cir[1]]+D[1];
	for(int i=2;i<=m;i++){
		sum+=D[i-1];
		h1[i]=max(h1[i-1],f[cir[i]]+sum);
		g1[i]=max(g1[i-1],f[cir[i]]+mx);
		mx=max(mx,f[cir[i]])+D[i];
	}
	sum=0,mx=f[cir[1]]+D[m];
	for(int i=m;i>1;i--){
		sum+=D[i];
		h2[i]=max(h2[i+1],f[cir[i]]+sum);
		g2[i]=max(g2[i+1],f[cir[i]]+mx);
		mx=max(mx,f[cir[i]])+D[i-1];
	}
	LL Ans=1ll<<60;
	for(int i=1;i<=m;i++) Ans=min(Ans,max(ans,max(max(g1[i],g2[i+1]),h1[i]+h2[i+1])));
	printf("%.1f\n",Ans/2.0);
}

在最短路问题中,这种计算 x x ,使得 A n s { x x A n s } Ans\in\{x|x\ge Ans\} 的思想很重要,应用也比较多。

猜你喜欢

转载自blog.csdn.net/C20181220_xiang_m_y/article/details/102596355