#include <iostream>
#include <string>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <math.h>
#include <numeric>
#include <complex>
#include <queue>
#include <deque>
#include <stack>
#include <map>
#include <unordered_map>
#include <set>
#include <bitset>
#include <functional>
#include <assert.h>
using ll=long long;
#define int ll
#define rng(i,a,b) for(int i=int(a);i<int(b);i++)
#define rep(i,b) rng(i,0,b)
#define gnr(i,a,b) for(int i=int(b)-1;i>=int(a);i--)
#define per(i,b) gnr(i,0,b)
#define eb emplace_back
#define all(x) x.begin(),x.end()
#define SZ(x) ((int)x.size())
#define fi first
#define se second
const int INF = 1e15;
const int MOD = 1e9+7;
using namespace std;

template<class T> using vc=vector<T>;
template<class T> using vvc=vc<vc<T>>;

using pi=pair<int,int>;
using vi=vc<int>;

template<class T>
inline bool chmax(T &a,T b){if(a<b){a=b;return true;}return false;}
template<class T>
inline bool chmin(T &a,T b){if(a>b){a=b;return true;}return false;}

int power(int x,int y,int p){
	int res=1;
	x=x%p;
	while(y){
		if(y&1)res=(res*x)%p;
		y=y>>1;
		x=(x*x)%p;
	}
	return res;
}

signed main() {
    cin.tie(0);
    ios::sync_with_stdio(0);
    cout<<fixed<<setprecision(20);

    int N; cin >> N;
    int MAX=1000000000;
    int A[N]; 
    rep(i,N) {
        cin >> A[i];
        if(A[i]==0) {
            cout << 0 << endl;
            return 0;
        }
    }
    int sum=1;
    int l=0;
    vc<pi> V(N);
    rep(i,N){
        V[i].fi=0;
        V[i].se=N-1;
    }
    for(int r=0;r<N;r++) {
        sum*=A[r];
        while(sum>MAX && l<=r) {
            sum/=A[l++];
            V[l-1].se=r-1;
        }
        V[r].fi=l;
    }
    int ans=1;
    rep(i,N) {
        ans*=power(A[i],(i-V[i].fi+1)*(V[i].se-i+1),MOD);
        ans%=MOD;
    }
    cout << ans << endl;
}