#include <bits/stdc++.h>
using namespace std;
const long long int N = 1e3 + 2;
long long int mod = 1e9 + 7;
long long int n;
long long int dp[N][N];
vector<long long int> v(N);
long long int fun(long long int i, long long int j) {
  if (i == j) {
    return 1;
  }
  long long int &ans = dp[i][j];
  if (ans != -1) {
    return ans;
  }
  if (i < j) {
    if (i + 1 == j) {
      ans = fun(i, v[i]) + 1;
    } else {
      ans = fun(i, j - 1) + fun(j - 1, v[j - 1]) + 1;
    }
  } else {
    ans = 1 + fun(j, i);
  }
  ans %= mod;
  return ans;
}
int32_t main() {
  scanf("%lld", &n);
  for (long long int i = 1; i <= n; ++i) {
    scanf("%lld", &v[i]);
  }
  memset(dp, -1, sizeof dp);
  long long int ans = fun(1, n + 1);
  printf("%lld", ans);
  return 0;
}
