BZOJ 3143: [Hnoi2013]游走（高斯消元+概率DP）

题面

https://www.lydsy.com/JudgeOnline/problem.php?id=3143

---

思路

期望的一种算法是先算概率，最后乘值， 这里我们先求期望次数，再编号。

显然的贪心就是给去的期望次数多的边编较小的编号。现在问题在于求经过每条边的期望次数。

我们求出点的，自然可以通过边两边的点求出边的期望。由于点数很小，可以用高斯消元来求。

随机任意游走，可以看做是稳定的，各个点的期望满足一个平衡关系。这就相当于列出一堆方程。$f[i]=\sum_{(i,j)\in E}\frac{f[j]}{deg[j]}$，直接解出每个$f$即可。

注意由于一开始在1号点，所以1的期望次数额外加1，n号点是终点没有贡献，所以n不出现在其他任何点的方程里，直接让$f[n]=0$。

将右边的$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;
}

---