#include <bits/stdc++.h>
using namespace std;
const int dx[4] = {-1, 0, 1, 0};
const int dy[4] = {0, 1, 0, -1};
long long POW(long long n, long long p) {
  if (p == 0) return 1;
  long long t = POW(n, p / 2);
  if (p & 1) return n * t * t;
  return t * t;
}
long long POW_MOD(long long n, long long p, long long M) {
  if (p == 0) return 1;
  long long t = POW(n, p / 2);
  if (p & 1) return ((((n % M) * (t % M)) % M) * (t % M)) % M;
  return ((t % M) * (t % M)) % M;
}
long long GCD(long long a, long long b) {
  if (b > a) swap(a, b);
  return b == 0 ? a : GCD(b, a % b);
}
long long C(long long n, long long r) {
  if (r == 0) return 1;
  if (r == 1) return n;
  if (r > (n / 2)) r = n - r;
  long long res = 1;
  for (int i = 1; i <= r; i++) {
    res *= (n - i + 1);
    res /= i;
  }
  return res;
}
bool out(int x, int y, int r, int c) {
  if (x < 0 || x >= r) return 1;
  if (y < 0 || y >= c) return 1;
  return 0;
}
long long dp[1005];
int p[1005];
int main() {
  int n;
  while (cin >> n) {
    for (int i = 1; i <= n; i++) cin >> p[i];
    memset(dp, 0, sizeof dp);
    dp[1] = 2;
    for (int i = 2; i <= n; i++) {
      for (int j = p[i]; j < i; j++) {
        dp[i] += dp[j];
        dp[i] %= 1000000007;
      }
      dp[i] += 2;
      dp[i] %= 1000000007;
    }
    long long res = 0;
    for (int i = 1; i <= n; i++) {
      res += dp[i];
      res %= 1000000007;
    }
    printf("%I64d", res);
    printf("\n");
  }
  return 0;
}
