考虑 $dp[i][j]$ 表示结束前 Petya 在 $i$ ，Vesya 在 $j$ 的概率。

那么很明显这东西是不能 dp 转移的。但是我们可以列很多个方程

$$dp[i][j] = p_j\sum_{v \in E(i)} \frac{1-p_v}{deg_v} dp[v][j] + p_i\sum_{u \in E(j)} \frac{1-p_u}{deg_u} + \sum_{v \in E(i)}\sum_{u \in E(j)} \frac{(1-p_v)(1-p_u)}{deg_vdeg_u}dp[v][u] + p_i p_j dp[i][j]$$

这个方程就是考虑，只有一边是从某个点走过来的，两边都是从某个点走过来的，以及都停留在原地的情况。

然后发现这个方程大小是 $400$ ，直接高消即可。

注意，不能从一个终止状态转移到其他状态，所以前面式子里面如果 $u=v$ 之类的，那么系数应当是 $0$ 。

考虑 $dp[a][b]$ ，可以理解为一开始就在 $a,b$ 的概率加上之后到达 $a,b$ 的概率啦。。

#include "iostream"
#include "algorithm"
#include "cstring"
#include "cstdio"
#include "cmath"
#include "vector"
#include "map"
#include "set"
#include "queue"
using namespace std;
#define MAXN 406
//#define int long long
#define rep(i, a, b) for (int i = (a), i##end = (b); i <= i##end; ++i)
#define per(i, a, b) for (int i = (a), i##end = (b); i >= i##end; --i)
#define pii pair<int,int>
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define vi vector<int>
#define all(x) (x).begin() , (x).end()
#define mem( a ) memset( a , 0 , sizeof a )
typedef long long ll;
#define eps 1e-7
int n , m , a , b;
int G[MAXN][MAXN] , deg[MAXN] , id[MAXN][MAXN] , tot = 0;
double p[MAXN];
int sgn( double a ) {
    return fabs( a ) < eps ? 0 : a > 0 ? 1 : -1;
}
double A[MAXN][MAXN];
void solve() {
    cin >> n >> m >> a >> b;
    int u , v;
    rep( i , 1 , m ) scanf("%d%d",&u,&v) , G[u][v] = G[v][u] = 1 , ++ deg[v] , ++ deg[u];
    rep( i , 1 , n ) scanf("%lf",p + i);
    rep( i , 1 , n ) rep( j , 1 , n ) id[i][j] = ++ tot;
    int dx;
    rep( i , 1 , n ) {
        rep( j , 1 , n ) {
            dx = id[i][j];
            rep( v , 1 , n ) if( G[i][v] && j != v ) {
                A[dx][id[v][j]] += p[j] * ( p[v] - 1 ) / deg[v];
            }
            rep( u , 1 , n ) if( G[j][u] && i != u ) {
                A[dx][id[i][u]] += p[i] * ( p[u] - 1 ) / deg[u];
            }
            rep( v , 1 , n ) if( G[i][v] )
                rep( u , 1 , n ) if( G[u][j] && u != v ) {
                    A[dx][id[v][u]] -= ( 1 - p[v] ) * ( 1 - p[u] ) / deg[v] / deg[u];
                }
            A[dx][dx] = 1 - ( ( i != j ) ? p[i] * p[j] : 0 );
        }
    }
    A[id[a][b]][tot + 1] = 1;
//    rep( i , 1 , tot ) { rep( j , 1 , tot + 1 ) cout << A[i][j] << ' '; puts(""); }
    for( int i = 1 ; i <= tot ; ++ i ) {
        int mxpos = i;
        for( int j = i + 1 ; j <= tot ; ++ j ) if( fabs( A[j][i] ) - fabs( A[mxpos][i] ) > eps )
            mxpos = j;
        swap( A[i] , A[mxpos] );
        double c = A[i][i];
        A[i][i] = 1;
        rep( j , i + 1 , tot + 1 ) A[i][j] /= c;
        for( int j = 1 ; j <= tot ; ++ j ) if( j != i ) {
            double cc = A[j][i]; A[j][i] = 0;
            for( int k = i + 1 ; k <= tot + 1 ; ++ k ) A[j][k] -= cc * A[i][k];
        }
    }
    for( int i = 1 ; i <= n ; ++ i ) printf("%lf ",A[id[i][i]][tot + 1]);
}

signed main() {
//    freopen("input","r",stdin);
//    freopen("fuckout","w",stdout);
//    int T;cin >> T;while( T-- ) solve();
    solve();
}