#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,avx2") #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define overload4(_1, _2, _3, _4, name, ...) name #define rep0(a) for (ll _ = 0; _ < ll(a); ++_) #define rep1(i, n) for (ll i = 0; i < ll(n); ++i) #define rep2(i, s, n) for (ll i = ll(s); i < ll(n); ++i) #define rep3(i, s, n, d) for(ll i = ll(s); i < ll(n); i+=d) #define rep(...) overload4(__VA_ARGS__,rep3,rep2,rep1,rep0)(__VA_ARGS__) #define rrep0(a) for (ll _ = (a)-1; _ >= ll(0); -- _) #define rrep1(i, n) for (ll i = ll(n)-1; i >= 0; i--) #define rrep2(i, n, t) for (ll i = ll(n)-1; i >= (ll)t; i--) #define rrep3(i, n, t, d) for (ll i = ll(n)-1; i >= (ll)t; i-=d) #define drep(...) overload4(__VA_ARGS__,rrep3,rrep2,rrep1,rrep0)(__VA_ARGS__) typedef long long ll; typedef unsigned long long ull; typedef long double LD; typedef double D; typedef pair P; typedef map M; // /* #include using namespace atcoder; //using namespace internal; using mint =modint998244353; //using mint =modint1000000007; //using mint=static_modint<4649>; #define ip(x) atcoder::internal::is_prime_constexpr(x) istream &operator>>(istream &is, mint &a) { int v; cin >> v; a = v; return is; } ostream &operator<<(ostream &os, const mint &a) { return os << a.val(); } auto v_pow(ll n,ll base){vector v(n,1);rep(i,n-1){v[i+1]*=base*v[i];}return v;} // */ template istream &operator>>(istream &is, vector &v) { for (auto &e : v) is >> e; return is; } template ostream &operator<<(ostream &os, const vector &v) { for (auto &e : v) os << e << ' '; return os; } template istream &operator>>(istream &is, pair &p) { return is >> p.first >> p.second; } template ostream &operator<<(ostream &os, const pair &p) { return os << '{' << p.first << ", " << p.second << '}'; } template istream &operator>>(istream &is, tuple &t) { return is >> get<0>(t) >> get<1>(t) >> get<2>(t); } template ostream &operator<<(ostream &os, const tuple &t) {return os << '{' << get<0>(t) << ", " << get<1>(t) << ", " << get<2>(t) << '}';} #define YES(n) cout << ((n) ? "YES" : "NO" ) << endl #define Yes(n) cout << ((n) ? "Yes" : "No" ) << endl template ll binary_search(F check,ll ok,ll ng) { /*binary_search(bool_func,okな値,ngな値)*///しっかり区間内に... assert(check(ok)); while(abs(ok-ng)>1){ auto x=(ng+ok)/2; tie(ok,ng)=(check(x)?make_pair(x, ng):make_pair(ok, x)); } return ok; } template double binary_search_real(F check,double ok,double ng,ll iter = 100) { rep(iter) { double x=(ok+ng)/2; tie(ok,ng)=(check(x)?make_pair(x,ng):make_pair(ok,x)); } return (ok+ng)/2; } void tatananonano() { ios::sync_with_stdio(false); std::cin.tie(nullptr); cout<< fixed << setprecision(10); } #define LL(...) ll __VA_ARGS__; IN(__VA_ARGS__) #define INT(...) int __VA_ARGS__; IN(__VA_ARGS__) #define STR(...) string __VA_ARGS__; IN(__VA_ARGS__) #define CHR(...) char __VA_ARGS__;IN(__VA_ARGS__) #define LDL(...) LD __VA_ARGS__;IN(__VA_ARGS__) #define vv(type, name, h, ...) \ vector> name(h, vector(__VA_ARGS__)) #define vvv(type, name, h, w, ...) \ vector>> name( \ h, vector>(w, vector(__VA_ARGS__))) //#define MINT(...) mint __VA_ARGS__;IN(__VA_ARGS__) template void scan(T &a) { cin >> a; } void IN() {} template void IN(Head &head, Tail &...tail) {scan(head);IN(tail...);} #define overload2(_1, _2, _3,name, ...) name #define out1(x) cout<,greater

>//小さい順 #define PQS PQ> #define fi first #define se second #define bit(n,k) ((n>>k)&1LL) #define popcount(n) __builtin_popcountll(n) template inline bool chmax(T& a,T b){if(a < b){a=b;return 1;}return 0;} template inline bool chmin(T& a,T b){if(a > b){a=b;return 1;}return 0;} bool in_rect(ll i,ll j,ll h,ll w) {return 0 <= i and i < h and 0 <= j and j < w;} typedef vector vec; typedef vector vs; typedef vector mat; const ll mod = 998244353; //const ll mod = 1000000007; const auto INF = (1LL<<(60)); template T ceil(T x, U y) {return (x > 0 ? (x + y - 1) / y : x / y);} template T floor(T x, U y) {return (x > 0 ? x / y : (x - y + 1) / y);} auto longest_increasing_subsequence(const vec &a, bool strict) { ll n=a.size(); vec lis,num(n); rep(i,n) { typename vec::iterator it; if(strict) it = lower_bound(begin(lis), end(lis), a[i]); else it = upper_bound(begin(lis), end(lis), a[i]); if(end(lis) == it) { lis.emplace_back(a[i]); num[i]=lis.size()-1;} else { *it = a[i]; num[i]=it-lis.begin();} } ll LIS=lis.size(); vec ans(LIS); drep(i,n){ if(LIS-1==num[i]){ ans[LIS-1]=i;//a[i] LIS--; } } return num; } vec shrink(vec a){ vec b=a; uni(b); vec res(a.size()); rep(i,a.size()){ res[i]=lb(all(b),a[i])-b.begin(); } return res; } int main(){ tatananonano(); LL(n); vec p(n); cin>>p; auto a=shrink(p); auto v=longest_increasing_subsequence(p,true); ll key=*max_element(all(v)); vector> fw(key+1,fenwick_tree(n+1)); vec dp(n+1); ll ans=0; rep(i,n){ if(v[i]!=0){dp[i+1]=fw[v[i]-1].sum(0,a[i]);} else{dp[i+1]=1;} fw[v[i]].add(a[i],dp[i+1]); if(v[i]==key){ans+=dp[i+1];} } cout<