#include <bits/stdc++.h>
using namespace std;
const int N = 3005;
const int M = 3e3 + 5;
const int MOD = 1e9 + 7;
const int CM = 998244353;
const int INF = 0x3f3f3f3f;
const double eps = 1e-6;
std::vector<int> g[N];
int dp[N][N];
int pre[N];
int n, D, m;
long long pow(long long a, long long b, long long p) {
  long long ans = 1;
  a %= p;
  while (b) {
    if (b & 1) ans = ans * a % p;
    a = a * a % p;
    b >>= 1;
  }
  return ans;
}
inline long long add(long long x, long long y) {
  return x + y >= MOD ? x + y - MOD : x + y;
}
inline long long dec(long long x, long long y) {
  return x - y < 0 ? x - y + MOD : x - y;
}
inline long long mul(long long x, long long y) { return x * y % MOD; }
long long cal2(long long k) {
  long long ans = 0;
  for (int i = 1; i <= m; i++) {
    long long res = 1;
    for (int j = 1; j <= m; j++) {
      if (i != j) res = mul(res, dec(i, j));
    }
    res = pow(res, MOD - 2, MOD);
    for (int j = 1; j <= m; j++) {
      if (i != j) res = mul(res, dec(k, j));
    }
    res = mul(res, pre[i]);
    ans = add(ans, res);
  }
  return ans;
}
void dfs(int u, int fath) {
  for (int j = 1; j <= m; j++) dp[u][j] = 1;
  for (int i = 0; i < g[u].size(); i++) {
    int v = g[u][i];
    dfs(v, u);
    pre[0] = 0;
    for (int j = 1; j <= m; j++) pre[j] = add(pre[j - 1], dp[v][j]);
    for (int j = 1; j <= m; j++) dp[u][j] = mul(dp[u][j], pre[j]);
  }
}
void solve() {
  scanf("%d%d", &n, &D);
  m = n + 1;
  for (int i = 1; i < n; i++) {
    int p;
    scanf("%d", &p);
    g[p].push_back(i + 1);
  }
  dfs(1, 0);
  pre[0] = 0;
  for (int i = 1; i <= m; i++) pre[i] = add(pre[i - 1], dp[1][i]);
  printf("%lld\n", cal2(D));
}
int main() {
  solve();
  return 0;
}
