#pragma GCC optimize("Ofast")
#define _USE_MATH_DEFINES
#include "bits/stdc++.h"

using namespace std;

using u8 = uint8_t;
using u16 = uint16_t;
using u32 = uint32_t;
using u64 = uint64_t;
using i8 = int8_t;
using i16 = int16_t;
using i32 = int32_t;
using i64 = int64_t;

constexpr char newl = '\n';
constexpr double eps = 1e-10;

#define FOR(i,a,b) for (int i = (a); i < (b); i++)
#define F0R(i,b) FOR(i,0,b)
#define RFO(i,a,b) for (int i = ((b)-1); i >=(a); i--)
#define RF0(i,b) RFO(i,0,b)
#define fi first
#define se second
#define debug(x) cout << #x << ": " << x << '\n';
#define rng(a) a.begin(),a.end()
#define rrng(a) a.rbegin(),a.rend()

template<typename T1, typename T2> inline void chmin(T1& a, T2 b) { if (a > b) a = b; }
template<typename T1, typename T2> inline void chmax(T1& a, T2 b) { if (a < b) a = b; }
template<class T> void Print(vector<T> v) {
  F0R(i, v.size()) {
    cout << v[i] << ' ';
  }
  cout << newl;
}

#if 1

constexpr int MOD = 1e9 + 7; // 1e9 + 7

template<int MOD>
struct Modint {
  long long val;
  constexpr Modint(long long v = 0) noexcept : val(v% MOD) {
    if (val < 0) val += MOD;
  }
  constexpr int getmod() { return MOD; }
  constexpr Modint operator - () const noexcept {
    return val ? MOD - val : 0;
  }
  constexpr Modint operator + (const Modint& r) const noexcept { return Modint(*this) += r; }
  constexpr Modint operator - (const Modint& r) const noexcept { return Modint(*this) -= r; }
  constexpr Modint operator * (const Modint& r) const noexcept { return Modint(*this) *= r; }
  constexpr Modint& operator += (const Modint& r) noexcept {
    val += r.val;
    if (val >= MOD) val -= MOD;
    return *this;
  }
  constexpr Modint& operator -= (const Modint& r) noexcept {
    val -= r.val;
    if (val < 0) val += MOD;
    return *this;
  }
  constexpr Modint& operator *= (const Modint& r) noexcept {
    val = val * r.val % MOD;
    return *this;
  }
  constexpr bool operator == (const Modint& r) const noexcept {
    return this->val == r.val;
  }
  constexpr bool operator != (const Modint& r) const noexcept {
    return this->val != r.val;
  }
  friend constexpr ostream& operator << (ostream& os, const Modint<MOD>& x) noexcept {
    return os << x.val;
  }
  friend constexpr istream& operator >> (istream& is, Modint<MOD>& x) noexcept {
    return is >> x.val;
  }
  
  constexpr Modint<MOD> pow(long long n) const noexcept {
    if (!n) return 1;
    auto t = pow(n >> 1);
    t *= t;
    if (n & 1) t *= *this;
    return t;
  }
  friend constexpr Modint<MOD> modpow(const Modint<MOD>& a, long long n) noexcept {
    return a.pow(n);
  }

  friend constexpr Modint<MOD> Ncr(Modint<MOD> a, int n) {
    Modint<MOD> p = 1, q = 1;

    while (n > 0) {
      p *= a;
      a -= 1;
      q *= n--;
    }

    return p / q;
  }

  // for prime mod
  constexpr Modint& operator /= (const Modint& r) noexcept {
    long long a = r.val, b = MOD, u = 1, v = 0;
    while (b) {
      long long t = a / b;
      a -= t * b; swap(a, b);
      u -= t * v; swap(u, v);
    }
    val = val * u % MOD;
    if (val < 0) val += MOD;
    return *this;
  }
  constexpr Modint operator / (const Modint& r) const noexcept { return Modint(*this) /= r; }
};
template<int MOD>
using mint = Modint<MOD>;

