CodeForces 113D

dp+高斯消元
两个人从a，b同时出发，随机走，当他们相遇的时候停止，问在A点停止的概率

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
const double eps=1e-14;
const int N=25*25;
int n,m,s,t;
int mp[N][N];
int g[N][N];
int tot;
int du[N];
double pp[N];

int getid(int i,int j){
    if(i>j)swap(i,j);
    return mp[i][j];
}

double A[N][N];
void getA(){
    tot=0;
    for(int i=0;i<n;i++){
        for(int j=i;j<n;j++){
            mp[j][i]=mp[i][j]=tot++;
        }
    }
    memset(A,0,sizeof(A));
    for(int i=0;i<n;i++){
        A[mp[i][i]][mp[i][i]]=1;
        for(int j=i+1;j<n;j++){
            A[mp[i][j]][mp[i][j]]+=1-pp[i]*pp[j];
            for(int k=0;k<n;k++)if(g[i][k]){
                int x=getid(i,j);
                int y=getid(k,j);
                A[x][y]-=(1-pp[i])/du[i]*pp[j];
            }
            for(int k=0;k<n;k++)if(g[j][k]){
                int x=getid(i,j);
                int y=getid(i,k);
                A[x][y]-=pp[i]*(1-pp[j])/du[j];
            }
            for(int k=0;k<n;k++)if(g[i][k]){
                for(int k1=0;k1<n;k1++)if(g[j][k1]){
                    int x=getid(i,j);
                    int y=getid(k,k1);
                    A[x][y]-=(1-pp[i])/du[i]*(1-pp[j])/du[j]
                        *g[i][k]*g[j][k1];
                }
            }
        }
    }
}


double B[N][N];
void getB(){
    memset(B,0,sizeof(B));
    for(int i=0;i<n;i++){
        B[mp[i][i]][i]=1;
    }
}

double X[N][N];
void Gauss(){
    int r=0;
    for(int i=0;i<tot;i++){
        r=i;
        int id=r;
        for(int j=r+1;j<tot;j++)if(fabs(A[j][i])>fabs(A[id][i])){
            id=j;
        }
        //if(fabs(A[id][i])<eps)continue;

        for(int j=0;j<tot;j++)swap(A[r][j],A[id][j]);
        for(int j=0;j<n;j++)swap(B[r][j],B[id][j]);

        for(int j=r+1;j<tot;j++)if(fabs(A[j][i])>eps){
            double col=A[j][i]/A[r][i];
            for(int k=i;k<tot;k++){
                A[j][k]-=A[r][k]*col;
            }
            for(int k=0;k<n;k++){
                B[j][k]-=B[r][k]*col;
            }
        }
        r++;
    }

    for(int i=0;i<n;i++){
        for(int j=tot-1;j>=0;j--){
            X[j][i]=B[j][i];
            for(int k=tot-1;k>j;k--){
                X[j][i]-=X[k][i]*A[j][k];
            }
            X[j][i]/=A[j][j];
        }
    }
}

int main(){
    #ifdef DouBi
    freopen("in.cpp","r",stdin);
    #endif // DouBi

    while(scanf("%d%d%d%d",&n,&m,&s,&t)!=EOF){
        memset(du,0,sizeof(du));
        memset(g,0,sizeof(g));
        s--;t--;
        for(int i=0;i<m;i++){
            int a,b;scanf("%d%d",&a,&b);a--;b--;
            du[a]++;du[b]++;
            g[a][b]=g[b][a]=1;
        }
        for(int i=0;i<n;i++){
            scanf("%lf",&pp[i]);
        }
        getA();
        getB();
        Gauss();
        for(int i=0;i<n;i++){
            printf("%.10lf%c",X[getid(s,t)][i],i==n-1?'\n':' ');
        }
    }
    return 0;
}