#include using namespace std; using Int = long long; struct FullyIndexableDictionary{ Int len,blk; vector bit,sum; FullyIndexableDictionary(){} FullyIndexableDictionary(Int len) :len(len),blk((len+31)>>5),bit(blk,0),sum(blk,0){} void set(Int k){ bit[k>>5]|=1<<(k&31); } void build(){ sum[0]=0; for(Int i=1;i>5]>>(k&31))&1); } Int rank(Int k){ return sum[k>>5]+__builtin_popcount(bit[k>>5]&((1<<(k&31))-1)); } Int rank(bool v,Int k){ return (v?rank(k):k-rank(k)); } Int select(bool v,Int k){ if(k<0||rank(v,len)<=k) return -1; Int low=0,high=len; while(low+1>1; if(rank(v,mid)>=k+1) high=mid; else low=mid; } return high-1; } Int select(bool v,Int i,Int l){ return select(v,i+rank(v,l)); } }; template struct WaveletMatrix{ Int len; FullyIndexableDictionary mat[MAXLOG]; Int zs[MAXLOG],buff1[MAXLOG],buff2[MAXLOG]; static const T npos=-1; Int freq_dfs(Int d,Int l,Int r,T val,T a,T b){ if(l==r) return 0; if(d==MAXLOG) return (a<=val&&val data){ len=data.size(); vector l(len),r(len); for(Int dep=0;dep>(MAXLOG-(dep+1)))&1; if(k) r[q++]=data[i],mat[dep].set(i); else l[p++]=data[i]; } zs[dep]=p; mat[dep].build(); swap(l,data); for(Int i=0;i>(MAXLOG-(dep+1)))&1; l=mat[dep].rank(bit,l)+zs[dep]*bit; r=mat[dep].rank(bit,r)+zs[dep]*bit; } return r-l; } // return the position of k-th v Int select(T v,Int k){ rank(v,len); for(Int dep=MAXLOG-1;dep>=0;dep--){ bool bit=(v>>(MAXLOG-(dep+1)))&1; k=mat[dep].select(bit,k,buff1[dep]); if(k>=buff2[dep]||k<0) return -1; k-=buff1[dep]; } return k; } Int select(T v,Int k,Int l){ return select(v,k+rank(v,l)); } // return k-th largest value in [l,r) T quantile(Int l,Int r,Int k){ if(r-l<=k||k<0) return -1; T res=0; for(Int dep=0;depk){ l=p+zs[dep]; r=q+zs[dep]; res|=(1<<(MAXLOG-(dep+1))); }else{ k-=(q-p); l-=p; r-=q; } } return res; } T rquantile(Int l,Int r,Int k){ return quantile(l,r,r-l-k-1); } // return number of points in [left, right) * [lower, upper) Int rangefreq(Int left,Int right,T lower,T upper){ return freq_dfs(0,left,right,0,lower,upper); } pair ll(Int l,Int r,T v){ Int res=0; for(Int dep=0;dep>(MAXLOG-(dep+1)))&1; if(bit) res+=r-l+mat[dep].rank(bit,l)-mat[dep].rank(bit,r); l=mat[dep].rank(bit,l)+zs[dep]*bit; r=mat[dep].rank(bit,r)+zs[dep]*bit; } return make_pair(res,r-l); } Int lt(Int l,Int r,T v){ auto p=ll(l,r,v); return p.first; } Int le(Int l,Int r,T v){ auto p=ll(l,r,v); return p.first+p.second; } T succ(Int l,Int r,T v){ Int k=le(l,r,v); return k==r-l?npos:rquantile(l,r,k); } T pred(Int l,Int r,T v){ Int k=lt(l,r,v); return k?rquantile(l,r,k-1):npos; } }; //INSERT ABOVE HERE signed main(){ Int n; cin>>n; vector x(n); for(Int i=0;i>x[i],x[i]--; WaveletMatrix wm(x); Int ans=n*(n-1); for(Int i=0;i