【HNOI2013】游走

Description

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

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

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

Input

输入第一行是正整数N和M，分别表示该图的顶点数和边数，接下来M行每行是整数u，v(1≤u,v≤N)，表示顶点u与顶点v之间存在一条边。

Output

输出仅包含一个实数，表示最小的期望值，保留3 位小数。

Sample Input

3 3

2 3

1 2

1 3

Sample Output

3.333

Data Constraint

30%的数据满足N≤10；

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

Hint

样例解释：边(1,2)编号为1，边(1,3)编号2，边(2,3)编号为3。

Solution

设DP,f[i]表示i号边的期望次数，但是发现很难转移（转移不了）
那么变成f[i]表示i号点的期望次数，那么

$$f[i]=\sum _{j与i相连} {{f[j]}\over r[j]}$$
r[j]为点j的度数，但是怎么求呢？
用高斯消元
其中f[n]=0因为到n就结束，f[n]不能给其它点贡献
f[1]在开始前的答案就为1，因为开始点在1号

Code

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 510
#define db double
using namespace std;
int n,m,d[N],c[N][N],q[N*1000][2];
db a[N][N],b[N],ans[N*1000];
void jh(int i,int j)
{
    fo(k,1,n) swap(a[i][k],a[j][k]);
}
void gaosi()
{
    fo(k,1,n-1) 
    {
        int l=k;
        fo(i,k+1,n-1) if(abs(a[i][k])>abs(a[l][k])) l=i;
        jh(l,k);
        fo(i,k+1,n-1)
        {
            if(a[i][k]==0) continue;
            b[i]-=b[k]*a[i][k]/a[k][k];
            fd(j,n,k) a[i][j]-=a[k][j]*a[i][k]/a[k][k];
        }
    }
    for(int k=n;k;k--)
    {
        fo(i,k+1,n) b[k]=b[k]-b[i]*a[k][i];
        b[k]/=a[k][k];
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,m)
    {
        int x,y;scanf("%d%d",&x,&y);q[i][0]=x;q[i][1]=y;
        d[x]++;d[y]++;
        c[x][y]++;c[y][x]++;
    }
    fo(i,1,n-1)
    {
        a[i][i]=-1,b[i]=0;
        fo(j,1,n)
        if(d[j]!=0) a[i][j]+=(1.0*c[i][j])/(d[j]);
    }
    a[n][n]=b[1]=-1;
    gaosi();
    fo(i,1,m)
    if(d[q[i][0]]!=0&&d[q[i][1]]!=0) 
    ans[i]=b[q[i][0]]/d[q[i][0]]+b[q[i][1]]/d[q[i][1]];
    sort(ans+1,ans+m+1);
    db an=0;
    fo(i,1,m) an+=ans[i]*(m-i+1);
    printf("%.3lf",an);
}