记录编号	312599	评测结果	AAAAAAAAAA
题目名称	[HNOI 2013]游走	最终得分	100
用户昵称	TenderRun	是否通过	通过
代码语言	C++	运行时间	3.826 s
提交时间	2016-09-26 21:21:49	内存使用	7.27 MiB

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
const int N=510,M=N*N;
int n,m,e[M][2],d[N];
long double A[N][N];
double ans,w[M],x[N];
void Guass_Elimination(){
	for(int i=1;i<n;i++){
		int p=i;
		for(int j=i+1;j<n;j++)
			if(fabs(A[p][i])<fabs(A[j][i]))
				p=j;
		if(p!=i)
			for(int j=1;j<=n;j++)
				swap(A[i][j],A[p][j]);
		long double tmp=A[i][i];
		for(int j=1;j<=n;j++)
			A[i][j]/=tmp;
		for(int j=1;j<n;j++)
			if(i!=j){
				tmp=A[j][i];
				for(int k=1;k<=n;k++)
					A[j][k]-=tmp*A[i][k];
			}
	}
	for(int i=1;i<n;i++)
		x[i]=A[i][n];
}
int main(){
	freopen("walk.in","r",stdin);
	freopen("walk.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d",&e[i][0],&e[i][1]);
		d[e[i][0]]+=1;d[e[i][1]]+=1;
	}A[1][n]=1;
	for(int i=1;i<n;i++)A[i][i]=1;
	for(int i=1;i<=m;i++){
		if(e[i][0]==n||e[i][1]==n)continue;
		A[e[i][0]][e[i][1]]+=-1.0/d[e[i][1]];
		A[e[i][1]][e[i][0]]+=-1.0/d[e[i][0]];
	}
	Guass_Elimination();
	for(int i=1;i<=m;i++){
		w[i]+=x[e[i][0]]/d[e[i][0]];
		w[i]+=x[e[i][1]]/d[e[i][1]];
	}
	sort(w+1,w+m+1);
	for(int i=1;i<=m;i++)
		ans+=w[i]*(m-i+1);
	printf("%.3lf\n",ans);	
	return 0;
}