bzoj3143: [Hnoi2013]游走 高斯消元求概率

设p[i]为点i经过次数的期望.

d[i]为点i的度数

这个对于1<i<n p[i]=sigma(p[j]/d[j]) (i,j之间有连边)  

特殊地，p[1]=1+sigma(p[j]/d[j]) (i,j之间有连边) 

而虽然n点会被经过好多次，但是n点对它连接的边是不会有贡献的，所以p[n]=0。

由于是连通图，必然有唯一解.

高斯消元解即可.

对于每一条边(v，u)它的期望经过次数为p[v]/d[i]+p[u]/d[v]

然后贪心赋边权，期望越大的边赋越小的值。。

竟然要把eps开到1e-9.。。。。。。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
#define eps 1e-10
double a[600][600],p[600];
int n,m,e,in[600],g[600][600];
struct node
{
    int u,v;
    double val;
}save[300000];
void guass()
{
    m=n+1;
    int row,col,i,j,k;
    for(row=1,col=1;row<n,col<m;row++,col++)
    {
        k=row;
        for(int i=row+1;i<=n;i++) if(fabs(a[i][col])>fabs(a[k][col])) k=i;

        if(k!=row)
        {
            for(i=col; i<=m; i++)  swap(a[k][i],a[row][i]);
        }
        for(i=row+1; i<=n; i++)
        {
            if(fabs(a[i][col])<eps)
                continue;
            double t=a[i][col]/a[row][col];
            a[i][col]=0.0;
            for(int j=col+1;j<=m;j++)
                a[i][j]-=a[row][j]*t;
        }
    }
    for(i=n;i>=1;i--)
    {
        p[i]=a[i][m];
        for(j=i+1;j<=n;j++)
            p[i]-=p[j]*a[i][j];
        p[i]/=a[i][i];
    }
}
bool cmp(node aa,node bb)
{
    return aa.val>bb.val;
}
int main()
{
    scanf("%d%d",&n,&e);
    int u,v;
    for(int i=1;i<=e;i++)
    {
        scanf("%d%d",&u,&v);
        save[i].u=u;save[i].v=v;
        in[u]++;
        in[v]++;
        g[u][in[u]]=v;
        g[v][in[v]]=u;
    }
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<=in[i];j++)
        {
            a[i][g[i][j]]=1.0/double(in[g[i][j]]);
        }
        a[i][i]=-1.0;
    }
    a[1][n+1]=-1.0;
    a[n][n]=1.0;
    guass();
    for(int i=1;i<=e;i++)
    {
        save[i].val=(p[save[i].v]/double(in[save[i].v]))+p[save[i].u]/double(in[save[i].u]);
    }
    sort(save+1,save+1+e,cmp);
    double ans=0;
    for(int i=1;i<=e;i++)
    {
        ans+=save[i].val*double(i);
    }
    printf("%.3f",ans);
    return 0;
}