D.The Child and Zoo
Of course our child likes walking in a zoo. The zoo has n areas, that are numbered from 1 to n. The i-th area contains ai animals in it. Also there are m roads in the zoo, and each road connects two distinct areas. Naturally the zoo is connected, so you can reach any area of the zoo from any other area using the roads.
Our child is very smart. Imagine the child want to go from area p to area q. Firstly he considers all the simple routes from p to q. For each route the child writes down the number, that is equal to the minimum number of animals among the route areas. Let’s denote the largest of the written numbers as f(p, q). Finally, the child chooses one of the routes for which he writes down the value f(p, q).
After the child has visited the zoo, he thinks about the question: what is the average value of f(p, q) for all pairs p, q (p ≠ q)? Can you answer his question?
Input
The first line contains two integers n and m (2 ≤ n ≤ 105; 0 ≤ m ≤ 105). The second line contains n integers: a1, a2, …, an (0 ≤ ai ≤ 105). Then follow m lines, each line contains two integers xi and yi (1 ≤ xi, yi ≤ n; xi ≠ yi), denoting the road between areas xi and yi.
All roads are bidirectional, each pair of areas is connected by at most one road.
Output
Output a real number — the value of .
The answer will be considered correct if its relative or absolute error doesn’t exceed 10 - 4.
Examples
Input
4 3
10 20 30 40
1 3
2 3
4 3
Output
16.666667
Input
3 3
10 20 30
1 2
2 3
3 1
Output
13.333333
Input
7 8
40 20 10 30 20 50 40
1 2
2 3
3 4
4 5
5 6
6 7
1 4
5 7
Output
18.571429
思路:
两点间的直接连边的数字肯定是两点点权的较小值
先对所有边赋值为端点的小点权、
然后把边从大到小排序,
然后遍历,如果当前边的两端不在一个集合里面(并查集判断)
说明当前边权是两点间的最大值,说明两个连通块间的点的路径上的最大值就是当前边权
即题目所说的f(x,y),因为两个连通块中的点可以两两匹配
则答案累加num[x]*num[y]再乘上边权
因为点对无序,(f(x,y)和f(y,x))算作两对,
所以答案要乘上2,
最后取平均值输出就行了
ps:
感觉好多图的点对问题都用到了并查集,但是每次都有点想不到那边去。
还需要锻炼。
code:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define int long long
const int maxm=1e5+5;
struct Node{
int a,b,c;
}e[maxm];
int pre[maxm];
int num[maxm];
int v[maxm];
int n,m;
void init(){
for(int i=0;i<maxm;i++){
pre[i]=i;
num[i]=1;
}
}
int ffind(int x){
return pre[x]==x?x:pre[x]=ffind(pre[x]);
}
bool cmp(Node a,Node b){
return a.c>b.c;
}
signed main(){
init();
cin>>n>>m;
for(int i=1;i<=n;i++){
cin>>v[i];
}
for(int i=1;i<=m;i++){
cin>>e[i].a>>e[i].b;
e[i].c=min(v[e[i].a],v[e[i].b]);
}
sort(e+1,e+1+m,cmp);
int ans=0;
for(int i=1;i<=m;i++){
int a=ffind(e[i].a);
int b=ffind(e[i].b);
if(a!=b){
ans+=num[a]*num[b]*e[i].c;
pre[a]=b;
num[b]+=num[a];
}
}
printf("%.6f\n",ans*2.0/(n*(n-1)));
return 0;
}