[JZOJ4884]图的半径/[JZOJ4572]Center

题目大意

给定一个$n$个点，$m$条边的无向带权图。在图上确定一个点，使得该点到图中最远点距离最近。点可以在路径上。答案保留两位小数。

$1\le n\le200,1\le m\le19900$，边权$w\le10^5$

---

题目分析

我们枚举边，计算点在每条边上时的最优答案。
令当前枚举的边是$(a,b)$，边长为$l$，我们设选取的点距离$a$的距离为$x$。则其答案为

$$\max_{i=1}^n\{\min({dis}_{a,i}+x,{dis}_{b,i}+l-x)\}$$
可以发现我们能够计算出对于一个点，$x$在什么范围的时候$dis_{a,i}<dis_{b,i}+l-x$，显然这个临界点是$x={dis_{b,i}-dis_{a,i}+l\over2}$。那么我们将所有点的临界$x$值计算出来并排序，然后按顺序依次处理，就能够将所有点分成距离$a$更近的点和距离$b$更近的点两个集合$A$与$B$。显然答案为
$$\max\left(\max\{A_i\}+x,\max\{B_i\}+l-x\right)$$
$x$能使两个元素取等号时答案最小，那么我们可以把这个$x$计算出来，同时还要检验一下这个$x$是否在令当前集合分组合法的范围内（即相邻的两个临界点内），如果不在，那么答案一定在两个临界点中的一个上。
我们需要一个能够支持加入、删除查询最大值的数据结构，使用set即可。当然，由于这个过程是一个集合不断删数，一个集合不断加数，正反两次预处理即可得到最大值。
最短路使用$Floyd$算法就好了。
时间复杂度$\mathrm O(nm\log_2n+n^3)$。

---

代码实现

#include <algorithm>
#include <iostream>
#include <climits>
#include <cstdio>
#include <cctype>
#include <set>

using namespace std;

typedef long double db;

inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) f=ch=='-'?-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

db max(db x,db y){return x>y?x:y;}
db min(db x,db y){return x<y?x:y;}

const int INF=INT_MAX/2;
const int N=205;
const int M=19950;
const int E=M<<1;

struct edge
{
    int x,y,l;
}e[M];

struct key
{
    db x;
    int id;
}pts[N];

bool operator<(key a,key b){return a.x<b.x;}

int dis[N][N];
int n,m,tot;
db ans;

void floyd()
{
    for (int k=1;k<=n;k++)
        for (int i=1;i<=n;i++)
            if (i!=k)
                for (int j=1;j<=n;j++)
                    if (i!=j&&j!=k)
                        dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}

db mx[2][N];

db getans(int x,db l,db r,db w)
{
    db a=mx[1][x],b=mx[0][x-1];
    db ret=(b-a+w)/2.0;
    if (l<=ret&&ret<=r) return a+ret;
    return min(max(a+l,b+w-l),max(a+r,b+w-r));
}

void calc()
{
    ans=INF;
    for (int i=1,x,y,l,cur;i<=m;i++)
    {
        x=e[i].x,y=e[i].y,l=e[i].l,tot=0;
        mx[0][0]=0.0,mx[1][tot+1]=0.0;
        for (int j=1;j<=n;j++)
        {
            db tmp=(dis[y][j]-dis[x][j]+l)/2.0;
            if (0<=tmp&&tmp<=l) pts[++tot].x=tmp,pts[tot].id=j;
            else if (tmp<0) mx[1][tot+1]=max(mx[1][tot+1],dis[y][j]);
            else mx[0][0]=min(mx[0][0],dis[x][j]);
        }
        if (!tot) continue;
        sort(pts+1,pts+1+tot);
        pts[0].x=0.0;
        for (int j=1;j<=tot;j++) mx[0][j]=max(mx[0][j-1],dis[y][pts[j].id]);
        for (int j=tot;j>=1;j--) mx[1][j]=max(mx[1][j+1],dis[x][pts[j].id]);
        for (int j=1;j<=tot;j++) if (pts[j].x!=pts[j-1].x) ans=min(ans,getans(j,pts[j-1].x,pts[j].x,l));
        ans=min(ans,getans(tot+1,pts[tot].x,l,l));
    }
}

int main()
{
    freopen("radius.in","r",stdin),freopen("radius.out","w",stdout);
    n=read(),m=read();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            dis[i][j]=INF;
    for (int i=1,x,y,z;i<=m;i++)
    {
        x=read(),y=read(),z=read();
        dis[x][y]=min(dis[x][y],z),dis[y][x]=min(dis[y][x],z);
        e[i].x=x,e[i].y=y,e[i].l=z;
    }
    for (int i=1;i<=n;i++) dis[i][i]=0;
    floyd(),calc();
    printf("%.2lf\n",(double)ans);
    fclose(stdin),fclose(stdout);
    return 0;
}