题目
Wayne 在玩儿一个很有趣的游戏。在游戏中,Wayne 建造了N 个城市,现在他想在这些城市间修一些公路,当然并不是任意两个城市间都能修,为了道路系统的美观,一共只有M 对城市间能修公路,即有若干三元组(Ui, Vi,Ci) 表示Ui 和Vi 间有一条长度为Ci 的双向道路。当然,游戏保证了,若所有道路都修建,那么任意两城市可以互相到达。
Wayne 拥有恰好N - 1 支修建队,每支队伍能且仅能修一条道路。当然,修建长度越大,修建的劳累度也越高,游戏设定是修建长度为C 的公路就会有C 的劳累度。当所有的队伍完工后,整个城市群必须连通,而这些修建队伍们会看看其他队伍的劳累情况,若劳累情况差异过大,可能就会引发骚动,不利于社会和谐发展。Wayne 对这个问题非常头疼,于是他想知道,这N - 1 支队伍劳累度的标准差最小能有多少。
分析
注意到N<=100,M <= 2000,数据很小。
于是,正解就是愉愉快快地大暴力。
不过,暴力也是没那么容易想出来滴。
首先,枚举这N-1个数的平均数p,因为要使
∑0<=i<n(ci−c¯)2n−1−−−−−−−−−−−−−−√
尽量小,显然要使
∑0<=i<n(ci−c¯)2
尽量小。那么以
|ci−p|
作为树的边权,做一遍最小生成树,并记录下最小生成树的边。
但是,我们知道最小生成树的边的
c¯
并不一定等于p。所以,把最小生成树的边重新求一遍平均值,再算答案。
真是又黄又暴力!
#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
using namespace std;
struct ddx
{
double f,g;
int x,y;
}a[12000];
int b[12000],n,m,d[12000],tot;
bool cmp(ddx x,ddx y)
{
return x.f<y.f;
}
int getfa(int x)
{
if(b[x]==x) return x;
int g=getfa(b[x]);
b[x]=g;
return g;
}
int main()
{ scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d%lf",&a[i].x,&a[i].y,&a[i].g);
}
double ans=maxlongint;
for(double i=0.1;i<=100;i+=0.25)
{
for(int j=1;j<=n;j++)
{
b[j]=j;
}
for(int j=1;j<=m;j++)
{
a[j].f=abs(a[j].g-i);
}
sort(a+1,a+m+1,cmp);
double sum=0;
tot=0;
for(int j=1;j<=m,tot<n-1;j++)
{
int x=getfa(a[j].x);
int y=getfa(a[j].y);
if(x!=y)
{
b[x]=y;
d[++tot]=j;
}
}
double g=0;
for(int j=1;j<=tot;j++)
{
g+=a[d[j]].g;
}
g/=n-1;
for(int j=1;j<=tot;j++)
{
sum+=(g-a[d[j]].g)*(g-a[d[j]].g);
}
ans=min(ans,sum);
}
printf("%.4lf\n",sqrt(ans/(n-1)));
}