Random Function and Tree

You have a rooted tree consisting of n vertices. Let’s number them with integers from 1 to n inclusive. The root of the tree is the vertex 1. For each i > 1 direct parent of the vertex i is p i. We say that vertex i is child for its direct parent p i.

You have initially painted all the vertices with red color. You like to repaint some vertices of the tree. To perform painting you use the function paint that you call with the root of the tree as an argument. Here is the pseudocode of this function:


count = 0 // global integer variable

rnd() { // this function is used in paint code
return 0 or 1 equiprobably
}

paint(s) {
if (count is even) then paint s with white color
else paint s with black color

count = count + 1

if rnd() = 1 then children = [array of vertex s children in ascending order of their numbers]
else children = [array of vertex s children in descending order of their numbers]

for child in children { // iterating over children array
if rnd() = 1 then paint(child) // calling paint recursively
}
}

As a result of this function, some vertices may change their colors to white or black and some of them may remain red.

Your task is to determine the number of distinct possible colorings of the vertices of the tree. We will assume that the coloring is possible if there is a nonzero probability to get this coloring with a single call of paint(1). We assume that the colorings are different if there is a pair of vertices that are painted with different colors in these colorings. Since the required number may be very large, find its remainder of division by 1000000007 (109 + 7).

Input

The first line contains a single integer n (2 ≤ n ≤ 105) — the number of vertexes in the tree.

The second line contains n - 1 integers p 2, p 3, …, p n (1 ≤ p i < i). Number p i is the parent of vertex i.

Output

Print a single integer — the answer to the problem modulo 1000000007 (109 + 7)

Examples

input

4
1 2 1

output

8

input

3
1 1

output

5

Note

All possible coloring patterns of the first sample are given below.

Solution:

#include <cstring>
#include <vector>
#include  <list>
#include <map>
#include <set>
#include <deque>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <memory.h>
#include <cassert>

using namespace std;

const int md = 1000000007;

inline void add(int &a, int b) {
a += b;
if (a >= md) {
a -= md;
}
}

inline int mul(int a, int b) {
return (long long)a * b % md;
}

const int N = 200010;

vector <int> g[N];
int f[N][2];
int aux[7][7], new_aux[7][7];

int main() {
int n;
scanf("%d", &n);
for (int i = 1; i <= n; i++) {
    g[i].clear();
  }
  for (int i = 2; i <= n; i++) {
    int foo;
    scanf("%d", &foo);
    g[foo].push_back(i);
  }
  for (int i = n; i >= 1; i--) {
for (int tot = 0; tot < 2; tot++) {
      for (int all = 0; all < 4; all++) {
        aux[tot][all] = 0;
      }
    }
    aux[0][0] = 1;
    int sz = g[i].size();
    for (int j = 0; j < sz; j++) {
      int u = g[i][j];
      for (int tot = 0; tot < 2; tot++) {
        for (int all = 0; all < 4; all++) {
          new_aux[tot][all] = aux[tot][all];
        }
      }
      for (int tot = 0; tot < 2; tot++) {
        for (int all = 0; all < 4; all++) {
          for (int cur = 0; cur < 2; cur++) {
            add(new_aux[tot ^ cur][all | (1 << cur)], mul(aux[tot][all], f[u][cur]));
          }
        }
      }
      for (int tot = 0; tot < 2; tot++) {
        for (int all = 0; all < 4; all++) {
          aux[tot][all] = new_aux[tot][all];
        }
      }
    }
    f[i][0] = f[i][1] = 0;
    for (int tot = 0; tot < 2; tot++) {
      for (int all = 0; all < 4; all++) {
        add(f[i][tot ^ 1], aux[tot][all]);
        if (all == 3 || (all == 2 && tot == 0)) {
          add(f[i][tot ^ 1], aux[tot][all]);
        }
      }
    }
  }
  printf("%d\n", (f[1][0] + f[1][1]) % md);
  return 0;
}