#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<vector>
#include<cmath>
#include<algorithm>
#include<map>
#include<queue>
#include<deque>
#include<iomanip>
#include<tuple>
#include<cassert>
using namespace std;
typedef long long int LL;
typedef pair<int,int> P;
typedef pair<LL,int> LP;
const int INF=1<<30;
const LL MAX=1e9+7;

void array_show(int *array,int array_n,char middle=' '){
	for(int i=0;i<array_n;i++)printf("%d%c",array[i],(i!=array_n-1?middle:'\n'));
}
void array_show(LL *array,int array_n,char middle=' '){
	for(int i=0;i<array_n;i++)printf("%lld%c",array[i],(i!=array_n-1?middle:'\n'));
}
void array_show(vector<int> &vec_s,int vec_n=-1,char middle=' '){
	if(vec_n==-1)vec_n=vec_s.size();
	for(int i=0;i<vec_n;i++)printf("%d%c",vec_s[i],(i!=vec_n-1?middle:'\n'));
}
void array_show(vector<LL> &vec_s,int vec_n=-1,char middle=' '){
	if(vec_n==-1)vec_n=vec_s.size();
	for(int i=0;i<vec_n;i++)printf("%lld%c",vec_s[i],(i!=vec_n-1?middle:'\n'));
}

class union_find_tree{
private:
	static constexpr int uft_N=100005;
	int uft_n;
	queue<int> uft_q1;
	vector<int> uft_parent;
	vector<int> uft_num;
public:
	void init(){
		uft_parent.assign(uft_n,-1);
		uft_num.assign(uft_n,1);
	}
	union_find_tree(int uft_num){
		assert(uft_num>=0);
		uft_n=uft_num;
		init();
	}
	union_find_tree(){
		uft_n=uft_N;
		init();
	}

	int check_parent(int uft_x){
		assert(uft_x>=0 && uft_x<uft_n);
		if(uft_parent[uft_x]!=-1){
			uft_q1.push(uft_x);
			return check_parent(uft_parent[uft_x]);
		}
		int uft_a;
		while(!uft_q1.empty()){
			uft_a=uft_q1.front(),uft_q1.pop();
			uft_parent[uft_a]=uft_x;
		}
		return uft_x;
	}

	int size(int pos){
		pos=check_parent(pos);
		return uft_num[pos];
	}

	bool connect(int uft_x,int uft_y){
		assert(uft_x>=0 && uft_x<uft_n);
		assert(uft_y>=0 && uft_y<uft_n);
		uft_x=check_parent(uft_x),uft_y=check_parent(uft_y);
		if(uft_x==uft_y)return true;
		if(uft_num[uft_x]>uft_num[uft_y])swap(uft_x,uft_y);
		uft_parent[uft_x]=uft_y;
		uft_num[uft_y]+=uft_num[uft_x];
		return false;
	}
};

long long int pow_mod(long long int p_a,long long int p_n,long long int p_p=1e9+7){
	//p_a^p_n mod p_p
	long long int p_b=1,p_t=1;
	for(;p_b<=p_n;p_b<<=1);
	for(p_b>>=1;p_b>0;p_b>>=1){
		p_t*=p_t;
		if(p_t>=p_p)p_t%=p_p;
		if(p_n&p_b)p_t*=p_a;
		if(p_t>=p_p)p_t%=p_p;
	}
	return p_t;
}

long long int gcd(long long int g_a,long long int g_b){
	if(g_a<g_b)gcd(g_b,g_a);
	if(g_b==0)return g_a;
	return gcd(g_b,g_a%g_b);
}

long long int divide(long long int d_a,long long int d_b,long long int d_p=1e9+7){
	// d_a/d_b mod d_p
	// prime:d_p is prime
	if(d_a>=d_p)d_a%=d_p;
	if(d_a<0)d_a+=d_p;
	if(d_b>=d_p)d_b%=d_p;
	if(d_b<0)d_b+=d_p;
	d_a*=pow_mod(d_b,d_p-2,d_p);
	return d_a%d_p;
}

int t[5500];
LL ta[5500];

int main(){
	int n,m;
	int i,j,k;
	LL a,b,c;
	LL s=0;
	cin>>n;
	union_find_tree ua(n);
	vector<int> v1;
	for(i=0,a=0;i<n;i++){
		cin>>t[i];
		t[i]--;
		if(t[i]!=-2){
			if(!ua.connect(i,t[i]))s++;
		}else{
			a++;
			v1.push_back(i);
		}
	}
	LL p,q,r;
	p=pow_mod(n-1,a);
	s*=p;
	s%=MAX;
	s+=(LL)v1.size()*pow_mod(n-1,a);
	s%=MAX;
	q=v1.size();
	ta[0]=1;
	for(i=0;i<q;i++){
		a=ua.size(v1[i]);
		for(j=q;j>0;j--){
			ta[j]+=ta[j-1]*a;
			if(ta[j]>=MAX)ta[j]%=MAX;
		}
	}
	ta[1]-=q;
	p=1;
	for(i=0;i<q;i++){
		if(i)p*=i;
		if(p>=MAX)p%=MAX;
		a=ta[i+1]*pow_mod(n-1,q-i-1)%MAX;
		a=a*p%MAX;
		s-=a;
		if(s<0)s+=MAX;
	}
	cout<<s<<endl;
}

Measured memory_bytes: 3760