BZOJ3143 [Hnoi2013]游走

本文版权归ljh2000和博客园共有，欢迎转载，但须保留此声明，并给出原文链接，谢谢合作。

本文作者：ljh2000
作者博客：http://www.cnblogs.com/ljh2000-jump/
转载请注明出处，侵权必究，保留最终解释权！

 

题目链接：BZOJ3143

 

正解：高斯消元

解题报告：

　　因为直接考虑每条边被经过的次数很不好考虑，而对于点统计次数就方便很多了。

　　令$p[i]$为从$i$出发开始游走的期望次数，则$p[u]=sum{frac{p[v]}{d[v]}},u,v$连通,$d$是点的度数。那么我可以得到一个$n$元$1$次方程组，注意$p[1]=1+sum{frac{p[v]}{d[v]}},1,v$连通,$n$不能走出去,所以无需计算贡献，然后高斯消元解出每个$p$即可。算出$p$之后，一条边的贡献就是所连接的两个点的$p/d$（拆成两个方向的角度很容易想清楚）。我开始想优化精度，就在高斯消元里面加了个$swap$操作发现除以$0$萎掉了…

　　以后看来不能随便加这种优化了...

 

//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);
const int MAXN = 520;
const int MAXM = 250011;
int n,m,d[MAXN];
bool mp[MAXN][MAXN];
LB tot,ans[MAXM],p[MAXN],f[MAXN][MAXN];//p[i]表示从i出发的期望次数
struct edge{ int x,y; }e[MAXM];
inline bool cmp(LB q,LB qq){ return q>qq; }
inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

inline void gauss(){
	LB now;
	for(int i=1;i<=n;i++) {
		for(int j=i+1;j<=n;j++) {
			now=f[j][i]/f[i][i];
			for(int l=i;l<=n+1;l++)
				f[j][l]-=now*f[i][l];
		}
	}

	for(int i=n;i>=1;i--) {
		p[i]=f[i][n+1];
		for(int j=i+1;j<=n;j++)
			p[i]-=p[j]*f[i][j];
		p[i]/=f[i][i];//还需解一个一元一次方程
	}
}

inline void work(){
	n=getint(); m=getint(); int x,y;
	for(int i=1;i<=m;i++) { e[i].x=x=getint(); e[i].y=y=getint(); mp[x][y]=mp[y][x]=1; d[x]++; d[y]++; }
	for(int i=1;i<n;i++) {//n走不出去...
		for(int j=1;j<=n;j++)
			if(mp[i][j])
				f[i][j]=1.0/d[j];
		f[i][i]=-1;
	}
	f[1][n+1]=-1;//p[1]=1+sigma(p[v])
	f[n][n]=1;

	gauss();

	for(int i=1;i<=m;i++) ans[i]=p[e[i].x]/d[e[i].x]+p[e[i].y]/d[e[i].y];
	sort(ans+1,ans+m+1,cmp);
	for(int i=1;i<=m;i++) tot+=ans[i]*i;
	printf("%.3Lf",tot);
}

int main()
{
    work();
    return 0;
}