【HNOI2013】游走

Description

一个无向连通图，顶点从1 编号到N，边从1 编号到M。

小Z 在该图上进行随机游走，初始时小Z 在1 号顶点，每一步小Z 以相等的概率随机选择当前顶点的某条边，沿着这条边走到下一个顶点，获得等于这条边的编号的分数。当小Z到达N 号顶点时游走结束，总分为所有获得的分数之和。

现在，请你对这M 条边进行编号，使得小Z 获得的总分的期望值最小。

100%的数据满足2≤N≤500且是一个无向简单连通图。

Analysis

贪心地对期望经过次数最多的边编号最小嘛，这个贪心显然正确
然后只需求出每条边的期望经过次数
边的期望好像不是很好求，但是可以转化成求点的期望
设$p[e(u,v)]$表示其期望经过次数，则$p[e(u,v)]=p_u/d_u+p_v/d_v$
$d_i$表示点 i <script type="math/tex" id="MathJax-Element-101">i</script>的度数
然后很容易得出一个方程组，用高斯消元解即可

Code

#include<cstdio>
#include<iostream>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,b,a) for(int i=b;i>=a;i--)
using namespace std;
const int N=510,M=N*N;
int tot,n,m,to[M],next[M],last[N],d[N],x[M],y[M];
double a[N][N],b[N],c[M];
void link(int u,int v)
{
    to[++tot]=v;
    next[tot]=last[u];
    last[u]=tot;
    d[u]++;
}
void Gauss()
{
    fo(i,1,n)
    {
        fo(j,i+1,n) a[i][j]/=a[i][i];
        b[i]/=a[i][i];
        fo(j,i+1,n)
        {
            fo(k,i+1,n) a[j][k]-=a[j][i]*a[i][k];
            b[j]-=a[j][i]*b[i];
        }
    }
    fd(i,n,1)
        fo(j,1,i-1) b[j]-=a[j][i]*b[i];
}
int main()
{
    scanf("%d %d",&n,&m);
    fo(i,1,m)
    {
        scanf("%d %d",&x[i],&y[i]);
        link(x[i],y[i]),link(y[i],x[i]);
    }
    fo(i,1,n-1)
    {
        a[i][i]=1;
        for(int j=last[i];j;j=next[j]) a[i][to[j]]=-1.0/d[to[j]];
    }
    b[1]=a[n][n]=1;
    Gauss();
    fo(i,1,m)
    c[i]=b[x[i]]/d[x[i]]+b[y[i]]/d[y[i]];
    sort(c+1,c+m+1);
    double ans=0;
    fo(i,1,m) ans+=(m-i+1)*c[i];
    printf("%.3lf",ans); 
    return 0;
}