bzoj1778: [Usaco2010 Hol]Dotp 驱逐猪猡（概率DP+高斯消元）

　　深夜肝题。。。有害身心健康QAQ

　　设f[i]为到达i的概率，d[i]为i的度数。

　　因为无限久之后炸弹爆炸的概率是1，所以最后在i点爆炸的概率实际上就是f[i]/sigma(f[])

　　列出方程组 f[i]=sigma(f[to]*(1-p/q)/d[to]+[i==1]*(1-p/q))

　　然后就可以高斯消元了

　　高斯消元的方法：自己的那一位是1，to的每一位上为-(1-p/q)/d[to]，n+1位上为0，这样就相当于x减去所有to为0。1的n+1上为1-p/q，因为炸弹还可能在自己这里不跑。

　　这题实际上还有另一个思路，利用期望来计算，这里就直接贴博客了。。。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
const double eps=1e-;
struct poi{int too,pre;}e[maxn*maxn*];
int n,m,p,q,to,x,y,tot,d[maxn],last[maxn];
double sum,a[maxn][maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
void add(int x,int y){e[++tot].too=y;e[tot].pre=last[x];last[x]=tot;}
void gauss()
{
    for(int i=;i<=n;i++)
    {
        for(to=i;to<=n;to++)if(fabs(a[to][i])>eps)break;
        if(to!=i)for(int j=;j<=n+;j++)swap(a[i][j],a[to][j]);
        double x=a[i][i];for(int j=;j<=n+;j++)a[i][j]/=x;
        for(int j=;j<=n;j++)
        if(i!=j)
        {
            x=a[j][i];
            for(int k=;k<=n+;k++)
            a[j][k]-=x*a[i][k];
        }
    }
}
int main()
{
    read(n);read(m);read(p);read(q);double pro=1.0*p/q;
    for(int i=;i<=m;i++)
    {
        read(x);read(y);
        d[x]++;d[y]++;
        add(x,y);add(y,x);
    }
    for(int i=;i<=n;i++)
    {
        a[i][i]=;
        for(int j=last[i];j;j=e[j].pre)
        a[i][e[j].too]-=(-pro)/d[e[j].too];
    }
    a[][n+]=-pro;
    gauss();
    for(int i=;i<=n;i++)sum+=a[i][n+];
    for(int i=;i<=n;i++)printf("%.9lf\n",a[i][n+]/sum);
    return ;
}