#include <bits/stdc++.h>
using namespace std;
using i64 = long long;
using ld = long double;
constexpr int maxn = 500000 + 10;
constexpr ld eps = 1e-10;
struct node
{
    int t, v, d, l, r;
    ld x, s, k;
} t[maxn << 2];
int cnt;
void push(int u)
{
    if (!u) return;
    t[u].s = t[t[u].l].s + t[t[u].r].s + t[u].x;
    t[u].t = t[t[u].l].t + t[u].v + t[t[u].r].t;
}
void pop(int o, int &x, int &y)
{
    if (t[o].l) pop(t[o].l, x, t[y = o].l);
    else x = o, y = t[o].r, t[x].l = t[x].r = 0;
    push(o);
}
void spl(int o, ld p, int &x, int &y)
{
    if (!o) return x = y = 0, void();
    if (t[o].k >= p) spl(t[o].r, p, t[x = o].r, y);
    else spl(t[o].l, p, x, t[y = o].l);
    push(o);
}
int mer(int u, int v)
{
    if (!u || !v) return u | v;
    if (t[u].d < t[v].d) swap(u, v);
    int x, y;
    spl(v, t[u].k, x, y);
    t[u].l = mer(t[u].l, x);
    t[u].r = mer(t[u].r, y);
    push(u);
    return u;
}
vector<int> g[maxn];
i64 w[maxn];
mt19937 rng(19491001);
ld ans[maxn];
int dfs(int u)
{
    int x = 0;
    for (int v : g[u])
    {
        int y = dfs(v);
        x = mer(x, y);
    }
    ld k = w[u];
    int v = 1;
    for (; x;)
    {
        int a, b;
        pop(x, a, b);
        if (t[a].k < k / v - eps)
        {
            x = mer(a, b);
            break;
        }
        x = b;
        k += t[a].x;
        v += t[a].v;
    }
    k /= v;
    int p = ++cnt;
    t[p].x = t[p].s = k * v;
    t[p].k = k;
    t[p].t = t[p].v = v;
    t[p].d = rng();
    x = mer(p, x);
    ans[u] = k;
    return x;
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    int n;
    cin >> n;
    for (int i=2;i<=n;++i)
    {
        int x;
        cin >> x;
        g[x].emplace_back(i);
    }
    for (int i=1;i<=n;++i) cin >> w[i];
    dfs(1);
    for (int i=1;i<=n;++i) cout << fixed << setprecision(15) << ans[i] << '\n';
    return 0;
}