BZOJ 1778: [Usaco2010 Hol]Dotp 驱逐猪猡 概率与期望+高斯消元

这个还挺友好的，自己相对轻松能想出来~
令 $f[i]$ 表示起点到点 $i$ 的期望次数，则 $ans[i]=f[i]\times \frac{p}{q}$

#include <cmath>   
#include <cstdio> 
#include <algorithm>  
#define N 305 
#define setIO(s) freopen(s".in","r",stdin) , freopen(s".out","w",stdout) 
using namespace std;  
int edges;  
double f[N][N];       
int deg[N],hd[N],to[N*N],nex[N*N]; 
void add(int u,int v) 
{ 
    nex[++edges]=hd[u],hd[u]=edges,to[edges]=v;  
}
void Gauss(int n)
{
    int i,j,k,now;
    for(i=1;i<=n;++i)
    {
        now=i;        
        for(j=i;j<=n;++j)
        {
            if(fabs(f[j][i])>fabs(f[now][i])) now=j; 
        }
        if(now!=i)
        {
            for(j=1;j<=n;++j) swap(f[i][j],f[now][j]);
        }
        if(f[i][i])
        {
            for(j=i+1;j<=n+1;++j) f[i][j]/=f[i][i];
            f[i][i]=1;
        }
        for(j=i+1;j<=n;++j)
        {
            double div=f[j][i];  
            for(k=i+1;k<=n+1;++k) f[j][k]-=div*f[i][k];                                        
            f[j][i]=0;
        }
    }
    for(i=n;i>=1;--i)
    {
        for(j=i+1;j<=n;++j)
        {
            f[i][n+1]-=f[j][n+1]*f[i][j];  
        }
    }
}        
int main() 
{ 
    // setIO("input"); 
    int n,m,p,q,i,j;   
    double in,out; 
    scanf("%d%d%d%d",&n,&m,&p,&q); 
    in=1.0*(double)(1.0*p/q),out=1.0-in;  
    for(i=1;i<=m;++i) 
    {
        int a,b; 
        scanf("%d%d",&a,&b),add(a,b),add(b,a),++deg[a],++deg[b]; 
    } 
    f[1][n+1]=1;   
    for(i=1;i<=n;++i) 
    {
        f[i][i]=1;   
        for(j=hd[i];j;j=nex[j]) 
        {
            int v=to[j];     
            f[i][v]=-((1.0/deg[v])*out);       
        }   
    } 
    Gauss(n);   
    for(i=1;i<=n;++i) 
    {
        printf("%.9f\n",f[i][n+1]*in); 
    }
    return 0; 
}