题面
https://www.lydsy.com/JudgeOnline/problem.php?id=3143
思路
期望的一种算法是先算概率,最后乘值, 这里我们先求期望次数,再编号。
显然的贪心就是给去的期望次数多的边编较小的编号。现在问题在于求经过每条边的期望次数。
我们求出点的,自然可以通过边两边的点求出边的期望。由于点数很小,可以用高斯消元来求。
随机任意游走,可以看做是稳定的,各个点的期望满足一个平衡关系。这就相当于列出一堆方程。 f[i]=∑(i,j)∈Ef[j]deg[j] ,直接解出每个 f 即可。
注意由于一开始在1号点,所以1的期望次数额外加1,n号点是终点没有贡献,所以n不出现在其他任何点的方程里,直接让
将右边的 f[i] 挪到左边使方程其次,只有 f[1] 有一个-1的常数。然后直接高斯消元。无向简单连通图好像没有自环和重边,不用考虑。
求出边的期望,然后排序标号即可。
代码
#include <bits/stdc++.h>
#define maxn 505
#define Eps 1e-10
using namespace std;
int n, m, D[maxn];
double A[maxn][maxn], tata;
struct ZWT{
int x, y;
double p;
bool operator < (const ZWT& OTHER) const{
return (OTHER.p - p) < -Eps;
}
}E[maxn*maxn];
void Gauss(){
for(int i = 1; i <= n; i++){
int x = i;
for(int j = i+1; j <= n; j++)
if(fabs(A[j][i]) > fabs(A[x][i])) x = j;
swap(A[x], A[i]);
for(int j = n+1; j >= i; j--) A[i][j] /= A[i][i];
for(int j = 1; j <= n; j++){
if(i == j) continue;
for(int k = n+1; k >= i; k--) A[j][k] -= A[j][i] * A[i][k];
}
}
}
int main(){
scanf("%d%d", &n, &m);
for(int i = 1; i <= m; i++){
scanf("%d%d", &E[i].x, &E[i].y);
D[E[i].x] ++;
D[E[i].y] ++;
}
for(int i = 1; i <= m; i++){
if(E[i].y != n) A[E[i].y][E[i].x] = 1.0 / D[E[i].x];
if(E[i].x != n) A[E[i].x][E[i].y] = 1.0 / D[E[i].y];
}
for(int i = 1; i <= n; i++) A[i][i] = -1.0;
A[1][n+1] = -1.0;
Gauss();
for(int i = 1; i <= m; i++)
E[i].p = A[E[i].x][n+1] / D[E[i].x] + A[E[i].y][n+1] / D[E[i].y];
sort(E+1, E+m+1);
for(int i = 1; i <= m; i++) tata += E[i].p * i;
printf("%.3lf\n", tata);
return 0;
}