bzoj3754 Tree之最小方差树

http://www.elijahqi.win/archives/2825
Description

Wayne在玩儿一个很有趣的游戏。在游戏中，Wayne建造了N个城市，现在他想在这些城市间修一些公路，当然并不是任意两个城市间都能修，为了道路系统的美观，一共只有M对城市间能修公路，即有若干三元组 (Ui,Vi,Ci)表示Ui和Vi间有一条长度为Ci的双向道路。当然，游戏保证了，若所有道路都修建，那么任意两城市可以互相到达。Wayne拥有恰好N-1支修建队，每支队伍能且仅能修一条道路。当然，修建长度越大，修建的劳累度也越高，游戏设定是修建长度为C的公路就会有C的劳累度。当所有的队伍完工后，整个城市群必须连通，而这些修建队伍们会看看其他队伍的劳累情况，若劳累情况差异过大，可能就会引发骚动，不利于社会和谐发展。Wayne对这个问题非常头疼，于是他想知道，这N1支队伍劳累度的标准差最小能有多少。
标准差的定为：设有N个数，分别为ai,它们的平均数为 ,那么标准差就是

Input

第一行两个正整数N,M
接下来M行，每行三个正整数Ui,Vi,Ci

Output

输出最小的标准差，保留四位小数。
Sample Input

3 3
1 2 1
2 3 2
3 1 3
Sample Output

0.5000
HINT

N<=100,M<=2000,Ci<=100
因为 c的范围非常小 所以可以求一下总共最小 还有总共最大 然后直接枚举总和 然后求一下最小的方差生成树 如果可以更新答案需要验证一下 是否和我枚举的答案一样

#include<cmath>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
int n,m,t;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
namespace sol1{
    const int N=2200;
    struct node{
        int x,y,z;double w1;
    }data[N];
    int fa[N];double ans;
    inline void init(){
        for (int i=1;i<=n;++i) fa[i]=i;
    }
    inline double sqr(double x){return x*x;}
    inline int find(int x){
        while(x!=fa[x]) x=fa[x]=fa[fa[x]];return x;
    }
    inline bool cmp3(const node &a,const node &b){return a.w1<b.w1;}
    inline void kruskal(int x){
        init();double ave=(double) x/(double)(n-1);
        for (int i=1;i<=m;++i) data[i].w1=sqr((double)data[i].z-ave);
        sort(data+1,data+m+1,cmp3);
        int all=0;double all1=0;int tot=0;
        for (int i=1;i<=m;++i){
            int x=find(data[i].x),y=find(data[i].y);
            if (x!=y){
                fa[x]=y;all+=data[i].z;all1+=data[i].w1;if (++tot==n-1) break;
            }
        }
        if (all==x) ans=min(ans,all1);
    }
    inline bool cmp1(const node &a,const node &b){
        return a.z<b.z;
    }
    inline bool cmp2(const node &a,const node &b){
        return a.z>b.z;
    }
    inline void gao(){
        ans=999999999999999999999999999999.0;
        for (int i=1;i<=m;++i) data[i].x=read(),data[i].y=read(),data[i].z=read();
        sort(data+1,data+m+1,cmp1);int sum1=0;
        for (int i=1;i<n;++i) sum1+=data[i].z;
        sort(data+1,data+m+1,cmp2);int sum2=0;
        for (int i=1;i<n;++i) sum2+=data[i].z;
        for (int i=sum1;i<=sum2;++i) kruskal(i);
        printf("%.4f\n",sqrt(ans/(double)(n-1)));
    }
}
int main(){
    freopen("t2.in","r",stdin);
    n=read();m=read();
    sol1::gao();
    return 0;
}