BZOJ3143 [Hnoi2013]游走

算出每条边的期望走过次数，期望大的给小的编号，答案就最小

每条边的期望走过次数可以由两个端点的期望次数算出来

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

*/