[HNOI2013]游走 【概率dp+高斯消元】

题目链接：https://www.luogu.org/problem/P3232

参考题解：http://old.orzsiyuan.com/articles/problem-HNOI-2013-Walk/

首先考虑求每条边的经过次数的期望，给期望经过次数多的边分配小编号。

边的期望经过次数可以由点的经过次数推出。

设dp[i]为i节点的期望经过次数，out[i]为i节点的出度。

那么对于(u,v)这条边，它的期望经过次数就是dp[u]*out[u]+dp[v]*out[v]。

而点的经过期望可以通过高斯消元进行求解。

注意N号点不会继续对其他点/边产生贡献。

#include <bits/stdc++.h>
#define rep(i, a, b) for(ll i = (a); i <= (b); i++)
#define pb push_back
#define ll long long
#define all(x) (x).begin(),(x).end()
#define pii pair<int,int>
using namespace std;
const int N = 550;
int n,m;
const double eps = 1e-7;
double a[N][N],ans[N];
bool gauss(int n) {
    rep(i, 1, n) {          //逐列消除
        int p = i;
        rep(j, i+1, n)          //找到该列最大系数所在的行p,减少误差(也可以找第一个非0系数,但误差会大一些)
            if(fabs(a[p][i])<fabs(a[j][i]))
                p = j;

        if(fabs(a[p][i])<eps) return 0;    //如果没有找到,说明无解/无穷解
        swap(a[i],a[p]);        //否则把p交换上来,用其进行消元

        double div = a[i][i];
        rep(j, i, n+1) a[i][j] /= div;  //把p的第i列系数置1
        rep(j, 1, n)  {                 //开始消其他行
            if(i==j) continue;
            div = a[j][i];
            rep(k, i, n+1)
                a[j][k] -= a[i][k]*div;
        }
    }
    rep(i, 1, n) ans[i] = a[i][n+1];
    return 1;
}
vector<int> nxt[N];
pii edge[N*N];
int main() {
    //freopen("a.txt","r",stdin);
    int n,m;
    scanf("%d%d",&n,&m);
    rep(i, 1, m) {
        int u,v;
        scanf("%d%d",&u,&v);
        edge[i].first = u;
        edge[i].second = v;
        nxt[u].pb(v);
        nxt[v].pb(u);
    }
    rep(u, 1, n) {
        a[u][u] += 1;
        for(auto v:nxt[u]) {
            if(v==n) continue;
            a[u][v] -= 1.0/nxt[v].size();
        }
        if(u!=1) a[u][n+1] = 0;
        else a[u][n+1] = 1;
    }
    gauss(n);
    vector<double> k;
    rep(i, 1, m) {
        int u = edge[i].first;
        int v = edge[i].second;
        double p = 0;
        if(u!=n) p += ans[u]*1.0/nxt[u].size();
        if(v!=n) p += ans[v]*1.0/nxt[v].size();
        k.pb(p);
    }
    sort(all(k));
    double sum=0;
    for(auto x:k) sum += x*(m--);
    printf("%.3f",sum);
    return 0;
}