Problem: Limak is a little grizzly bear. He will once attack Deerland but now he can only destroy trees in role-playing games. Limak starts with a tree with one vertex. The only vertex has index 1 and is a root of the tree.

Sometimes, a game chooses a subtree and allows Limak to attack it. When a subtree is attacked then each of its edges is destroyed with probability <image>, independently of other edges. Then, Limak gets the penalty — an integer equal to the height of the subtree after the attack. The height is defined as the maximum number of edges on the path between the root of the subtree and any vertex in the subtree.

You must handle queries of two types.

  * 1 v denotes a query of the first type. A new vertex appears and its parent is v. A new vertex has the next available index (so, new vertices will be numbered 2, 3, ...). 
  * 2 v denotes a query of the second type. For a moment let's assume that the game allows Limak to attack a subtree rooted in v. Then, what would be the expected value of the penalty Limak gets after the attack? 



In a query of the second type, Limak doesn't actually attack the subtree and thus the query doesn't affect next queries.

Input

The first line of the input contains one integer q (1 ≤ q ≤ 500 000) — the number of queries.

Then, q lines follow. The i-th of them contains two integers typei and vi (1 ≤ typei ≤ 2). If typei = 1 then vi denotes a parent of a new vertex, while if typei = 2 then you should print the answer for a subtree rooted in vi.

It's guaranteed that there will be at least 1 query of the second type, that is, the output won't be empty.

It's guaranteed that just before the i-th query a vertex vi already exists.

Output

For each query of the second type print one real number —the expected value of the penalty if Limak attacks the given subtree. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 6.

Namely: let's assume that your answer is a, and the answer of the jury is b. The checker program will consider your answer correct if <image>.

Examples

Input

7
1 1
1 1
2 1
1 2
1 3
2 2
2 1


Output

0.7500000000
0.5000000000
1.1875000000


Input

8
2 1
1 1
1 2
1 3
1 4
2 1
1 4
2 1


Output

0.0000000000
0.9375000000
0.9687500000

Note

Below, you can see the drawing for the first sample. Red circles denote queries of the second type.

<image>

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
#pragma GCC target( \
    "sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
char buf[1 << 21], *p1 = buf, *p2 = buf, buffer[1 << 21];
int P1 = -1;
const int P2 = (1 << 21) - 1;
inline void flush() { fwrite(buffer, 1, P1 + 1, stdout), P1 = -1; }
inline void putc(const char &x) {
  if (P1 == P2) flush();
  buffer[++P1] = x;
}
using namespace std;
template <typename T>
inline void read(T &x) {
  x = 0;
  char c =
      (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 21, stdin), p1 == p2)
           ? EOF
           : *p1++);
  for (; c < '0' || c > '9';
       c = (p1 == p2 &&
                    (p2 = (p1 = buf) + fread(buf, 1, 1 << 21, stdin), p1 == p2)
                ? EOF
                : *p1++))
    ;
  for (; c >= '0' && c <= '9';
       c = (p1 == p2 &&
                    (p2 = (p1 = buf) + fread(buf, 1, 1 << 21, stdin), p1 == p2)
                ? EOF
                : *p1++))
    x = (x << 1) + (x << 3) + (c ^ 48);
}
template <typename T>
inline void write(T x) {
  int y = 10, len = 1;
  while (y <= x) y *= 10, len++;
  while (len--) y /= 10, putc(x / y + 48), x %= y;
}
int q, tot, fa[500005];
double F[500005][65];
int main() {
  read(q), tot = 1;
  for (int i = 1; i <= 63; i++) F[1][i] = 1;
  for (int op, x; q--;) {
    read(op), read(x);
    if (op == 2) {
      double res = 62;
      for (int i = 1; i < 63; i++) res -= F[x][i];
      printf("%.10lf\n", res);
      continue;
    }
    vector<int> v;
    fa[++tot] = x, x = tot;
    for (int i = 1; i <= 63; i++) F[x][i] = 1;
    for (int i = 1; i <= 64 && x; i++) v.push_back(x), x = fa[x];
    for (int i = v.size() - 2; i; i--)
      F[v[i + 1]][i + 1] /= (F[v[i]][i] + 1) / 2;
    for (int i = 0; i + 1 < (int)v.size(); i++)
      F[v[i + 1]][i + 1] *= (F[v[i]][i] + 1) / 2;
  }
  return 0;
}