[51Nod2558] 选址

link

考虑二分答案 $F$ ，那么现在的问题变成是否对于覆盖并有交集。

考虑边 $(u,v)$ ，若覆盖并在 $(u,v,w)$ 线段中，设点 $i$ 走到 $u$ 号后还能走 $F1$ , 走到 $v$ 还能走 $F2$ ，则现在要求的是一个子问题:求在 $n$ 个 $(0,F1),(w-F2+1,w)$ 中判断是否有交集，若存在点 $x$ 使得 $n$ 个线段都能被覆盖时，$x$ 肯定为 $F1$ 或 $w-F2+1$ ，所以直接暴力枚举即可，时间复杂度 $O(n^4\log 值域)$。

而如何优化时间复杂度，考虑如何快速计算贡献，直接离散化后差分计算即可。

时间复杂度 $O(n^3\log 值域)$。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
using namespace std;
inline int read(){
    int f=,ans=;char c=getchar();
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){ans=ans*+c-'';c=getchar();}
    return f*ans;
}
const int MAXN=;
struct Edge{
    int u,v,w;
}x[MAXN*MAXN];
int n,m,dis[MAXN][MAXN];
double f1[MAXN][MAXN*MAXN],f2[MAXN][MAXN*MAXN],l,r,minn=LLONG_MAX,eps=1e-,num[MAXN<<],tmp[MAXN<<];
int M;
int G1[MAXN<<][MAXN<<],G2[MAXN<<][MAXN<<];
int NU[MAXN<<],cnt[MAXN],Cnt;
inline int Q(double res){return lower_bound(tmp+,tmp+M+,res)-tmp;}
inline bool Query(int id,double G){
    memset(NU,,sizeof(NU));memset(cnt,,sizeof(cnt));
    int Num=;
    num[++Num]=;tmp[Num]=num[Num];
    num[++Num]=x[id].w;tmp[Num]=num[Num];
    for(register int i=;i<=n;i++){
        double t1=f1[i][id],t2=x[id].w-f2[i][id];
        Cnt++;
        if(t1>=) num[++Num]=min(t1,(double)x[id].w),tmp[Num]=num[Num];
        if(f2[i][id]>=) num[++Num]=max((double),t2),tmp[Num]=num[Num];
    }
    sort(tmp+,tmp+Num+);
    M=unique(tmp+,tmp+Num+)-tmp-;
    for(register int i=;i<=Num;i++) Cnt++,num[i]=lower_bound(tmp+,tmp+M+,num[i])-tmp;
    for(register int i=;i<=n;i++){
        Cnt++;
        double t1=f1[i][id],t2=x[id].w-f2[i][id];
        if(t1>=){
            int l=,r=Q(min(t1,(double)x[id].w))+;
            G1[l][++NU[l]]=i;G2[l][NU[l]]=;
             G1[r][++NU[r]]=i;G2[r][NU[r]]=-;
        }
        if(f2[i][id]>=){
            int l=Q(max((double),t2));
            G1[l][++NU[l]]=i;G2[l][NU[l]]=;
        }
    }
    int Sum=;
    for(register int i=;i<=M;i++){
        for(register int j=;j<=NU[i];j++){
            Cnt++;
            int f=G1[i][j],opt=G2[i][j];
            if(opt==){
                if(cnt[f]==) Sum++;
                cnt[f]++;
            }else{
                if(cnt[f]==) Sum--;
                cnt[f]--;
            }
        }if(Sum==n) return ;
    }return ;
}
inline bool check(double G){
    for(register int u=;u<=n;u++){
        for(register int i=;i<=m;i++){
            f1[u][i]=G-dis[u][x[i].u];
            f2[u][i]=G-dis[u][x[i].v];
        }
    }
    for(register int i=;i<=m;i++)
        if(Query(i,G)) return ;
    return ;
}
int main(){
//    freopen("make.in","r",stdin);
    n=read(),m=read();
    memset(dis,/,sizeof(dis));
    for(int i=;i<=m;i++){
        int u=read(),v=read(),w=read();
        r+=w;
        x[i].u=u,x[i].v=v,x[i].w=w;
        dis[u][v]=dis[v][u]=w;
    }
    for(int i=;i<=n;i++) dis[i][i]=;
    for(int p=;p<=n;p++)
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                dis[i][j]=min(dis[i][j],dis[i][p]+dis[p][j]);
    while(l<=r){
        double mid=(l+r)/2.0;
        if(check(mid)) minn=min(minn,mid),r=mid-eps;
        else l=mid+eps;
    }printf("%.9lf\n",minn);
}