ソースコード
#include <bits/stdc++.h>
using namespace std;
#define F first
#define S second
#define R cin>>
#define Z class
#define ll long long
#define ln cout<<'\n'
#define in(a) insert(a)
#define pb(a) push_back(a)
#define pd(a) printf("%.10f\n",a)
#define mem(a) memset(a,0,sizeof(a))
#define all(c) (c).begin(),(c).end()
#define iter(c) __typeof((c).begin())
#define rrep(i,n) for(ll i=(ll)(n)-1;i>=0;i--)
#define REP(i,m,n) for(ll i=(ll)(m);i<(ll)(n);i++)
#define rep(i,n) REP(i,0,n)
#define tr(it,c) for(iter(c) it=(c).begin();it!=(c).end();it++)
template<Z A>void pr(A a){cout<<a;ln;}
template<Z A,Z B>void pr(A a,B b){cout<<a<<' ';pr(b);}
template<Z A,Z B,Z C>void pr(A a,B b,C c){cout<<a<<' ';pr(b,c);}
template<Z A,Z B,Z C,Z D>void pr(A a,B b,C c,D d){cout<<a<<' ';pr(b,c,d);}
template<Z A>void PR(A a,ll n){rep(i,n){if(i)cout<<' ';cout<<a[i];}ln;}
ll check(ll n,ll m,ll x,ll y){return x>=0&&x<n&&y>=0&&y<m;}
const ll MAX=1e9+7,MAXL=1LL<<61,dx[4]={-1,0,1,0},dy[4]={0,1,0,-1};
typedef pair<int,int> P;
int N=1<<18;
class StarrySkyTree{
public:
ll Mi[555555],A[555555];
void init(){memset(Mi,0,sizeof(Mi)),memset(A,0,sizeof(A));}
void add(int a,int b,ll x,int k=0,int l=0,int r=N) {
if(r<=a||b<=l) return;
if(a<=l&&r<=b){
A[k]+=x;
A[k]%=MAX;
while(k){
k=(k-1)/2;
Mi[k]=min(Mi[k*2+1]+A[k*2+1],Mi[k*2+2]+A[k*2+2]);
Mi[k]%=MAX;
}
return;
}
add(a,b,x,k*2+1,l,(l+r)/2);
add(a,b,x,k*2+2,(l+r)/2,r);
}
ll getMin(int a,int b,int k=0,int l=0,int r=N) {
if(r<=a||b<=l)return LLONG_MAX;if(a<=l&&r<=b)return (Mi[k]+A[k])%MAX;
ll left=getMin(a,b,k*2+1,l,(l+r)/2),right=getMin(a,b,k*2+2,(l+r)/2,r);
return (min(left,right)+A[k])%MAX;
}
};
StarrySkyTree t;
void Main() {
t.init();
int n;
R n;
int a[n];
rep(i,n) R a[i];
sort(a,a+n);
ll d[111111];
mem(d);
rep(i,n) d[a[i]]++;
vector<int> v;
rep(i,111111) {
if(d[i]) v.pb(d[i]);
}
t.add(0,1,1);
rep(i,v.size()) {
ll x=t.getMin(i,i+1)*v[i]%MAX;
t.add(i+1,v.size()+1,x);
}
pr(t.getMin(v.size(),v.size()+1));
}
int main(){ios::sync_with_stdio(0);cin.tie(0);Main();return 0;}