#pragma region library #if defined(__GNUC__) && !defined(__llvm__) && !defined(__clang__) #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #pragma GCC target("avx2") #include #endif #if __has_include() #include #endif #pragma region #define whole(f, x, ...) ([&](auto& v) { return (f)(std::begin(v), std::end(v), ## __VA_ARGS__); })(x) #define rep32(i, l, r) for (int i = (int)(l); (i) < (int)(r); ++(i)) #define rep64(i, l, r) for (long long i = (long long)(l); i < (long long)(r); ++(i)) #define revrep32(i, r, l) for (int i = (int)(r); (i) >= (int)(l); --(i)) #define revrep64(i, r, l) for (long long i = (long long)(r); i >= (long long)(l); --(i)) using u32 = unsigned int; using usize = size_t; using index_t = size_t; using i64 = long long; using u64 = unsigned long long; constexpr size_t operator "" _uz(unsigned long long v) { return static_cast(v); } template using Heap = std::priority_queue; template using RevHeap = std::priority_queue, std::greater>; template using SHeap = std::priority_queue, F>; constexpr int dx[] = { 1, 0, -1, 0 }; constexpr int dy[] = { 0, 1, 0, -1 }; constexpr int INF32 = 1001001001; constexpr long long INF64 = 1501501501501501501ll; constexpr int mod = 1000000007; constexpr int mod2 = 998244353; #pragma endregion template constexpr bool chmin(T& a, const U& b) noexcept { return a > b ? a = b, true : false; } template constexpr bool chmax(T& a, const U& b) noexcept { return a < b ? a = b, true : false; } template constexpr type len(const T& v) noexcept { return (type)(v.size()); } template std::vector> zip(std::vector& t, std::vector&... ts) { const size_t n = t.size(); std::vector> res; for (size_t i = 0; i < n; i++) res.emplace_back(t[i], ts[i]...); return res; } template std::vector> idzip(std::vector& t, std::vector&... ts) { const size_t n = t.size(); std::vector> res; for (size_t i = 0; i < n; i++) res.emplace_back(i, t[i], ts[i]...); return res; } namespace nlb { template constexpr T Pow(T A, U B) noexcept { T res = 1; while (B) { if (B & 1) { res *= A; } A *= A; B >>= 1; } return res; } template constexpr T modpow(T A, U B, M MOD = 1000000007) noexcept { A %= MOD; T res = 1; while (B) { if (B & 1) { res *= A; res %= MOD; } A *= A; A %= MOD; B >>= 1; } return res; } template constexpr T inverse(T A, const M MOD = 1000000007) noexcept { T B = MOD, U = 1, V = 0; while (B) { T t = A / B; A -= t * B; std::swap(A, B); U -= t * V; std::swap(U, V); } U %= MOD; return U < 0 ? U += MOD, U : U; } template constexpr bool isprime(const T N) noexcept { if (N == 1) return false; for (T i = 2; i * i <= N; i++) if (N % i == 0) return false; return true; } inline std::tuple extgcd(const long long a, const long long b) { if (b == 0) return { a,1,0 }; long long g, x, y; std::tie(g, x, y) = extgcd(b, a % b); return std::make_tuple(g, y, x - a / b * y); } inline std::pair Chinese_Rem(const std::vector& b, const std::vector& m) { long long r = 0, M = 1; for (int i = 0; i < (int)b.size(); ++i) { long long p, q, d; std::tie(d, p, q) = extgcd(M, m[i]); if ((b[i] - r) % d != 0) return std::make_pair(0, -1); long long tmp = (b[i] - r) / d * p % (m[i] / d); r += M * tmp; M *= m[i] / d; } return std::make_pair((r % M + M) % M, M); } } class ProCon_all_init { static constexpr bool float_fixed = true; static constexpr int ios_perc = 15; static constexpr bool ios_fast = true; static constexpr bool auto_flush = false; public: ProCon_all_init() { if (ios_fast) { std::cin.tie(nullptr); std::cout.tie(nullptr); std::ios::sync_with_stdio(false); } if (float_fixed) { std::cout << std::fixed << std::setprecision(ios_perc); } if (auto_flush) { std::cout << std::unitbuf; } } } ProCon_; #pragma endregion int main(void) { // input i64 N; std::cin >> N; std::vector A(N); for (auto& el : A) std::cin >> el; // solve // 分割統治法をします // 中心を含む or 含まない // l が中心にちかづくするにつれて r は離れる std::function&)> solver = [&](std::vector v)->i64 { int n = len(v); if (n < 2) return 0; if (n == 2) return 1; int half = n / 2; int half2 = n - half; std::vector lvec, rvec; rep32(i, 0, half) lvec.emplace_back(v[i]); rep32(i, half, n) rvec.emplace_back(v[i]); std::multiset S; for (auto& el : lvec) S.emplace(el); i64 R = 0; i64 Answer = 0; rep32(i, 0, half) { // 尺取り while (len(S) < 2 && R < half2) { S.emplace(rvec[R++]); } while (R < half2) { chmax(R, 0); S.emplace(rvec[R]); auto itr1 = S.begin(); auto itr2 = itr1; itr2++; auto itr3 = S.end(); itr3--; if (!(*itr1 + *itr2 >= *itr3)) { S.erase(S.find(rvec[R])); break; } R++; } Answer += R; S.erase(S.find(lvec[i])); } return Answer + solver(lvec) + solver(rvec); }; std::cout << solver(A) << std::endl; }