Problem: One day as Petya and his friend Vasya were having one of their numerous trips, they decided to visit a museum castle. The museum has a specific shape: it consists of n rooms connected with m corridors so that one can access any room from any other one.

After the two friends had a little walk around the museum, they decided to split and watch the pieces of art each of them found interesting. They agreed to meet in one of the rooms at six p.m. However, they forgot one quite essential thing: they didn't specify the place to meet and when the time came, they started to rush about the museum looking for each other (they couldn't call each other as roaming made a call's cost skyrocket).

Yet, even despite the whole rush, they couldn't get enough of the pieces of art, that's why each of them has the following strategy: each minute he make a decision where to go — with probability pi he doesn't move to any other place during this minute (i.e. he stays in the room). With probability 1 - pi he equiprobably choose one of the adjacent rooms and went there along the corridor. Here i is the ordinal number of the current room. Building was expensive in ancient times, that's why each corridor connected two different rooms, and any two rooms had no more than one corridor between them. 

The boys act simultaneously. As the corridors are dark, it is impossible to meet there; however, one can walk along the corridors in both directions (besides, the two boys can be going through the same corridor simultaneously without meeting). The boys act like that until they meet each other. More formally, the two friends meet when at some moment of time both of them decided to appear in the same room.

For each room find the probability that the boys will meet there considering that at 6 p.m. they are positioned in rooms a and b correspondingly.

Input

The first line contains four integers: n (1 ≤ n ≤ 22), representing the numbers of rooms; m <image>, representing the number of corridors; a, b (1 ≤ a, b ≤ n), representing the numbers of Petya's and Vasya's starting rooms correspondingly.

Next m lines contain pairs of numbers — the numbers of rooms connected by a corridor. Next n lines contain probabilities pi (0.01 ≤ pi ≤ 0.99) with the accuracy of up to four digits after the decimal point — the probability to stay in room i.

It is guaranteed that every room can be reached from every other room by corridors.

Output

In the only line print n space-separated numbers, the i-th number should represent the probability that the friends meet in the i-th room with absolute or relative error of no more than 10 - 6.

Examples

Input

2 1 1 2
1 2
0.5
0.5


Output

0.5000000000 0.5000000000 

Input

4 4 1 2
1 2
2 3
3 4
4 1
0.5
0.5
0.5
0.5


Output

0.3333333333 0.3333333333 0.1666666667 0.1666666667 

Note

In the first sample the museum is symmetric. That means the probabilities to meet in rooms 1 and 2 are equal. And their sum equals to one. So, each probability equals 0.5.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
namespace IO {
template <class T>
inline void readin(T& x) {
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  if (f) x = -x;
}
template <class T>
inline T readin(T x) {
  x = 0;
  char c;
  bool f = 0;
  while ((c = getchar()) < '0' || '9' < c) f |= (c == '-');
  for (x = (c ^ 48); '0' <= (c = getchar()) && c <= '9';
       x = (x << 1) + (x << 3) + (c ^ 48))
    ;
  return f ? -x : x;
}
template <class T>
void fwrit(const T x) {
  if (x < 0) return (void)(putchar('-'), fwrit(-x));
  if (x > 9) fwrit(x / 10);
  putchar(x % 10 ^ 48);
}
template <class T>
inline T Max(const T x, const T y) {
  return x < y ? y : x;
}
template <class T>
inline T Min(const T x, const T y) {
  return x < y ? x : y;
}
template <class T>
inline T fab(const T x) {
  return x > 0 ? x : -x;
}
inline int gcd(const int a, const int b) { return b ? gcd(b, a % b) : a; }
inline void getInv(int inv[], const int lim, const int MOD) {
  inv[0] = inv[1] = 1;
  for (int i = 2; i <= lim; ++i)
    inv[i] = 1ll * inv[MOD % i] * (MOD - MOD / i) % MOD;
}
inline long long mulMod(const long long a, const long long b,
                        const long long mod) {
  return ((a * b - (long long)((long double)a / mod * b + 1e-8) * mod) % mod +
          mod) %
         mod;
}
}  // namespace IO
using namespace IO;
const int maxn = 22;
const int maxm = maxn * (maxn + 1) / 2;
const double eps = 1e-9;
struct edge {
  int to, nxt;
} e[maxm * 2 + 5];
int tail[maxn + 5], ecnt, d[maxn + 5];
inline void add_edge(const int u, const int v) {
  e[++ecnt] = edge{v, tail[u]};
  tail[u] = ecnt;
  e[++ecnt] = edge{u, tail[v]};
  tail[v] = ecnt;
  ++d[u], ++d[v];
}
int n, m, a, b, N;
inline int getid(const int x, const int y) { return (x - 1) * n + y; }
double p[maxn + 5];
double f[maxn * maxn + 5][maxn * maxn + 5], x[maxn * maxn + 5];
inline void Init() {
  cin >> n >> m >> a >> b;
  N = n * n;
  int u, v;
  for (signed i = (1), i_end_ = (m); i <= i_end_; ++i) {
    cin >> u >> v;
    add_edge(u, v);
  }
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i) cin >> p[i];
}
inline void Build_f() {
  f[getid(a, b)][N + 1] = 1;
  int u;
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i)
    for (signed j = (1), j_end_ = (n); j <= j_end_; ++j) {
      u = getid(i, j);
      f[u][u] = 1;
      if (i != j) f[u][u] -= p[i] * p[j];
      for (int t = tail[i], x; t; t = e[t].nxt)
        if ((x = e[t].to) != j) f[u][getid(x, j)] -= (1 - p[x]) / d[x] * p[j];
      for (int t = tail[j], x; t; t = e[t].nxt)
        if (i != (x = e[t].to)) f[u][getid(i, x)] -= p[i] * (1 - p[x]) / d[x];
      for (int t1 = tail[i], x; t1; t1 = e[t1].nxt)
        for (int t2 = tail[j], y; t2; t2 = e[t2].nxt)
          if ((x = e[t1].to) != (y = e[t2].to))
            f[u][getid(x, y)] -= (1 - p[x]) / d[x] * (1 - p[y]) / d[y];
    }
}
inline void Gauss() {
  for (signed i = (1), i_end_ = (N); i <= i_end_; ++i) {
    int pos = i;
    for (signed j = (i + 1), j_end_ = (N); j <= j_end_; ++j)
      if (fab(f[j][i]) > fab(f[pos][i])) pos = j;
    if (pos != i) swap(f[pos], f[i]);
    if (fab(f[i][i]) < eps) continue;
    double tmp = f[i][i];
    for (signed j = (i), j_end_ = (N + 1); j <= j_end_; ++j) f[i][j] /= tmp;
    for (signed j = (i + 1), j_end_ = (N); j <= j_end_; ++j) {
      tmp = f[j][i];
      for (signed k = (i), k_end_ = (N + 1); k <= k_end_; ++k)
        f[j][k] -= tmp * f[i][k];
    }
  }
  for (signed i = (N), i_end_ = (1); i >= i_end_; --i) {
    x[i] = f[i][N + 1];
    for (signed j = (i - 1), j_end_ = (1); j >= j_end_; --j)
      f[j][N + 1] -= f[j][i] * x[i];
  }
}
signed main() {
  Init();
  Build_f();
  Gauss();
  for (signed i = (1), i_end_ = (n); i <= i_end_; ++i)
    cout << fixed << setprecision(9) << x[getid(i, i)] << ' ';
  return 0;
}