#include #include #include #include using usize = std::size_t; template using Vec = std::vector; Vec> compute_inverse(const usize n, const Vec p) { Vec> sum(n + 1, Vec(n + 1, 0)); for (usize i = 0; i < n; ++i) { ++sum[i + 1][p[i] + 1]; } for (usize i = 0; i < n; ++i) { for (usize j = 0; j <= n; ++j) { sum[i + 1][j] += sum[i][j]; } } for (usize i = 0; i <= n; ++i) { for (usize j = 0; j < n; ++j) { sum[i][j + 1] += sum[i][j]; } } Vec> inv(n, Vec(n + 1, 0)); for (usize i = 0; i < n; ++i) { usize s = 0; for (usize j = i; j < n; ++j) { s += sum[j][n] - sum[j][p[j]] - sum[i][n] + sum[i][p[j]]; inv[i][j + 1] = s; } } return inv; } class Manager { // f(x) = inv[i][x] + offset (i < x) // max range: [?, limit] struct Func { usize i; usize offset; usize limit; Func(usize _i, usize _o, usize _l) : i(_i), offset(_o), limit(_l) {} }; const usize n; const Vec> inv; Vec funcs; public: Manager(usize _n, Vec> _i) : n(_n), inv(std::move(_i)), funcs() {} void add(const usize i, const usize offset) { if (funcs.empty()) { funcs.emplace_back(i, offset, n); return; } while (funcs.back().limit <= i) { funcs.pop_back(); } if (const Func &f = funcs.back(); offset <= inv[f.i][i + 1] + f.offset) { return; } while (!funcs.empty()) { const Func &f = funcs.back(); if (inv[f.i][f.limit] + f.offset < inv[i][f.limit] + offset) { funcs.pop_back(); } else { usize ok = i + 1; usize ng = f.limit; while (ng - ok > 1) { const usize mid = (ok + ng) / 2; if (inv[f.i][mid] + f.offset < inv[i][mid] + offset) { ok = mid; } else { ng = mid; } } funcs.emplace_back(i, offset, ok); return; } } funcs.emplace_back(i, offset, n); } usize max(const usize x) { while (funcs.back().limit < x) { funcs.pop_back(); } const Func &f = funcs.back(); return inv[f.i][x] + f.offset; } void clear() { funcs.clear(); } }; Vec solve(const usize n, Vec p) { const Vec> inv = compute_inverse(n, std::move(p)); Manager manager(n, inv); Vec> dp(n + 1, Vec(n + 1)); for (usize i = 1; i <= n; ++i) { dp[1][i] = inv[0][i]; } for (usize k = 2; k <= n; ++k) { manager.add(k - 1, dp[k - 1][k - 1]); for (usize i = k; i <= n; ++i) { dp[k][i] = manager.max(i); if (i != n) { manager.add(i, dp[k - 1][i]); } } manager.clear(); } Vec ans(n + 1, 0); for (usize i = 1; i <= n; ++i) { ans[i] = inv[0][n] - dp[i][n]; } return ans; } int main() { usize n; std::cin >> n; Vec p(n); for (usize &ele : p) { std::cin >> ele; --ele; } const Vec ans = solve(n, std::move(p)); for (usize i = 1; i <= n; ++i) { std::cout << ans[i] << '\n'; } }