#pragma GCC optimize("Ofast") //#pragma GCC optimize("unroll-loops") //#pragma GCC push_options //#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native") #include #include #include using namespace::std; __attribute__((constructor))void init(){cin.tie(0);ios::sync_with_stdio(false);cout< #include #include #include // #include // #include // namespace mp = boost::multiprecision; // typedef mp::number> cdouble; // typedef mp::cpp_int cint; templateusing pbds=__gnu_pbds::tree,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>; templateusing pbds_map=__gnu_pbds::tree,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>; templateusing hash_map=__gnu_pbds::gp_hash_table; templateusing pqueue =__gnu_pbds::priority_queue,__gnu_pbds::rc_binomial_heap_tag>; typedef long long lint; #define INF (1LL<<32) #define IINF (1<<30) #define EPS (1e-10) #define MOD 1000000007LL //#define MOD 998244353LL typedef vector vec; typedef vector> mat; typedef vector>> mat3; typedef vector svec; typedef vector> smat; templateinline void numout(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<inline void numout2(T t){for(auto i:t)numout(i);} templateinline void output(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<inline void output2(T t){for(auto i:t)output(i);} templateinline void _output(T t){bool f=0;for(lint i=0;iinline void _output2(T t){for(lint i=0;i=0;--i) #define rrepi(i,a,b) for(lint i=lint(b)-1;i>=lint(a);--i) #define irep(i) for(lint i=0;;++i) #define all(n) begin(n),end(n) #define dist(a,b,c,d) sqrt(pow(a-c,2)+pow(b-d,2)) inline lint gcd(lint A,lint B){return B?gcd(B,A%B):A;} inline lint lcm(lint A,lint B){return A/gcd(A,B)*B;} // inline cint cgcd(cint A,cint B){return B?cgcd(B,A%B):A;} // inline cint clcm(cint A,cint B){return A/cgcd(A,B)*B;} //inline bool chmin(auto& s,const auto& t){bool res=s>t;s=min(s,t);return res;} //inline bool chmax(auto& s,const auto& t){bool res=s dx={1,0,-1,0,1,1,-1,-1}; const vector dy={0,1,0,-1,1,-1,1,-1}; #define SUM(v) accumulate(all(v),0LL) //auto call=[&](auto f,auto... args){return f(f,args...);}; class mint { using u64 = std::uint_fast64_t; public: u64 a; constexpr mint(const u64 x = 0)noexcept:a(x % MOD){} constexpr u64 &value()noexcept{return a;} constexpr const u64 &value() const noexcept {return a;} constexpr mint operator+(const mint rhs)const noexcept{return mint(*this) += rhs;} constexpr mint operator-(const mint rhs)const noexcept{return mint(*this)-=rhs;} constexpr mint operator*(const mint rhs) const noexcept {return mint(*this) *= rhs;} constexpr mint operator/(const mint rhs) const noexcept {return mint(*this) /= rhs;} constexpr mint &operator+=(const mint rhs) noexcept { a += rhs.a; if (a >= MOD)a -= MOD; return *this; } constexpr mint &operator-=(const mint rhs) noexcept { if (a= MOD)a -= MOD; return *this; } constexpr mint operator--(int n) noexcept { if (a<1)a += MOD; a -= 1; return *this; } constexpr mint &operator/=(mint rhs) noexcept { u64 exp=MOD-2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } mint pow(long long n){ long long x=a; mint ret = 1; while(n>0) { if(n&1)(ret*=x); (x*=x)%=MOD; n>>=1; } return ret; } mint inv()const{ int _a=a,b=MOD,u=1,v=0,t; while(b>0) { t=_a/b; swap(_a-=t*b,b); swap(u-=t*v,v); } return mint(u); } mint comb(lint b){ using lint=long long; static bool init=1; static lint fac[3000001],ifac[3000001]; if(init){ init=0; fac[0]=1; ifac[0]=1; auto mod_pow=[&](lint x,lint n){ lint ans = 1; while(n != 0){ if(n&1)ans=ans*x%MOD; x=x*x%MOD; n=n>>1; } return ans; }; for(int i=0;i<3000000;i++){ fac[i+1]=fac[i]*(i+1)%MOD; ifac[i+1]=ifac[i]*mod_pow(i+1, MOD-2)%MOD; } } if(a==0&&b==0)return 1; if((lint)a>(istream& lhs,mint& rhs) noexcept { lhs >> rhs.a; return lhs; } }; inline mint fact(lint a){return a?fact(a-1)*a:1;} template struct dynamic_lazy_segment{ private: T et=0; E ee=0; lint n=1LL<<32; lint size=1LL<<32; struct node{ T val; E lazy; node *lch,*rch; node(dynamic_lazy_segment &seg):lch(nullptr),rch(nullptr){val=seg.et;lazy=seg.ee;} }; node* root=new node(*this); public: dynamic_lazy_segment(lint size):n(1),size(size){ while(nval; //全部区間内 if(a<=l&&r<=b){ t->lazy=g(t->lazy,x); eval(t,r-l,k); return h(t->val,t->lazy,r-l); } //一部区間内 if(!t->lch)t->lch=new node(*this); if(!t->rch)t->rch=new node(*this); return t->val=f(update(a,b,x,2*k+1,l,(l+r)/2,t->lch),update(a,b,x,2*k+2,(l+r)/2,r,t->rch)); } //0-indexed [a,b) T get(const lint& a,const lint& b,const lint& k,const lint& l,const lint& r,node* t){ eval(t,r-l,k); //区間外 if(r<=a||b<=l)return et; //全部区間内 if(a<=l&&r<=b)return t->val; //一部区間内 if(!t->lch)t->lch=new node(*this); if(!t->rch)t->rch=new node(*this); return f(get(a,b,2*k+1,l,(l+r)/2,t->lch),get(a,b,2*k+2,(l+r)/2,r,t->rch)); } inline void eval(node* t,const lint& len,const lint& k){ if(t->lazy==ee) return; if(k*2+1lch==nullptr)t->lch=new node(*this); if(t->rch==nullptr)t->rch=new node(*this); t->lch->lazy=g(t->lch->lazy,t->lazy); t->rch->lazy=g(t->rch->lazy,t->lazy); } t->val=h(t->val,t->lazy,len); t->lazy=ee; } T f(T a,T b){ return a+b; } T h(T a,E b,lint len){ return a+b*len; } E g(E a,E b){ return b+a; } //auto update(auto a,auto b){return b==ee?a:b;} }; int main(){ lint n; cin>>n; vec a(n); rep(i,n)cin>>a[i]; vec dp; vectordp2(n); rep(i,n){ auto itr=lower_bound(all(dp),a[i]); if(itr==dp.end()){ dp.push_back(a[i]); dp2[i]=dp.size()-1; } else { *itr=a[i]; dp2[i]=itr-dp.begin(); } } hash_map>m; rep(i,n){ if(dp2[i]==0)m[0].update(a[i],a[i]+1,1); else{ m[dp2[i]].update(a[i],a[i]+1,m[dp2[i]-1].get(-INF,a[i]).a); } } cout<