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>
#pragma GCC optimize("Ofast")
#pragma GCC target( \
    "sse,sse2,sse3,ssse3,sse4,sse4.1,sse4.2,popcnt,abm,mmx,avx,tune=native")
using namespace std;
const int MOD = 1000000007;
const int UNDEF = -1;
const int INF = 1 << 30;
template <typename T>
inline bool chkmax(T &aa, T bb) {
  return aa < bb ? aa = bb, true : false;
}
template <typename T>
inline bool chkmin(T &aa, T bb) {
  return aa > bb ? aa = bb, true : false;
}
int rint();
char rch();
long long rlong();
class MatrixInverse {
 public:
  bool zero(double x) { return fabs(x) < 1e-9; }
  bool inv(vector<vector<double> > &M) {
    int N = M.size();
    for (int i = 0; i < N; i++) M[i].resize(2 * N);
    for (int i = 0; i < N; i++) {
      M[i][i + N] = 1;
    }
    for (int r = 0; r < N; r++) {
      int k = r;
      while (zero(M[k][r])) {
        k++;
        if (k == N) return false;
      }
      swap(M[r], M[k]);
      double lv = M[r][r];
      for (int j = 0; j < 2 * N; j++) M[r][j] /= lv;
      for (int i = 0; i < N; i++) {
        if (i != r) {
          lv = M[i][r];
          for (int j = 0; j < 2 * N; j++) M[i][j] -= lv * M[r][j];
        }
      }
    }
    return true;
  }
  vector<double> solution(vector<vector<double> > &M, vector<double> &inp) {
    int n = M.size();
    vector<double> sol(n);
    for (int r = 0; r < n; r++) {
      for (int c = 0; c < n; c++) {
        sol[r] += M[r][c + n] * inp[c];
      }
    }
    return sol;
  }
};
const int mn = 22;
vector<int> g[mn];
double vp[mn];
int gn;
vector<vector<double> > M;
void ad(int sx, int sy, int tx, int ty, double val) {
  int n = gn;
  M[sx * n + sy][tx * n + ty] += val;
}
double vans[mn];
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  int n, m, a, b;
  cin >> n >> m >> a >> b;
  --a;
  --b;
  gn = n;
  if (n == 1) {
    printf("1\n");
    exit(0);
  }
  for (int i = 0; i < m; i++) {
    int x, y;
    cin >> x >> y;
    --x;
    --y;
    g[x].push_back(y);
    g[y].push_back(x);
  }
  for (int i = 0; i < n; i++) cin >> vp[i];
  M;
  M.resize(n * n);
  for (int i = 0; i < n * n; i++) M[i].resize(n * n);
  for (int sx = 0; sx < n; sx++) {
    for (int sy = 0; sy < n; sy++) {
      if (sx == sy) {
        ad(sx, sy, sx, sy, 1);
        continue;
      }
      ad(sx, sy, sx, sy, vp[sx] * vp[sy] - 1);
      {
        double prob = (1 - vp[sx]) * (vp[sy]) / (g[sx].size());
        for (auto &tx : g[sx]) {
          ad(sx, sy, tx, sy, prob);
        }
      }
      {
        double prob = (vp[sx]) * (1 - vp[sy]) / (g[sy].size());
        for (auto &ty : g[sy]) {
          ad(sx, sy, sx, ty, prob);
        }
      }
      {
        double prob =
            (1 - vp[sx]) * (1 - vp[sy]) / (g[sx].size() * g[sy].size());
        for (auto &tx : g[sx]) {
          for (auto &ty : g[sy]) {
            ad(sx, sy, tx, ty, prob);
          }
        }
      }
    }
  }
  MatrixInverse matinv;
  matinv.inv(M);
  for (int r = 0; r < n; r++) {
    vector<double> v(n * n);
    v[r * n + r] = 1;
    vector<double> ans = matinv.solution(M, v);
    printf("%.10f ", ans[a * n + b]);
  }
  printf("\n");
}
static char stdinBuffer[1024];
static char *stdinDataEnd = stdinBuffer + sizeof(stdinBuffer);
static const char *stdinPos = stdinDataEnd;
void readAhead(size_t amount) {
  size_t remaining = stdinDataEnd - stdinPos;
  if (remaining < amount) {
    memmove(stdinBuffer, stdinPos, remaining);
    size_t sz = fread(stdinBuffer + remaining, 1,
                      sizeof(stdinBuffer) - remaining, stdin);
    stdinPos = stdinBuffer;
    stdinDataEnd = stdinBuffer + remaining + sz;
    if (stdinDataEnd != stdinBuffer + sizeof(stdinBuffer)) *stdinDataEnd = 0;
  }
}
int rint() {
  readAhead(16);
  int x = 0;
  bool neg = false;
  while (*stdinPos == ' ' || *stdinPos == '\n') ++stdinPos;
  if (*stdinPos == '-') {
    ++stdinPos;
    neg = true;
  }
  while (*stdinPos >= '0' && *stdinPos <= '9') {
    x *= 10;
    x += *stdinPos - '0';
    ++stdinPos;
  }
  return neg ? -x : x;
}
char rch() {
  readAhead(16);
  while (*stdinPos == ' ' || *stdinPos == '\n') ++stdinPos;
  char ans = *stdinPos;
  ++stdinPos;
  return ans;
}
long long rlong() {
  readAhead(32);
  long long x = 0;
  bool neg = false;
  while (*stdinPos == ' ' || *stdinPos == '\n') ++stdinPos;
  if (*stdinPos == '-') {
    ++stdinPos;
    neg = true;
  }
  while (*stdinPos >= '0' && *stdinPos <= '9') {
    x *= 10;
    x += *stdinPos - '0';
    ++stdinPos;
  }
  return neg ? -x : x;
}