图的半径

题目大意

给定无向图，求图的半径。
图的半径定义：找到图上一点（可以在边上），使得其余所有点到该点最短路最大的那个最小，这个最小的最大最短路就是图的半径。该点称为中心。

观察

先做floyd跑双源最短路
我们枚举一条边j连向k，长度为t，认为中心在这条边上。
那么答案$min(max_{i=1}^n(min(dis[j][i]+x,dis[k][i]+t-x))）)$
x表示中心和j的距离，取值范围[0,t]
我们观察每个点到中心最短路关于x的函数，显然是个严格单峰函数或严格单调函数。
如果它是单峰函数，意味着什么？
我们定义A集合表示所有走到中心的最短路是经过j的点，B集合则是选择走k。
x=0时每个点都属于某集合，那么单峰函数就是说x到达某个值后会切换集合。
切换集合的次数只有n次，我们可以求出x是什么时出现某点切换集合，称会切换集合的为关键点。
顺次枚举n个关键点，两个关键点间每个点所属集合关系是确定的。对于A集合维护dis[j][i]的最大值Amx，B集合也维护个Bmx，那么要最大化$max(Amx+x,Bmx+t-x)$
相等时最优化，可以解出x的最优取值。当然别忘记如果x最优情况不在目前两个关键点区间内，那么最优值就是在两个关键点上。
Amx和Bmx的维护可以使用数据结构。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef double db;
const int maxn=200+10,maxm=50000+10;
const db eps=0.001;
int dis[maxn][maxn],from[maxm],go[maxm],dist[maxm],id[maxn],wz[maxn];
bool bz[maxn];
struct dong{
    db x;
    int id;
} a[maxn],b[maxn];
int i,j,k,t,n,m,top,tot,u,v,c;
db ans,p,q,x,l,r,c1,c2;
bool cmp(dong a,dong b){
    return a.x<b.x;
}
int main(){
    freopen("radius.in","r",stdin);freopen("radius.out","w",stdout);
    scanf("%d%d",&n,&m);
    fo(i,1,n)
        fo(j,1,n)
            if (i!=j) dis[i][j]=1000000000;
    fo(i,1,m){
        scanf("%d%d%d",&from[i],&go[i],&dist[i]);
        dis[from[i]][go[i]]=min(dis[from[i]][go[i]],dist[i]);
        dis[go[i]][from[i]]=min(dis[go[i]][from[i]],dist[i]);
    }
    fo(k,1,n)
        fo(i,1,n)
            if (i!=k)
                fo(j,1,n)
                    if (j!=k)
                        if (i!=j)
                            dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    ans=1000000000;
    fo(i,1,m){
        q=0;
        u=from[i];v=go[i];c=dist[i];
        a[1].x=a[1].id=a[2].id=0;a[top=2].x=c;
        tot=0;
        fo(j,1,n){
            if (dis[j][u]<dis[j][v]+c){
                p=(db)(dis[j][v]+c-dis[j][u])/2;
                if (0<p&&p<c){
                    a[++top].x=p;
                    a[top].id=j;
                }
                b[++tot].x=dis[j][u];
                b[tot].id=j;
            }
            else if (q<dis[j][v]) q=dis[j][v];
        }
        sort(a+1,a+top+1,cmp);
        sort(b+1,b+tot+1,cmp);
        fo(j,1,tot) bz[wz[b[j].id]=j]=0;
        fo(j,1,top-1){
            if (a[j].id){
                bz[wz[a[j].id]]=1;
                if (q<dis[a[j].id][v]) q=dis[a[j].id][v];
            }
            while (tot&&bz[tot]) tot--;
            c1=b[tot].x;
            c2=q;
            x=(c2+c-c1)/2;
            if (a[j].x<=x&&x<=a[j+1].x) ans=min(ans,c1+x);
            else ans=min(ans,min(max(c1,c2+c),max(c1+c,c2)));
        }
    }
    printf("%.2lf\n",ans);
}