[Wc2010]重建计划

Description

Input

第一行包含一个正整数N，表示X国的城市个数. 第二行包含两个正整数L和U，表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案，每行三个正整数Ai,Bi,Vi分别表示道路(Ai,Bi),其价值为Vi 其中城市由1..N进行标号

Output

输出最大平均估值，保留三位小数

Sample Input

4
2 3
1 2 1
1 3 2
1 4 3

Sample Output

2.500

HINT

N<=100000,1<=L<=U<=N-1,Vi<=1000000

题解在这里

ZYYS

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 struct Node
 {
   int next,to;
   double d,c;
 }edge[];
 int num,head[],size[],maxsize[],minsize,root;
 int rt[],tot,n,fa[],L,U,dep[];
 double mx[],dist[];
 bool vis[];
 void add(int u,int v,double d)
 {
   num++;
   edge[num].next=head[u];
   head[u]=num;
   edge[num].to=v;
   edge[num].d=d;
 }
 void get_size(int x,int pa)
 {int i;
   size[x]=;maxsize[x]=;
   for (i=head[x];i;i=edge[i].next)
     {
       int v=edge[i].to;
       if (v==pa||vis[v]) continue;
       get_size(v,x);
       size[x]+=size[v];
       if (size[v]>maxsize[x]) maxsize[x]=size[v];
     }
 }
 void get_root(int x,int pa,int r)
 {int i;
   maxsize[x]=max(maxsize[x],size[r]-size[x]);
   if (maxsize[x]<minsize)
     {
       minsize=maxsize[x];
       root=x;
     }
   for (i=head[x];i;i=edge[i].next)
     {
       int v=edge[i].to;
       if (v==pa||vis[v]) continue;
       get_root(v,x,r);
     }
 }
 void pre_divide(int x)
 {int i;
   minsize=2e9;
   get_size(x,);
   get_root(x,,x);
   rt[++tot]=root;
   vis[root]=;
   for (i=head[root];i;i=edge[i].next)
     {
       int v=edge[i].to;
       if (vis[v]==)
     pre_divide(v);
     }
 }
 void pre(double x)
 {int i;
   for (i=;i<=n;i++)
     {
       edge[i*-].c=edge[i*-].d-x;
       edge[i*].c=edge[i*].d-x;
       vis[i]=;mx[i]=-2e9;
     }
 }
 bool pd(int root)
 {int i;
   int maxdep=;
   int q[],h=,t=,I,Q[];
   for (I=head[root];I;I=edge[I].next)
     {
       int v=edge[I].to;
       if (vis[v]) continue;
       h=;t=;
       fa[v]=root;dist[v]=edge[I].c;
       dep[v]=;
       q[]=v;
       while (h<t)
     {
       h++;
       int u=q[h];
       if (dep[u]==U) break;
       for (i=head[u];i;i=edge[i].next)
         {
           int v=edge[i].to;
           if (fa[u]!=v&&vis[v]==)
         {
           fa[v]=u;
           dep[v]=dep[u]+;
           dist[v]=dist[u]+edge[i].c;
           t++;
           q[t]=v;
         }
         }
     }
       int hh=,tt=,now=maxdep;
       for (i=;i<=t;i++)
     {
       int x=q[i];
       while (dep[x]+now>=L&&now>=)
         {
           while (hh<tt&&mx[Q[tt]]<mx[now]) tt--;
           tt++;Q[tt]=now;now--;
         }
       while (hh<tt&&Q[hh+]+dep[x]>U) hh++;
       if (hh<tt&&mx[Q[hh+]]+dist[x]>=) return ;
     }
           maxdep=max(maxdep,dep[q[t]]);
       for (i=;i<=t;i++)
         {
           fa[q[i]]=;
           if (dist[q[i]]>mx[dep[q[i]]]) mx[dep[q[i]]]=dist[q[i]];
         }
     }
   for (i=;i<=maxdep;i++)
     mx[i]=-2e9;
   return ;
 }
 bool check(int root,int &num)
 {int i;
   vis[root]=;
   if (pd(root)) return ;
   for (i=head[root];i;i=edge[i].next)
     {
       int v=edge[i].to;
       if (vis[v]==)
     {
       num++;
       if (check(rt[num],num)) return ;
     }
     }
   return ;
 }
 int main()
 {int i,u,v;
   double d,r,eps=1e-;
   cin>>n;
   cin>>L>>U;
   for (i=;i<=n-;i++)
     {
       scanf("%d%d%lf",&u,&v,&d);
       add(u,v,d);add(v,u,d);
       r=max(r,d);
     }
   pre_divide();
   double l=;
   while (r-l>eps)
     {
       double mid=(l+r)/2.0;
       pre(mid);
       int tot=;
       if (check(rt[],tot)) l=mid;
       else r=mid;
     }
   printf("%.3lf",(l+r)/2.0);
 }