#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")


struct Heap {
  Heap *l, *r;
  // Int key;
  Int p, q;
  // Heap(Int key_ = 0) : l(nullptr), r(nullptr), key(key_) {}
  Heap() : l(nullptr), r(nullptr), p(0), q(0) {}
};
bool operator<(const Heap &a, const Heap &b) {
  return ((__int128)a.p * b.q < (__int128)b.p * a.q);
}
// max heap
Heap *meld(Heap *a, Heap *b) {
  if (!a) return b;
  if (!b) return a;
  // if (a->key < b->key) swap(a, b);
  if (*a < *b) swap(a, b);
  a->r = meld(a->r, b);
  swap(a->l, a->r);
  return a;
}
void pop(Heap *&a) {
  a = meld(a->l, a->r);
}


int N;
vector<int> P;
vector<Int> W;

vector<vector<int>> graph;

namespace brute {
double dfs(int u, int p, double thr) {
  double ret = W[u] - thr;
  for (const int v : graph[u]) if (p != v) {
    const double res = dfs(v, u, thr);
    ret += max(res, 0.0);
  }
  return ret;
}
vector<double> run() {
  const Int maxW = *max_element(W.begin(), W.end());
  vector<double> ans(N);
  for (int u = 0; u < N; ++u) {
    double lo = W[u], hi = maxW;
    for (int iter = 0; iter < 70; ++iter) {
      const double mid = (lo + hi) / 2.0;
      const auto res = dfs(u, -1, mid);
      ((res >= 0.0) ? lo : hi) = mid;
    }
    ans[u] = lo;
  }
  return ans;
}
}  // brute

int main() {
  for (; ~scanf("%d", &N); ) {
    P.assign(N, -1);
    for (int u = 1; u < N; ++u) {
      scanf("%d", &P[u]);
      --P[u];
    }
    W.resize(N);
    for (int u = 0; u < N; ++u) {
      scanf("%lld", &W[u]);
    }
    
    graph.assign(N, {});
    for (int u = 1; u < N; ++u) {
      graph[P[u]].push_back(u);
      graph[u].push_back(P[u]);
    }
    
    const auto ans = brute::run();
    for (int u = 0; u < N; ++u) {
      printf("%.10f\n", ans[u]);
    }
break;
  }
  return 0;
}