【GDOI2014模拟】Tree 题解+代码

Description

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

Input

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

Output

输出仅一行表示最小的标准差，保留4 位小数。

Sample Input

3 3
1 2 1
2 3 2
3 1 3

Sample Output

0.5000

Data Constraint

对于20% 的数据，M <= 20。
对于另外30% 的数据，Ci <= 10。
对于100% 的数据，N <= 100，M <= 2000，Ci <= 100。

Solution

简化题目，变成了求一个最小标准差生成树。
标准差是什么？就是反应一堆数据的波动大小的。它不是每个数据都要与平均数相减吗？而且才小于100.那就枚举一个平均数呗。把每条边的值与枚举出来的这个数相减， 求一个最小生成树即可。但是这个枚举出的平均数不一定是正确的，那就在选出边之后再重新求答案就行了。
平均数不能每次+1枚举，+0.25或者+0.1就行了。证明自己想。小心被卡精度

code

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define N 2010
#define ll long long
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define sqr(x) ((x)*(x))
using namespace std;
int n,bz[N],f[N],m;
struct note{
    int x,y,d;
    double b;
};
note a[N];
bool cnt(note x,note y)
{
    return x.b<y.b;
}
int gf(int x)
{
    if (f[x]==0) return x;
    f[x]=gf(f[x]);return f[x];
}
double sb()
{
    memset(f,0,sizeof(f));
    fo(i,1,m)
    {
        bz[i]=0;
        int x=gf(a[i].x),y=gf(a[i].y);
        if (x!=y) bz[i]=1,f[x]=y;
    }
    double p=0;
    fo(i,1,m) if (bz[i]) p+=a[i].d;
    p/=(n-1);
    double ans=0;
    fo(i,1,m) if (bz[i]) ans+=sqr(a[i].d-p);
    ans=ans/(n-1);
    return sqrt(ans);
}
int main()
{
    scanf("%d%d",&n,&m);int mx=0;double ans=1234567;
    fo(i,1,m) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].d),mx=max(mx,a[i].d);
    for(double i=0;i<=mx;i+=0.25)
    {
        fo(j,1,m) a[j].b=abs(a[j].d-i);sort(a+1,a+m+1,cnt);
        ans=min(ans,sb());
    }
    printf("%.4f\n",ans);
}