#include // clang-format off using namespace std; using Int = long long; #define REP2(i, n) for (Int i = 0, max_i = (n); i < max_i; i++) #define REP3(i, a, b) for (Int i = (a), max_i = (b); i < max_i; i++) #define OVERLOAD2(_1, _2, _3, name, ...) name #define REP(...) OVERLOAD2(__VA_ARGS__, REP3, REP2)(__VA_ARGS__) struct SetupIO { SetupIO() { cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(13); } } setup_io; #ifndef _MY_DEBUG #define dump(...) #endif // clang-format on /** * author: knshnb * created: Thu Feb 20 21:05:46 JST 2020 **/ template T pow(T x, int n, const T UNION = 1) { T ret = UNION; while (n) { if (n & 1) ret *= x; x *= x; n >>= 1; } return ret; } template struct ModInt { int x; static unordered_map to_inv; ModInt() : x(0) {} ModInt(long long x_) { if ((x = x_ % MD + MD) >= MD) x -= MD; } ModInt& operator+=(ModInt that) { if ((x += that.x) >= MD) x -= MD; return *this; } ModInt& operator*=(ModInt that) { x = (unsigned long long)x * that.x % MD; return *this; } ModInt& operator-=(ModInt that) { if ((x -= that.x) < 0) x += MD; return *this; } ModInt& operator/=(ModInt that) { x = (unsigned long long)x * that.inv().x % MD; return *this; } ModInt operator-() const { return -x < 0 ? MD - x : -x; } ModInt operator+(ModInt that) const { return ModInt(*this) += that; } ModInt operator*(ModInt that) const { return ModInt(*this) *= that; } ModInt operator-(ModInt that) const { return ModInt(*this) -= that; } ModInt operator/(ModInt that) const { return ModInt(*this) /= that; } bool operator==(ModInt that) const { return x == that.x; } bool operator!=(ModInt that) const { return x != that.x; } ModInt inv() const { return to_inv.count(this->x) ? to_inv[this->x] : (to_inv[this->x] = pow(*this, MD - 2).x); } friend ostream& operator<<(ostream& s, ModInt a) { s << a.x; return s; } friend istream& operator>>(istream& s, ModInt& a) { long long tmp; s >> tmp; a = ModInt(tmp); return s; } friend string to_string(ModInt a) { return to_string(a.x); } }; template unordered_map ModInt::to_inv; using mint = ModInt<1000000007>; vector fact, fact_inv; void init_factorial(int n) { fact = vector(n + 1, 1); fact_inv = vector(n + 1); for (int i = 0; i < n; i++) fact[i + 1] = fact[i] * (i + 1); fact_inv[n] = mint(1) / fact[n]; for (int i = n - 1; i >= 0; i--) fact_inv[i] = fact_inv[i + 1] * (i + 1); // for (int i = 0; i < n + 1; i++) assert(fact[i] * fact_inv[i] == 1); } mint comb(int n, int r) { return fact[n] * fact_inv[r] * fact_inv[n - r]; } template struct SegmentTree { const F op; const T e; SegmentTree(F op_, T e_) : op(op_), e(e_) {} int n; vector t; void set_by_identity(int n_) { n = n_; t.clear(), t.resize(2 * n, e); } void set_by_vector(const vector& a) { n = a.size(); t.clear(), t.resize(2 * n, e); for (int i = 0; i < n; i++) t[i + n] = a[i]; build(); } void build() { for (int i = n - 1; i; --i) t[i] = op(t[2 * i], t[2 * i + 1]); } T& operator[](int i) { return t[i + n]; } // [l, r) T query(int l, int r) const { assert(0 <= l && l <= r && r <= n); T resl = e, resr = e; for (l += n, r += n; l < r; l >>= 1, r >>= 1) { if (l & 1) resl = op(resl, t[l++]); if (r & 1) resr = op(t[--r], resr); } return op(resl, resr); } // iをaに変更 void update(int i, const T& x) { assert(0 <= i && i < n); i += n; t[i] = op(t[i], x); while (i >>= 1) t[i] = op(t[2 * i], t[2 * i + 1]); } friend string to_string(const SegmentTree& seg) { return to_string(vector(seg.t.begin() + seg.n, seg.t.end())); } }; template auto make_segment_tree(F op, T e) { return SegmentTree(op, e); } template struct Compressor { vector sorted; Compressor(const vector& a) : sorted(a) { sort(sorted.begin(), sorted.end()); sorted.erase(unique(sorted.begin(), sorted.end()), sorted.end()); } int get(T x) { auto it = lower_bound(sorted.begin(), sorted.end(), x); assert(*it == x); return it - sorted.begin(); } }; signed main() { Int n; cin >> n; vector a(n); REP(i, n) cin >> a[i]; Compressor cp(a); // {max, max_num} auto seg = make_segment_tree>( [](pair a, pair b) -> pair { mint num = 0; if (a.first >= b.first) num += a.second; if (a.first <= b.first) num += b.second; return {max(a.first, b.first), num}; }, {-1e18, 0}); seg.set_by_vector(vector>(n + 1, {0, 0})); Int ma = -1; REP(i, n) { Int x = cp.get(a[i]) + 1; auto prv = seg.query(0, x); prv.first++; if (prv.second == 0) prv.second += 1; seg.update(x, prv); ma = max(ma, seg[x].first); } mint ans = 0; REP(i, n + 1) { if (seg[i].first == ma) ans += seg[i].second; } cout << ans << endl; }