Problem: A maze is represented by a tree (an undirected graph, where exactly one way exists between each pair of vertices). In the maze the entrance vertex and the exit vertex are chosen with some probability. The exit from the maze is sought by Deep First Search. If there are several possible ways to move, the move is chosen equiprobably. Consider the following pseudo-code:
    
    
      
    DFS(x)  
        if x == exit vertex then  
            finish search  
        flag[x] <- TRUE  
        random shuffle the vertices' order in V(x) // here all permutations have equal probability to be chosen  
        for i <- 1 to length[V] do  
            if flag[V[i]] = FALSE then  
                count++;  
                DFS(y);  
        count++;  
    

V(x) is the list vertices adjacent to x. The flag array is initially filled as FALSE. DFS initially starts with a parameter of an entrance vertex. When the search is finished, variable count will contain the number of moves.

Your task is to count the mathematical expectation of the number of moves one has to do to exit the maze.

Input

The first line determines the number of vertices in the graph n (1 ≤ n ≤ 105). The next n - 1 lines contain pairs of integers ai and bi, which show the existence of an edge between ai and bi vertices (1 ≤ ai, bi ≤ n). It is guaranteed that the given graph is a tree.

Next n lines contain pairs of non-negative numbers xi and yi, which represent the probability of choosing the i-th vertex as an entrance and exit correspondingly. The probabilities to choose vertex i as an entrance and an exit equal <image> and <image> correspondingly. The sum of all xi and the sum of all yi are positive and do not exceed 106.

Output

Print the expectation of the number of moves. The absolute or relative error should not exceed 10 - 9.

Examples

Input

2
1 2
0 1
1 0


Output

1.00000000000000000000


Input

3
1 2
1 3
1 0
0 2
0 3


Output

2.00000000000000000000


Input

7
1 2
1 3
2 4
2 5
3 6
3 7
1 1
1 1
1 1
1 1
1 1
1 1
1 1


Output

4.04081632653

Note

In the first sample the entrance vertex is always 1 and the exit vertex is always 2.

In the second sample the entrance vertex is always 1 and the exit vertex with the probability of 2/5 will be 2 of with the probability if 3/5 will be 3. The mathematical expectations for the exit vertices 2 and 3 will be equal (symmetrical cases). During the first move one can go to the exit vertex with the probability of 0.5 or to go to a vertex that's not the exit vertex with the probability of 0.5. In the first case the number of moves equals 1, in the second one it equals 3. The total mathematical expectation is counted as 2 / 5 × (1 × 0.5 + 3 × 0.5) + 3 / 5 × (1 × 0.5 + 3 × 0.5)

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int oo = 0x3f3f3f3f;
const double eps = 1e-9;
vector<int> adj[100100];
int sp[100100], fp[100100];
int ed[100100][2], N;
double dp[100100][2];
int pro[100100][2], sfp, ssp;
int useless[100100][2];
void calc1(int ii, int nr) {
  int &pr = pro[ii][nr];
  int &st = useless[ii][nr];
  if (pr != -1) return;
  st = 0;
  int n = ed[ii][nr];
  pr = fp[n];
  for (int i = (0); i < (int((adj[n]).size())); i++)
    if (adj[n][i] != ii) {
      int k = 0;
      if (ed[adj[n][i]][0] == n) k = 1;
      calc1(adj[n][i], k);
      pr += pro[adj[n][i]][k];
      st += useless[adj[n][i]][k] + 2;
    }
  for (int i = (0); i < (int((adj[n]).size())); i++)
    if (adj[n][i] != ii) {
      int k = 0;
      if (ed[adj[n][i]][0] == n) k = 1;
      pro[adj[n][i]][k ^ 1] = sfp - pro[adj[n][i]][k];
      useless[adj[n][i]][k ^ 1] = 2 * (N - 1) - useless[adj[n][i]][k] - 2;
    }
}
double help[10][100100];
int helpid[100100][2];
void calc(int ii, int nr) {
  int &pr = pro[ii][nr];
  int &st = useless[ii][nr];
  double &res = dp[ii][nr];
  if (res > -0.5) return;
  res = 0;
  int n = ed[ii][nr];
  if (pr != fp[n]) {
    if (helpid[n][0] == -1) {
      help[2][n] = 0.0;
      for (int i = (0); i < (int((adj[n]).size())); i++)
        if (adj[n][i] != ii) {
          int k = 0;
          if (ed[adj[n][i]][0] == n) k = 1;
          calc(adj[n][i], k);
          res += (pro[adj[n][i]][k]) / ((double)pr) * (dp[adj[n][i]][k] + 1);
          help[2][n] += (pro[adj[n][i]][k]) * (dp[adj[n][i]][k] + 1);
        }
      helpid[n][0] = ii;
      helpid[n][1] = nr;
    } else {
      if (helpid[n][0] != -2) {
        calc(helpid[n][0], helpid[n][1] ^ 1);
        help[2][n] += (pro[helpid[n][0]][helpid[n][1] ^ 1]) *
                      (dp[helpid[n][0]][helpid[n][1] ^ 1] + 1);
        helpid[n][0] = -2;
      }
      res += (help[2][n] - pro[ii][nr ^ 1] * (dp[ii][nr ^ 1] + 1)) /
             ((double)(pr));
    }
    double relevantUseless = help[0][n] - (useless[ii][nr ^ 1] + 2) * 0.5;
    res += relevantUseless * (pr - fp[n]) / ((double)pr);
    double relevantUseless2 =
        help[1][n] + (useless[ii][nr ^ 1] + 2) * 0.5 * (double)pro[ii][nr ^ 1];
    res += relevantUseless2 / ((double)pr);
  }
}
int main() {
  scanf("%d", &N);
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) scanf("%d", &(ed[i][j]));
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) --ed[i][j];
  for (int i = (0); i < (N); i++) scanf("%d%d", sp + i, fp + i);
  for (int i = (0); i < (N); i++) adj[i].clear();
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) adj[ed[i][j]].push_back(i);
  ssp = 0;
  for (int i = (0); i < (N); i++) ssp += sp[i];
  sfp = 0;
  for (int i = (0); i < (N); i++) sfp += fp[i];
  memset(pro, -1, sizeof(pro));
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) calc1(i, j);
  for (int i = (0); i < (N); i++) {
    help[0][i] = 0.0;
    help[1][i] = 0.0;
    for (int j = (0); j < (int((adj[i]).size())); j++) {
      int k = 0;
      if (ed[adj[i][j]][0] == i) k = 1;
      help[0][i] += 0.5 * (useless[adj[i][j]][k] + 2);
      help[1][i] -=
          0.5 * (useless[adj[i][j]][k] + 2) * (double)(pro[adj[i][j]][k]);
    }
  }
  memset(helpid, -1, sizeof(helpid));
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) dp[i][j] = -1;
  for (int i = (0); i < (N - 1); i++)
    for (int j = (0); j < (2); j++) calc(i, j);
  double res = 0.0;
  for (int i = (0); i < (N); i++) {
    double cur = 0.0;
    for (int j = (0); j < (int((adj[i]).size())); j++) {
      int k = 0;
      if (ed[adj[i][j]][0] == i) k = 1;
      cur += (pro[adj[i][j]][k] / (double)(sfp)) * (dp[adj[i][j]][k] + 1);
      cur += (sfp - pro[adj[i][j]][k] - fp[i]) / (double)(sfp)*0.5 *
             (useless[adj[i][j]][k] + 2);
    }
    res += cur * sp[i];
  }
  res /= ssp;
  printf("%.12lf\n", res);
  return 0;
}