class UnionFind {
public:
  UnionFind(int size) : size(size) {

    par = vector<int>(size, -1);
  }

  int GetRoot(int v) {
    int root = v;
    while (par[root] >= 0) {
      root = par[root];
    }
    if (root != v) {
      par[v] = root;
    }
    return root;
  }

  bool IsUnion(int a, int b) {
    int ra = GetRoot(a), rb = GetRoot(b);

    return ra == rb && ra >= 0;
  }

  int GetUnionSize(int a) {
    return -par[GetRoot(a)];
  }

  void Unite(int a, int b) {
    int ra = GetRoot(a), rb = GetRoot(b);

    if (ra == rb) return;

    if (GetUnionSize(ra) > GetUnionSize(rb)) {
      par[ra] += par[rb];
      par[rb] = ra;
    }
    else {
      par[rb] += par[ra];
      par[ra] = rb;
    }
  }

protected:

  // 親がいれば、親のindexを入れる
  // 親がいなければ、グループのサイズを負にして入れておく（一人なら-1）
  vector<int> par;

  int size;
};

// INSERT ABOVE HERE
signed main() {
  cin.tie(0);
  ios_base::sync_with_stdio(false);

  int N;
  cin >> N;

  UnionFind uf(N);

  vector<int> ps(N);

  int k = 0;
  F0R(i, N) {
    cin >> ps[i];
    if (ps[i] > 0) {
      ps[i]--;
      uf.Unite(i, ps[i]);
    }
    else {
      k++;
    }
  }

  map<int,bool> roots;
  F0R(i, N) {
    roots[uf.GetRoot(i)] |= ps[i] < 0;
  }

  // -1 を含む根付木
  vector<int> ns;

  // -1 を含まない木（既に閉路）
  vector<int> ms;

  for (auto r:roots) {
    int s = uf.GetUnionSize(r.first);
    if (r.second) ns.push_back(s);
    else ms.push_back(s);
  }

  //Print(ns);
  
  // 辺の数 = N - サイクル数（閉路の数）

  // i個の頂点集合 S1 から Si までで1つのサイクルを作るパターンは、何回あるか
  // （ただしこの集合は-1である頂点を含むものに限定する）

  // -1 である頂点数を k とする
  
  //  i
  //  ∏ Sv * (i-1)! * (N-1)^(k-i)
  // v=1

  // ただしv=1の計算のとき、Svはそれぞれ-1する（自分自身には接続できないため）

  // (i-1)! * (N-1)^(k-i) の部分は i に対して固定

  // dp[i] = i個の根付木（頂点集合）を選ぶ時のΠSvの和

  // dp[i][j] i個の根付木を選ぶとき、j個目の根付木まで見たときの和

  vector<mint<MOD>> cp(ns.size() + 1);

  mint<MOD> cycles = 0;

  if (k > 0) {
    mint<MOD> sum = 0;
    F0R(i, ns.size()) {
      sum += ns[i] - 1; // 自分自身には繋げられないため
    }
    cycles = sum * 1 * mint<MOD>(N - 1).pow(k - 1);

    F0R(i, ns.size()) {
      cp[i + 1] = cp[i] + ns[i];
    }

    mint<MOD> fact = 1;
    FOR(i, 1, ns.size()) {
      //Print(cp);

      vector<mint<MOD>> dp(ns.size() + 1);

      FOR(j, i, ns.size()) {
        dp[j + 1] = dp[j] + cp[j] * ns[j];
      }

      fact *= i;

      cycles += dp.back() * fact * mint<MOD>(N - 1).pow(k - i - 1);
      //debug(dp.back().val);

      cp = move(dp);
    }
  }
  mint<MOD> result = mint<MOD>(N - ms.size()) * mint<MOD>(N - 1).pow(k);
  result -= cycles;

  cout << result;
}
#endif


Measured memory_bytes: 3964