#include<cstdio>
#define MOD 1000000007

using namespace std;

long long vec[1300],mat[600][1300],in[1300],n,l,c,i,j,k,t,sum;

int main(){
    scanf("%lld",&n);
    for(i=0;i<n;i++){
        scanf("%lld",&in[i]);
    }
    vec[0]=1;
    c=1;
    for(i=n-1;i>0;i--){
        if(in[i]<in[i-1]){
            for(j=0;j*2-1<=c;j++)mat[j][0]=vec[j];
            for(j=1;j<c;j++){
                for(k=1;k<=j;k++){
                    mat[k][j]=0;
                    if(j+k>c)break;
                    t=1;
                    for(l=0;l<=k;l++){
                        if(j-k<l&&j!=k)break;
                        mat[k][j]+=(mat[l][j-k])*t;
                        mat[k][j]%=MOD;
                        t*=k-l;
                        t%=MOD;
                    }
                }
            }
            sum=0;
            t=1;
            for(j=0;j<=c;j++){
                sum+=vec[j]*t;
                sum%=MOD;
                t*=c-j;
                t%=MOD;
            }
            for(j=0;j<=n;j++)vec[j]=0;
            vec[c]=sum;
            for(j=c-1;j>0;j--){
                vec[j]=0;
                t=1;
                for(k=0;k<=j;k++){
                    if(c-j<k)break;
                    vec[j]+=mat[k][c-j]*t;
                    vec[j]%=MOD;
                    t*=j-k;
                    t%=MOD;
                }
            }
            vec[0]=0;
            c=1;
        }
        else c++;
    }
    for(j=0;j*2-1<=c;j++)mat[j][0]=vec[j];
            for(j=1;j<c;j++){
                for(k=1;k<=j;k++){
                    mat[k][j]=0;
                    if(j+k>c)break;
                    t=1;
                    for(l=0;l<=k;l++){
                        if(j-k<l&&j!=k)break;
                        mat[k][j]+=(mat[l][j-k])*t;
                        mat[k][j]%=MOD;
                        t*=k-l;
                        t%=MOD;
                    }
                }
            }
            sum=0;
            t=1;
            for(j=0;j<=c;j++){
                sum+=vec[j]*t;
                sum%=MOD;
                t*=c-j;
                t%=MOD;
            }
            for(j=0;j<=n;j++)vec[j]=0;
            vec[c]=sum;
            for(j=c-1;j>0;j--){
                vec[j]=0;
                t=1;
                for(k=0;k<=j;k++){
                    if(c-j<k)break;
                    vec[j]+=mat[k][c-j]*t;
                    vec[j]%=MOD;
                    t*=j-k;
                    t%=MOD;
                }
            }
            vec[0]=0;
    sum=0;
    for(i=0;i<=c;i++){
        sum+=vec[i];
        sum%=MOD;
    }
    printf("%lld\n",sum);
    return 0;
}