BZOJ 3143 [Hnoi2013]游走（概率DP）

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。 

思路：首先肯定能想到贪心给贡献小的边大边权，给贡献小的边小边权，而边的贡献能有点的贡献得到：

即f[u][v]=ans[u]+ans[v]；

现在问题就转变成了求点的贡献了，而每个点的贡献又能通过每个点的期望求出：

即ans[i]=exp[i]/du[i]；

而对于每个点的期望又有：

i=1：

exp[1]=1+sigma(exp[j]/du[j]);（j与1有变相连）ps：因为1是入点所以期望+1。

i!=1：

exp[i]=sigma(exp[j]/du[j]);（j与i有变相连）

ps：由于n点无贡献 所以j不包含n。

稍微移项得到n-1元方程组：

高斯消元解方程组即可得到点的期望

再依次求得点的贡献->边的贡献->最终答案就很容易了

 ac代码：

#include <bits/stdc++.h>

using namespace std;

int n,m;
double Map[505][505];
double ans[505];
int G[505][505];
int u[250005],v[250005];
int du[505];
double f[250005];

void guess()
{
    for(int i=1;i<=n;i++)
    {
        int r=i;
        for(int j=i+1;j<=n;j++)
        {
            if(fabs(Map[r][i])<fabs(Map[j][i]))
            {
                r=j;
            }
        }
        if(i!=r) swap(Map[i],Map[r]);
        double div=Map[i][i];
        for(int j=i;j<=n+1;j++)
        {
            Map[i][j]/=div;
        }
        for(int j=i+1;j<=n;j++)
        {
            div=Map[j][i];
            for(int k=i;k<=n+1;k++)
            {
                Map[j][k]-=Map[i][k]*div;
            }
        }
    }
    ans[n]=Map[n][n+1];
    for(int i=n-1;i>=1;i--)
    {
        ans[i]=Map[i][n+1];
        for(int j=i+1; j<=n; j++)
        {
            ans[i]-=Map[i][j]*ans[j];
        }
    }
}

int main()
{
    cin >> n >> m;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&u[i],&v[i]);
        G[u[i]][v[i]]=1;
        G[v[i]][u[i]]=1;
        du[u[i]]++;
        du[v[i]]++;
    }
    for(int i=1;i<=n;i++)
    {
        if(i==1)
        {
            Map[i][n+1]=1;
        }
        else
        {
            Map[i][n+1]=0;
        }
        for(int j=1;j<=n;j++)
        {
            if(i==j)
            {
                Map[i][j]=1;
                continue;
            }
            if(j==n)
            {
                Map[i][j]=0;
                continue;
            }
            if(G[i][j])
            {
                Map[i][j]=-1.0/(double)du[j];
            }
        }
    }
    guess();
    for(int i=1;i<n;i++)
    {
        ans[i]/=du[i];
    }
    ans[n]=0;
    for(int i=1;i<=m;i++)
    {
        f[i]=ans[u[i]]+ans[v[i]];
    }
    sort(f+1,f+m+1);
    double sum=0;
    for(int i=1;i<=m;i++)
    {
        sum+=(m-i+1)*f[i];
    }
    printf("%.3f",sum);
    return 0;
}