JZOJ4884. 【NOIP2016提高A组集训第12场11.10】图的半径

Description

mhy12345学习了树的直径，于是开始研究图的半径，具体来说，我们需要在图中选定一个地方作为中心，其中这个中心有可能在路径上。
而这个中心的选址需要能够使得所有节点达到这个中心的最短路里面最大值最小(也就是说离中心最远的结点距离尽可能小)，求出这个最大值的最小值，作为图的半径。

Data Constraint

N<=200，W<=100000，M<=19900
数据有梯度

Solution

对于一条边(x,y)，我们可以先假设所有点都到达x点，每次将这些点中到达x最大的点转移到y点，并更新y点上的最大值，这样我们就可一直到到两个点(x,y)的最大值分别是什么，进而求得中间点，时间复杂度O($NMlogN$)。

code

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
const int maxn=50005,maxn1=205;
struct code{
    int a,b;
}g1[maxn1];
int v[maxn],d[maxn],first[maxn],last[maxn],next[maxn],value[maxn],f[maxn1][maxn1];
int n,i,t,j,m,k,l,x,y,z,g[maxn1][maxn1],a,num,b,p,q,ans,xx;
bool bz[maxn],bz1;
double ans1;
void lian(int x,int y,int z){
    last[++num]=y;next[num]=first[x];first[x]=num;value[num]=z;
}
bool cmp(code x,code y){
    return x.a>y.a || x.a==y.a && x.b<y.b;
}
int main(){
    freopen("radius.in","r",stdin);freopen("radius.out","w",stdout);
    //freopen("data.in","r",stdin);
    scanf("%d%d",&n,&m);
    memset(g,127,sizeof(g));l=g[1][1];
    for (i=1;i<=n;i++)
        for (j=1;j<=n;j++)
            if (j!=i) f[i][j]=1e9;
    for (i=1;i<=m;i++)
        scanf("%d%d%d",&x,&y,&z),f[x][y]=f[y][x]=g[x][y]=g[y][x]=min(z,f[x][y]);
    ans=1e9;
    for (k=1;k<=n;k++)
        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
                f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    for (i=1;i<=n;i++){
        t=0;
        for (j=1;j<=n;j++)
                t=max(t,f[i][j]);
        ans=min(ans,t*2);
    }
    t=0;
    for (x=1;x<=n;x++)  
        for (y=x+1;y<=n;y++)
            if (g[x][y]!=l){
                p=0;
                for (i=1;i<=n;i++) g1[i].a=f[i][y],g1[i].b=f[i][x];
                sort(g1+1,g1+n+1,cmp);
                t=0;
                for (i=1;i<n;i++){
                    if (g1[i].b>p) p=g1[i].b;a=p;
                    b=g1[i+1].a;
                    if (a>b) swap(a,b);
                    if (a+g[x][y]<=b){
                        if (b*2<ans) ans=b*2;
                    }
                    else if (b+a+g[x][y]<ans) ans=b+a+g[x][y];
                }
            }
    ans1=ans*1.0/2;
    printf("%.2lf\n",ans1);
}