【BZOJ 3143】 [Hnoi2013]游走 期望DP|高斯消元|贪心

每个 val表示 这个点经过几次

这样就明白了 我开始迷惑好长时间

#include <cstdio>
#include <cmath>
#include <iostream>
#include <algorithm>
using namespace std;
const int MAXN=500+10; 
const double EPS=1e-6;
double val[MAXN*MAXN],f[MAXN][MAXN],ans[MAXN*MAXN];
int n,m,g[MAXN][MAXN],d[MAXN],s[MAXN*MAXN],t[MAXN*MAXN];
void gauss()
{
	for(int i=1;i<n;i++)
	{
		int t=i;
		for(int j=i;j<n;j++)
			if(abs(f[i][j])>EPS){t=i;break;}
		if(i!=t)
			for(int j=i;j<n;j++) swap(f[i][j],f[t][j]);
		for(int j=i+1;j<n;j++)
		{
			double tmp=f[j][i]/f[i][i];
			for(int k=i;k<=n;k++)
				f[j][k]-=tmp*f[i][k];
		}
	}
	for(int i=n-1;i>=1;i--)
	{
		val[i]=f[i][n]/f[i][i];
		for(int j=i-1;j>=1;j--)
			f[j][n] -=val[i]*f[j][i];		
	}
} 
int main()
{
	cin >>n >>m;
	for(int i=1;i<=m;i++) 
	{
		scanf("%d %d",&s[i],&t[i]);
		g[s[i]][t[i]]++;
		g[t[i]][s[i]]++;
		d[s[i]]++;
		d[t[i]]++;
	}
	for(int i=1;i<n;i++)
		for(int j=1;j<n;j++)
		{
			if(i==j) f[i][j]=1;
			else f[i][j]=-1.0*g[j][i]/d[j]; 
		}
	f[1][n]=1; gauss();
	for(int i=1;i<=m;i++)
		ans[i]=val[s[i]]/d[s[i]]+val[t[i]]/d[t[i]];
	sort(ans+1,ans+m+1);
	double out=0;
	for(int i=1;i<=m;i++) out+=ans[i]*(m-i+1);
	printf("%.3lf",out);
	return 0;
} 

没注意内存re3发 也是没谁了