#include typedef long long ll; typedef unsigned long long ull; #define FOR(i,a,b) for(int (i)=(a);i<(b);i++) #define REP(i,n) FOR(i,0,n) #define RANGE(vec) (vec).begin(),(vec).end() template class BIT { public: std::vector data; // [1,n] BIT() {} BIT(int n) { init(n); } void init(int n) { data.resize((1LL<<(int)ceil(log2(n))), 0); } // a[i] += x O(log(n)) void add(int i, int x) { int maxN = data.size()+1; // 2 冪 for (int k = i+1; k <= maxN; k += (k & -k)) data[k-1] += x; } // a[0]+...+a[i-1] O(log(n)) T sum(int i) { T s = 0; for (int k = i+1; k > 0; k -= (k & -k)) s += data[k-1]; return s; } }; using namespace std; // // 愚直なやり方だと O(N^3) で TLE // class TrendAndCountermeasures_PineDecorationSequence2 { public: void solve(void) { int N; cin>>N; map conv; // mapping 用 vector A(N); REP(i,N) { cin>>A[i]; conv[A[i]] = 0; } // A は一致しているかどうかが問題になるだけなので // 0...N-1 をマッピングして max(A) < N にする int k = 0; for (auto kv : conv) conv[kv.first] = k++; REP(i,N) A[i] = conv[A[i]]; vector L(N,0); vector R(N,0); BIT bit1(N); BIT bit2(N); REP(i,N) { ++R[A[i]]; bit2.add(A[i],1); } ll ret = 0; REP(i,N) { bit2.add(A[i],-1); ret += bit1.sum(A[i]-1)*bit2.sum(A[i]-1); ret += (i - bit1.sum(A[i])) * ((N-i-1) - bit2.sum(A[i])); bit1.add(A[i],1); } ll diff = 0; REP(i,N) { diff -= (ll)L[A[i]] * R[A[i]]; ret -= diff; ++L[A[i]]; --R[A[i]]; diff += (ll)L[A[i]] * R[A[i]]; } cout<solve(); delete obj; return 0; } #endif