Problem p03754 (C++)

Buggy submission (status: Runtime Error):
// #pragma GCC target("avx")  // CPU 処理並列化
// #pragma GCC optimize("O3")  // CPU 処理並列化
// #pragma GCC optimize("unroll-loops")  // 条件処理の呼び出しを減らす
#include <algorithm>
#include <bitset>
#include <cassert>
#include <complex>
#include <deque>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <time.h>
#include <unordered_set>
#include <vector>
using namespace std;
const long long mod = 1000000007;
const long long inf = mod * mod;
const long long d2 = (mod + 1) / 2;
const double EPS = 1e-10;
const double INF = 1e+10;
const double PI = acos(-1.0);
const int C_SIZE = 3100000;
const int UF_SIZE = 210000;
namespace {
long long fact[C_SIZE];
long long finv[C_SIZE];
long long inv[C_SIZE];
long long Comb(int a, int b) {
  if (a < b || b < 0)
    return 0;
  return fact[a] * finv[b] % mod * finv[a - b] % mod;
}
void init_C(int n) {
  fact[0] = finv[0] = inv[1] = 1;
  for (int i = 2; i < n; i++) {
    inv[i] = (mod - (mod / i) * inv[mod % i] % mod) % mod;
  }
  for (int i = 1; i < n; i++) {
    fact[i] = fact[i - 1] * i % mod;
    finv[i] = finv[i - 1] * inv[i] % mod;
  }
}
long long pw(long long a, long long b) {
  if (a < 0LL)
    return 0;
  if (b < 0LL)
    return 0;
  long long ret = 1;
  while (b) {
    if (b % 2)
      ret = ret * a % mod;
    a = a * a % mod;
    b /= 2;
  }
  return ret;
}
long long pw_mod(long long a, long long b, long long M) {
  if (a < 0LL)
    return 0;
  if (b < 0LL)
    return 0;
  long long ret = 1;
  while (b) {
    if (b % 2)
      ret = ret * a % M;
    a = a * a % M;
    b /= 2;
  }
  return ret;
}
int pw_mod_int(int a, int b, int M) {
  if (a < 0)
    return 0;
  if (b < 0)
    return 0;
  int ret = 1;
  while (b) {
    if (b % 2)
      ret = (long long)ret * a % M;
    a = (long long)a * a % M;
    b /= 2;
  }
  return ret;
}
int ABS(int a) { return max(a, -a); }
long long ABS(long long a) { return max(a, -a); }
double ABS(double a) { return max(a, -a); }
int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
int UF[UF_SIZE];
void init_UF(int n) {
  for (int i = 0; i < n; i++)
    UF[i] = -1;
}
int FIND(int a) {
  if (UF[a] < 0)
    return a;
  return UF[a] = FIND(UF[a]);
}
void UNION(int a, int b) {
  a = FIND(a);
  b = FIND(b);
  if (a == b)
    return;
  if (UF[a] > UF[b])
    swap(a, b);
  UF[a] += UF[b];
  UF[b] = a;
}
} // namespace
// ここから編集しろ
struct TreeDP {
  typedef double vint;
  static const int maxN = 110000;
  struct edge {
    int to;
    int rev;
    // ここに何かをもたせたいときはその都度編集する
    edge(int a, int b) {
      to = a;
      rev = b;
    }
    edge() { to = rev = 0; }
  };
  vector<edge> g[maxN];
  vector<int> lp;
  vector<int> rp;
  vector<vint> L[maxN];
  vector<vint> R[maxN];
  int N;
  vint dfs(int a, int b) {
    L[a][0] = R[a][0] = 0;
    while (lp[a] < g[a].size() && b != lp[a]) {
      L[a][lp[a] + 1] =
          L[a][lp[a]] + (1 + dfs(g[a][lp[a]].to, g[a][lp[a]].rev));
      lp[a]++;
    }
    while (rp[a] < g[a].size() && b != g[a].size() - 1 - rp[a]) {
      R[a][rp[a] + 1] =
          R[a][rp[a]] + (1 + dfs(g[a][g[a].size() - 1 - rp[a]].to,
                                 g[a][g[a].size() - 1 - rp[a]].rev));
      rp[a]++;
    }
    if (b == -1) {
      if (g[a].size() == 0)
        return 0;
      return L[a][g[a].size()] / g[a].size();
    } else {
      if (g[a].size() == 1)
        return 0;
      return (L[a][b] + R[a][g[a].size() - b - 1]) / (g[a].size() - 1);
    }
  }

  // 何か重みのある辺が追加したいときはここを変える
  // 各辺について2度以上呼ばないようにする
  void add_edge(int from, int to) {
    g[from].push_back(edge(to, g[to].size()));
    g[to].push_back(edge(from, g[from].size() - 1));
  }
  vector<vint> solve(int n) {
    N = n;
    for (int i = 0; i < n; i++) {
      L[i] = vector<vint>(g[i].size() + 1);
      R[i] = vector<vint>(g[i].size() + 1);
    }
    lp = vector<int>(n);
    rp = vector<int>(n);
    vector<vint> ret(n);
    for (int i = 0; i < n; i++) {
      ret[i] = dfs(i, -1);
    }
    return ret;
  }
};

int main() {
  int a;
  scanf("%d", &a);
  TreeDP td;
  for (int i = 0; i < a - 1; i++) {
    int p, q;
    scanf("%d%d", &p, &q);
    p--;
    q--;
    td.add_edge(p, q);
  }
  vector<double> ret = td.solve(a);
  for (int i = 0; i < a; i++)
    printf("%.12f\n", ret[i]);
}


