算出每条边的期望走过次数,期望大的给小的编号,答案就最小
每条边的期望走过次数可以由两个端点的期望次数算出来
1个点的期望走过次数可以由与他相邻的点的期望表示出来
比较特殊的是n号点不会对任何点的期望产生恭喜,因为到n就停了,还有1号点的期望次数要+1因为他是起点
然后高斯消元就可以解出每个点的期望,然后每条边的期望也就能算出来了
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<vector>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
#include<map>
using namespace std;
#define MAXN 510
#define MAXM 1010
#define eps 1e-8
#define ll long long
#define INF 1000000000
#define MOD 1000000007
double a[MAXN][MAXN];
int n,m;
bool mp[MAXN][MAXN];
int d[MAXN];
double e[MAXN*MAXN];
int tot;
double ans;
void gs(){
int i,j,k;
for(i=1;i<n;i++){
if(fabs(a[i][i])<eps){
for(j=i+1;j<n;j++){
if(fabs(a[j][i])>eps){
for(k=1;k<=n+1;k++){
swap(a[i][k],a[j][k]);
}
break;
}
}
}
if(fabs(a[i][i])>eps){
for(j=1;j<n;j++){
if(i!=j){
double t=a[j][i]/a[i][i];
for(k=1;k<=n+1;k++){
a[j][k]-=a[i][k]*t;
}
}
}
}
}
}
int main(){
int i,j,x,y;
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++){
scanf("%d%d",&x,&y);
d[x]++;
d[y]++;
mp[x][y]=mp[y][x]=1;
}
for(i=1;i<n;i++){
a[i][i]=1;
for(j=1;j<n;j++){
if(mp[i][j]){
a[i][j]-=1.0/d[j];
}
}
}
a[1][n+1]=1;
gs();
for(i=1;i<n;i++){
a[i][n+1]/=a[i][i];
}
for(i=1;i<=n;i++){
for(j=i+1;j<=n;j++){
if(mp[i][j]){
e[++tot]=a[i][n+1]/d[i]+a[j][n+1]/d[j];
}
}
}
sort(e+1,e+tot+1);
for(i=1;i<=m;i++){
ans+=e[i]*(m-i+1);
}
printf("%.3lf\n",ans);
return 0;
}
/*
2 1
1 2
*/