码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ3754】Tree之最小方差树

时间:2015-08-12 14:42:10      阅读:129      评论:0      收藏:0      [点我收藏+]

标签:最小生成树   方差   乱搞   

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

1 2 1

2 3 2

3 1 3
Sample Output

0.5000
HINT

N<=100,M<=2000,Ci<=100

Source

都知道方差是衡量一组数据均匀不均匀的
我们把方差看成一些点到一条平行于x轴的直线的距离的平方
这时候可以想到去枚举这条直线的y坐标,范围当然跟c的范围一样
然后每枚举到一个值对这个值和边权记一下差的绝对值
按这个绝对值排序然后做最小生成树,记下边权和平均值求方差
显然这样得到的方案比较优
要注意的是枚举平行于x轴直线坐标时候不能简单1-100只取整数,要每次加0.5/0.1/0.25/0.2之类的要不会WA.(反正我取加1WA了之后直接改成了0.1)
Code:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<queue>
#define MAXN 2010
using namespace std;
int n,m;
double minn=1e20;
/*
    跪TA!
    跪TA!
    跪TA!
    跪TA!
    跪TA!
    跪TA! 
*/
double sum,Sum;
bool vis[MAXN];
struct edge
{
    int u,v,w;
    double delta;
    bool operator<(const edge& a)const
    {
        if (delta!=a.delta) return delta<a.delta;
        return w<a.w;
    }
}e[MAXN];
int f[MAXN];
void in(int &x)
{
    char ch=getchar();x=0;
    while (!(ch>=‘0‘&&ch<=‘9‘)) ch=getchar();
    while (ch>=‘0‘&&ch<=‘9‘)    x=x*10+ch-‘0‘,ch=getchar();
}
int find(int x)
{
    if (f[x]==x)    return x;
    return f[x]=find(f[x]);
}
void Union(int x,int y)
{
    int a=find(x),b=find(y),flag=rand()%2;
    if (flag)   f[a]=b;
    else    f[b]=a;
}
int main()
{
    freopen("loliconcon.in","r",stdin);
    freopen("loliconcon.out","w",stdout);
    in(n);in(m);
    for (int i=1;i<=m;i++)  in(e[i].u),in(e[i].v),in(e[i].w);
    for (double i=1;i<=100;i+=0.25)
    {
        sum=0;Sum=0;
        memset(vis,0,sizeof(vis));
        for (int j=1;j<=n;j++)  f[j]=j;//init
        for (int j=1;j<=m;j++)  e[j].delta=abs(e[j].w-i);
        sort(e+1,e+m+1);
        for (int j=1;j<=m;j++)//Kruskal
            if  (find(e[j].u)!=find(e[j].v))
            {
                Union(e[j].u,e[j].v);
                sum+=e[j].w;
                vis[j]=1;
            }
        sum/=(double)(n-1);
        for (int j=1;j<=m;j++)
            if (vis[j]) Sum+=(e[j].w-sum)*(e[j].w-sum);
        Sum/=(double)(n-1);Sum=sqrt(Sum);
        minn=min(minn,Sum);
    }
    printf("%.4f\n",minn);
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

【BZOJ3754】Tree之最小方差树

标签:最小生成树   方差   乱搞   

原文地址:http://blog.csdn.net/creationaugust/article/details/47443755

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!