Fixed submission:
// #pragma GCC target("avx")  // CPU 処理並列化
// #pragma GCC optimize("O3")  // CPU 処理並列化
// #pragma GCC optimize("unroll-loops")  // 条件処理の呼び出しを減らす
#include <algorithm>
#include <bitset>
#include <cassert>
#include <complex>
#include <deque>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <time.h>
#include <unordered_set>
#include <vector>
using namespace std;
const long long mod = 1000000007;
const long long inf = mod * mod;
const long long d2 = (mod + 1) / 2;
const double EPS = 1e-10;
const double INF = 1e+10;
const double PI = acos(-1.0);
const int C_SIZE = 3100000;
const int UF_SIZE = 210000;
namespace {
long long fact[C_SIZE];
long long finv[C_SIZE];
long long inv[C_SIZE];
long long Comb(int a, int b) {
  if (a < b || b < 0)
    return 0;
  return fact[a] * finv[b] % mod * finv[a - b] % mod;
}
void init_C(int n) {
  fact[0] = finv[0] = inv[1] = 1;
  for (int i = 2; i < n; i++) {
    inv[i] = (mod - (mod / i) * inv[mod % i] % mod) % mod;
  }
  for (int i = 1; i < n; i++) {
    fact[i] = fact[i - 1] * i % mod;
    finv[i] = finv[i - 1] * inv[i] % mod;
  }
}
long long pw(long long a, long long b) {
  if (a < 0LL)
    return 0;
  if (b < 0LL)
    return 0;
  long long ret = 1;
  while (b) {
    if (b % 2)
      ret = ret * a % mod;
    a = a * a % mod;
    b /= 2;
  }
  return ret;
}
long long pw_mod(long long a, long long b, long long M) {
  if (a < 0LL)
    return 0;
  if (b < 0LL)
    return 0;
  long long ret = 1;
  while (b) {
    if (b % 2)
      ret = ret * a % M;
    a = a * a % M;
    b /= 2;
  }
  return ret;
}
int pw_mod_int(int a, int b, int M) {
  if (a < 0)
    return 0;
  if (b < 0)
    return 0;
  int ret = 1;
  while (b) {
    if (b % 2)
      ret = (long long)ret * a % M;
    a = (long long)a * a % M;
    b /= 2;
  }
  return ret;
}
int ABS(int a) { return max(a, -a); }
long long ABS(long long a) { return max(a, -a); }
double ABS(double a) { return max(a, -a); }
int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
int UF[UF_SIZE];
void init_UF(int n) {
  for (int i = 0; i < n; i++)
    UF[i] = -1;
}
int FIND(int a) {
  if (UF[a] < 0)
    return a;
  return UF[a] = FIND(UF[a]);
}
void UNION(int a, int b) {
  a = FIND(a);
  b = FIND(b);
  if (a == b)
    return;
  if (UF[a] > UF[b])
    swap(a, b);
  UF[a] += UF[b];
  UF[b] = a;
}
} // namespace
// ここから編集しろ
struct TreeDP {
  typedef double vint;
  static const int maxN = 210000;
  struct edge {
    int to;
    int rev;
    // ここに何かをもたせたいときはその都度編集する
    edge(int a, int b) {
      to = a;
      rev = b;
    }
    edge() { to = rev = 0; }
  };
  vector<edge> g[maxN];
  vector<int> lp;
  vector<int> rp;
  vector<vint> L[maxN];
  vector<vint> R[maxN];
  int N;
  vint dfs(int a, int b) {
    L[a][0] = R[a][0] = 0;
    while (lp[a] < g[a].size() && b != lp[a]) {
      L[a][lp[a] + 1] =
          L[a][lp[a]] + (1 + dfs(g[a][lp[a]].to, g[a][lp[a]].rev));
      lp[a]++;
    }
    while (rp[a] < g[a].size() && b != g[a].size() - 1 - rp[a]) {
      R[a][rp[a] + 1] =
          R[a][rp[a]] + (1 + dfs(g[a][g[a].size() - 1 - rp[a]].to,
                                 g[a][g[a].size() - 1 - rp[a]].rev));
      rp[a]++;
    }
    if (b == -1) {
      if (g[a].size() == 0)
        return 0;
      return L[a][g[a].size()] / g[a].size();
    } else {
      if (g[a].size() == 1)
        return 0;
      return (L[a][b] + R[a][g[a].size() - b - 1]) / (g[a].size() - 1);
    }
  }

  // 何か重みのある辺が追加したいときはここを変える
  // 各辺について2度以上呼ばないようにする
  void add_edge(int from, int to) {
    g[from].push_back(edge(to, g[to].size()));
    g[to].push_back(edge(from, g[from].size() - 1));
  }
  vector<vint> solve(int n) {
    N = n;
    for (int i = 0; i < n; i++) {
      L[i] = vector<vint>(g[i].size() + 1);
      R[i] = vector<vint>(g[i].size() + 1);
    }
    lp = vector<int>(n);
    rp = vector<int>(n);
    vector<vint> ret(n);
    for (int i = 0; i < n; i++) {
      ret[i] = dfs(i, -1);
    }
    return ret;
  }
};

int main() {
  int a;
  scanf("%d", &a);
  TreeDP td;
  for (int i = 0; i < a - 1; i++) {
    int p, q;
    scanf("%d%d", &p, &q);
    p--;
    q--;
    td.add_edge(p, q);
  }
  vector<double> ret = td.solve(a);
  for (int i = 0; i < a; i++)
    printf("%.12f\n", ret[i]);
}