【xsy1303】生成树 乱搞

题目大意：给你n个点，问这n个点构成的最小标准差生成树的值

这题题解里面写的都是什么代码？？

你用O(n^4)的复杂度，枚举出两条边，然后求平均数。

对于剩下的边，我们求出这些边与平均数的差，然后求最下差生成树。

然后就过了？？？

不是很懂

 1 #include<bits/stdc++.h>
 2 #define M 22
 3 #define sqr(x) ((x)*(x))
 4 using namespace std;
 5 
 6 double x[M]={0},y[M]={0},dis[M][M]={0}; int n;
 7 
 8 int m=0;
 9 
10 struct edge{
11     int u,v; double w;
12     friend bool operator <(edge a,edge b){return a.w<b.w;}
13 }a[M*M]; bool sel[M*M]={0};
14 
15 int f[M]={0}; int get(int x){return x==f[x]?x:f[x]=get(f[x]);}
16 
17 map<double,int> mp; 
18 
19 double solve(double avg){
20     if(mp[avg]) return 459354345739875984; mp[avg]=1;
21     memset(sel,0,sizeof(sel));
22     for(int i=1;i<=n;i++) f[i]=i;
23     m=0;
24     for(int i=1;i<=n;i++)
25     for(int j=i+1;j<=n;j++)
26     a[++m]={i,j,sqr(dis[i][j]-avg)};
27     sort(a+1,a+m+1);
28     double sum=0;
29     for(int i=1;i<=m;i++){
30         int U=get(a[i].u),V=get(a[i].v);
31         if(U==V) continue;
32         sel[i]=1;
33         f[U]=V;
34     }
35     for(int i=1;i<=m;i++) if(sel[i]){
36         sum+=dis[a[i].u][a[i].v];
37     }
38     avg=sum/(n-1); sum=0;
39     for(int i=1;i<=m;i++) if(sel[i]){
40         sum+=sqr(avg-dis[a[i].u][a[i].v]);
41     }
42     return sum;
43 }
44     
45 
46 int Main(){
47     scanf("%d",&n);
48     for(int i=1;i<=n;i++) scanf("%lf",x+i);
49     for(int j=1;j<=n;j++) scanf("%lf",y+j);
50     for(int i=1;i<=n;i++)
51     for(int j=1;j<=n;j++)
52     dis[i][j]=sqrt(sqr(x[i]-x[j])+sqr(y[i]-y[j]));
53     
54     double minn=123456789000000;
55     for(int i1=1;i1<=n;i1++)
56     for(int j1=1;j1<=n;j1++) if(i1!=j1){
57         double dis1=dis[i1][j1];
58         for(int i2=1;i2<=n;i2++)
59         for(int j2=1;j2<=n;j2++) if(i2!=j2){
60             double dis2=dis[i2][j2];
61             minn=min(minn,solve((dis1+dis2)/2));
62         }
63     }
64     printf("%.3lf\n",sqrt(minn/(n-1)));
65 }
66 int main(){
67     int cas; cin>>cas;
68     while(cas--) Main();
69 }