小T打算在城市C开设一家外送快餐店。送餐到某一个地点的时间与外卖店到该地点之间最短路径长度是成正比的,小T希望快餐店的地址选在离最远的顾客距离最近的地方。 快餐店的顾客分布在城市C的N 个建筑中,这N 个建筑通过恰好N 条双向道路连接起来,不存在任何两条道路连接了相同的两个建筑。任意两个建筑之间至少存在一条由双向道路连接而成的路径。小T的快餐店可以开设在任一建筑中,也可以开设在任意一条道路的某个位置上(该位置与道路两端的建筑的距离不一定是整数)。 现给定城市C的地图(道路分布及其长度),请找出最佳的快餐店选址,输出其与最远的顾客之间的距离。
第一行包含一个整数N,表示城市C中的建筑和道路数目。
接下来N行,每行3个整数,Ai,Bi,Li(1≤i≤N;Li>0),表示一条道路连接了建筑Ai与Bi,其长度为Li 。
仅包含一个实数,四舍五入保留恰好一位小数,表示最佳快餐店选址距离最远用户的距离。
注意:你的结果必须恰好有一位小数,小数位数不正确不得分。
数据范围
对于 10%的数据,N<=80,Li=1;
对于 30%的数据,N<=600,Li<=100;
对于 60% 的数据,N<=2000,Li<=10^9;
对于 100% 的数据,N<=10^5,Li<=10^9
//MADE BY QT666
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#define int long long
#define lson x<<1
#define rson x<<1|1
using namespace std;
typedef long long ll;
const int N=300050;
int fa[N],dfn[N],tt,head[N],to[N],nxt[N],w[N],cnt;
int len,cir[N],v[N],dis[N],bj[N],n,d[N],p[N],f[N],ans;
int pre[N],suf[N],tc1[N],tc2[N];
void lnk(int x,int y,int z){
to[++cnt]=y,nxt[cnt]=head[x],w[cnt]=z,head[x]=cnt;
to[++cnt]=x,nxt[cnt]=head[y],w[cnt]=z,head[y]=cnt;
}
void getrt(int rt,int x){
for(int u=x;u!=fa[rt];u=fa[u])
cir[++len]=u,v[len+1]=dis[u],bj[u]=1;
}
void Tarjan(int x,int ff){
dfn[x]=++tt;fa[x]=ff;
for(int i=head[x];i;i=nxt[i]){
int y=to[i];
if(!dfn[y]) dis[y]=w[i],Tarjan(y,x);
else if(dfn[y]>dfn[x]) v[1]=w[i],getrt(x,y);
}
}
void dfs(int x,int ff,int top){
for(int i=head[x];i;i=nxt[i]){
int y=to[i];if(y==ff||bj[y]) continue;
dfs(y,x,top);
p[top]=max(p[top],f[x]+f[y]+w[i]);f[x]=max(f[x],f[y]+w[i]);
}
}
main(){
freopen("foodshop.in","r",stdin);
freopen("foodshop.out","w",stdout);
scanf("%lld",&n);
for(int i=1;i<=n;i++){int x,y,z;scanf("%lld%lld%lld",&x,&y,&z);lnk(x,y,z);}
Tarjan(1,0);
for(int i=1;i<=len;i++) dfs(cir[i],cir[i],cir[i]),d[cir[i]]=f[cir[i]];
int ans=1ll<<60,sum,pr;sum=0,pr=0;
for(int i=2;i<=len;i++){
pre[i]=max(pre[i-1],sum+d[cir[i-1]]);
tc1[i]=max(tc1[i-1],sum+d[cir[i-1]]+pr);
pr=max(pr,d[cir[i-1]]-sum);
sum+=v[i];
}
pre[1]=max(pre[len],sum+d[cir[len]]);
tc1[1]=max(tc1[len],sum+d[cir[len]]+pr);
sum=0;sum+=v[1];suf[len]=sum+d[cir[len]];
pr=d[cir[len]]-sum;
for(int i=len-1;i>=2;i--){
suf[i]=max(suf[i+1],sum+d[cir[i]]+v[i+1]);
tc2[i]=max(tc2[i+1],sum+d[cir[i]]+v[i+1]+pr);
pr=max(pr,d[cir[i]]-(sum+v[i+1]));
sum+=v[i+1];
}
for(int i=1;i<=len;i++){
ans=min(ans,max(pre[i]+suf[i],max(tc1[i],tc2[i])));
}
int tmp=0;
for(int i=1;i<=len;i++){
tmp=max(tmp,p[cir[i]]);
ans=max(ans,p[cir[i]]);
}
printf("%.1f\n",1.0*ans/2);
return 0;
}