[BZOJ3143][HNOI2013]游走（概率期望+高斯消元）

题目：

我是超链接

题解：

与概率挂钩而且互相之间都有关系的我们考虑高斯消元（当然还有这么小的数据范围）
如果我们再设f[i]表示到i点的期望，那么由于权值不确定，我们有2n个未知数，然而并不能找到2n个方程
怎么办呢？我们设f[i]表示到i点的概率，不难发现$f[i]=\sum_{i,j相连}{f[j]\over du[j]}$
注意到n点就不会再走了，要特别处理
然后我们获得了n个方程可以解出f
那么对于一条边，我们经过ta的概率是$P(u,v)={f[u]\over du[u]}+{f[v]\over du[v]}$
然后我们按照边的概率从大到小赋值编号就好了，最后期望=概率*权值
有个细节就是最后要把f[n]设为0，因为到了n就不会往外走了
不是说好的无向简单连通图吗为什么这么多边。

代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=1000005;
int tot,nxt[N*2],point[505],v[N*2],du[505],x[N],y[N],n;double a[505] [505],b[505],f[505],p[N];
void addline(int x,int y)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x;
}
void gauss()
{
    for (int i=1;i<=n;i++)
    {
        int num=i;
        for (int j=i+1;j<=n;j++)
          if (a[num][i]<a[j][i]) num=j;
        if (num!=i)
        {
            for (int j=i;j<=n;j++)
              swap(a[num][j],a[i][j]);
            swap(b[num],b[i]);
        }
        for (int j=i+1;j<=n;j++)
        {
            double t=a[j][i]/a[i][i];
            for (int k=i;k<=n;k++) a[j][k]-=t*a[i][k];
            b[j]-=b[i]*t;
        }
    }
    for (int i=n;i>=1;i--)
    {
        f[i]=b[i]/a[i][i];
        for (int j=1;j<i;j++) b[j]-=f[i]*a[j][i];
    }
}
int main()
{
    int m;scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d",&x[i],&y[i]);
        addline(x[i],y[i]); du[x[i]]++; du[y[i]]++;
    }
    for (int xx=1;xx<=n;xx++)
    {
        a[xx][xx]=1.0;
        for (int i=point[xx];i;i=nxt[i])
          if (v[i]!=n) a[xx][v[i]]=-1.0/du[v[i]];
    }
    b[1]=1.0;
    gauss();f[n]=0;
    for (int i=1;i<=m;i++) p[i]=f[x[i]]/du[x[i]]+f[y[i]]/du[y[i]];
    sort(p+1,p+m+1);double ans=0;
    for (int i=1;i<=m;i++) ans+=p[i]*(m-i+1.0);
    printf("%.3lf",ans);
}