#include using namespace std; using lint = long long int; using pint = pair; using plint = pair; struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_; #define ALL(x) (x).begin(), (x).end() #define SZ(x) ((lint)(x).size()) #define POW2(n) (1LL << (n)) #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i=i##_begin_;i--) #define REP(i, n) FOR(i,0,n) #define IREP(i, n) IFOR(i,0,n) template void ndarray(vector &vec, int len) { vec.resize(len); } template void ndarray(vector &vec, int len, Args... args) { vec.resize(len); for (auto &v : vec) ndarray(v, args...); } template bool mmax(T &m, const T q) { if (m < q) {m = q; return true;} else return false; } template bool mmin(T &m, const T q) { if (m > q) {m = q; return true;} else return false; } template pair operator+(const pair &l, const pair &r) { return make_pair(l.first + r.first, l.second + r.second); } template pair operator-(const pair &l, const pair &r) { return make_pair(l.first - r.first, l.second - r.second); } template istream &operator>>(istream &is, vector &vec){ for (auto &v : vec) is >> v; return is; } ///// This part below is only for debug, not used ///// template ostream &operator<<(ostream &os, const vector &vec){ os << "["; for (auto v : vec) os << v << ","; os << "]"; return os; } template ostream &operator<<(ostream &os, const deque &vec){ os << "deq["; for (auto v : vec) os << v << ","; os << "]"; return os; } template ostream &operator<<(ostream &os, const set &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template ostream &operator<<(ostream &os, const unordered_set &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template ostream &operator<<(ostream &os, const multiset &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template ostream &operator<<(ostream &os, const unordered_multiset &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template ostream &operator<<(ostream &os, const pair &pa){ os << "(" << pa.first << "," << pa.second << ")"; return os; } template ostream &operator<<(ostream &os, const map &mp){ os << "{"; for (auto v : mp) os << v.first << "=>" << v.second << ","; os << "}"; return os; } template ostream &operator<<(ostream &os, const unordered_map &mp){ os << "{"; for (auto v : mp) os << v.first << "=>" << v.second << ","; os << "}"; return os; } #define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl; ///// END ///// template struct LazySegTree { int N; int head; vector val; vector dval; VAL I_val; DVAL I_dval; using vv2v = function; using d2d = function; using d2v = function; vv2v merge_val; d2d add_dval; d2v refl_dval; LazySegTree(vector val_init, VAL val_default, DVAL dval_default, vv2v f, d2d dadd, d2v dreflect) : N(val_init.size()), I_val(val_default), I_dval(dval_default), merge_val(f), add_dval(dadd), refl_dval(dreflect) { int N_tmp = 1; while (N_tmp < N) N_tmp <<= 1; val = vector(N_tmp*2, I_val); dval = vector(N_tmp*2, I_dval); head = N_tmp - 1; for (int i=0; i=0; pos--) val[pos] = merge_val(val[pos*2+1], val[pos*2+2]); } void resolve_dval(int pos, int l, int r) { // posで遅延を解消して子孫に押し付ける refl_dval(val[pos], dval[pos]); if (pos < head) { add_dval(dval[pos*2+1], dval[pos]); add_dval(dval[pos*2+2], dval[pos]); } dval[pos] = I_dval; } vector> history; void update(int begin, int end, DVAL dval_q) { update(begin, end, dval_q, 0, 0, head+1); history.emplace_back(pint(begin, end), dval_q); } void roll_back() { for (auto p : history) { update(p.first.first, p.first.second, -p.second, 0, 0, head + 1); } history.clear(); } void update(int begin, int end, DVAL dval_q, int pos, int l, int r) { // 後でリファクタリング if (begin <= l && r <= end) { // 担当区間全部使う add_dval(dval[pos], dval_q); resolve_dval(pos, l, r); } else if (begin < r && l < end) { // 少なくともどこかで交差 resolve_dval(pos, l, r); update(begin, end, dval_q, pos*2+1, l, (l+r)/2); update(begin, end, dval_q, pos*2+2, (l+r)/2, r); val[pos] = merge_val(val[pos*2+1], val[pos*2+2]); } else resolve_dval(pos, l, r); } VAL getvalue(int begin, int end) { return getvalue(begin, end, 0, 0, head+1); } VAL getvalue(int begin, int end, int pos, int l, int r) { resolve_dval(pos, l, r); if (begin <= l && r <= end) return val[pos]; else if (begin> N; vector A(N); cin >> A; map> B; REP(i, N) B[A[i]].emplace_back(i); lint ret = 0; LazySegTree st(vector(N * 2 + 10), 0, 0, [](lint a, lint b) { return a + b; }, // 本情報のmerge [](lint &a, lint b) { a += b; }, // 遅延情報の子孫への伝播 [](lint &a, lint b) { a += b; }); // 遅延情報の本情報への反映 for (auto p : B) { p.second.push_back(N); int now = -1; int v = N + 5; st.update(v, v + 1, 1); for (auto x : p.second) { while (now + 1 < x) { now++; v--; st.update(v, v + 1, 1); lint addv = st.getvalue(0, v); ret += addv; if (!addv) break; } int down = x - now - 1; if (down > 0) { st.update(v - down, v, 1); now += down; v -= down; } if (x != N) { v++; st.update(v, v + 1, 1); ret += st.getvalue(0, v); } now = x; } st.roll_back(); } cout << ret << endl; }