【bzoj3143】[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之间存在一条边。 输入保证30%的数据满足N≤10，100%的数据满足2≤N≤500且是一个无向简单连通图。

Output

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

Sample Input

3 3

2 3

1 2

1 3

Sample Output

3.333
HINT

边(1,2)编号为1，边(1,3)编号2，边(2,3)编号为3。

题解
http://blog.csdn.net/vmurder/article/details/44542575
Orz dalao

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
#define inf 100000000
#define mod 10000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
double a[505][505],ANS,ans[505],w[250005];
int du[505],n,m;
struct node{int x,y;}e[250005];
void Gauss(int n)
{
    for (int i=1;i<n;i++)
    {
        int k=0;
        for (int j=i;j<=n;j++)
            if (fabs(a[j][i])>fabs(a[k][i])) k=j;
        for (int j=i;j<=n+1;j++) swap(a[i][j],a[k][j]);
        for (int j=i+1;j<=n;j++)
        {
            double tmp=-a[j][i]/a[i][i];
            for (int k=i;k<=n+1;k++)
                a[j][k]+=a[i][k]*tmp;
        }
    }
    for (int i=n;i;i--)
    {
        for (int j=n;j>i;j--)
            a[i][n+1]-=a[i][j]*ans[j];
        ans[i]=a[i][n+1]/a[i][i];
    }
}
int main()
{
    n=read();m=read();
    for (int i=1;i<=m;i++)
    {
        e[i].x=read();e[i].y=read();
        du[e[i].x]++;du[e[i].y]++;
    }
    for (int i=1;i<n;i++)a[i][i]=-1;
    for (int i=1;i<=m;i++)
    {
        a[e[i].x][e[i].y]+=1.0/du[e[i].y];
        a[e[i].y][e[i].x]+=1.0/du[e[i].x];
    }
    for (int i=1;i<=n;i++)a[n][i]=0;
    a[1][n+1]=-1;a[n][n]=1;
    Gauss(n);
    for (int i=1;i<=m;i++) w[i]=ans[e[i].x]/du[e[i].x]+ans[e[i].y]/du[e[i].y];
    sort(w+1,w+m+1);
    for (int i=1;i<=m;i++) ANS+=(m-i+1)*w[i];
    printf("%.3lf",ANS);
    return 0;
}