#include <vector>
#include <string> 
#include <set> 
#include <map> 
#include <unordered_set>
#include <unordered_map>
#include <queue> 
#include <bitset> 
#include <stack>
#include <list>

#include <numeric> 
#include <algorithm> 
#include <random>
#include <chrono>

#include <cstdio>
#include <fstream>
#include <iostream> 
#include <sstream> 
#include <iomanip>
#include <climits>

#include <cctype>
#include <cmath> 
#include <ctime>
#include <cassert>

using namespace std;

#define ULL unsigned long long
#define LL long long
#define PII pair <int, int>
#define VB vector <bool>
#define VI vector <int>
#define VLL vector <LL>
#define VD vector <double>
#define VS vector <string>
#define VPII vector <pair <int, int> >
#define VVI vector < VI >
#define VVB vector < VB >
#define SI set < int >
#define USI unordered_set <int>
#define MII map <int, int>
#define UMII unordered_map <int, int>
#define MS multiset
#define US unordered_set
#define UM unordered_map
#define UMS unordered_multiset
#define UMM unordered_multimap

#define FORN(i, n) for(int i = 0; i < (n); ++i)
#define FOR(i, a, b) for(int i = (a); i <= (b); ++i)
#define FORD(i, a, b) for(int i = (a); i >= (b); --i)

#define SZ size()
#define BG begin() 
#define EN end() 
#define CL clear()
#define X first
#define Y second
#define RS resize
#define PB push_back
#define MP make_pair
#define ALL(x) x.begin(), x.end()
#define INS insert
#define ER erase
#define CNT count

#define IN_FILE "a.in"
#define OUT_FILE "a.out"

template <typename T>
void PR(T var1)
{
	cout << var1 <<endl;
}
template <typename T, typename... Types>
void PR(T var1, Types... var2)
{
	cout << var1;
	PR(var2...);
}

#define MOD 1000000007

int n;
VI a, pow2, cnt(100001), part;
map<PII, int> memo;

int go(int s, int k)
{
	if (memo.CNT({s, k})) return memo[{s, k}];
	if (k == 1) return memo[{s, k}] = (s <= part[1]);
	if (s == 0) return memo[{s, k}] = 1;
	int ans = go(s - 1, k);
	if (s > part[k]) ans -= go(s - part[k] - 1, k - 1);
	ans += go(s, k - 1);
	if (ans < MOD) ans += MOD;
	return memo[{s, k}] = ans % MOD;
}

int main()
{
 	cin >> n;
	a.RS(n + 1);
	FOR(i, 1, n) 
	{
		cin >> a[i];
		++cnt[a[i]];
	}
	FOR(i, 0, 100000)
		if (cnt[i]) part.PB(cnt[i]);

	int sum = 0;
	FOR(i, 0, (n - part[0])) 
	{
		sum += go(i, part.SZ - 1);
		sum %= MOD;
	}

	PR(sum);

	return 0;
}