#include // created [2019/12/15] 00:08:40 constexpr unsigned int mod = 1000000007; template constexpr T inf_v = std::numeric_limits::max() / 4; template constexpr Real pi_v = Real{3.141592653589793238462643383279502884}; using i32 = int32_t; using i64 = int64_t; using u32 = uint32_t; using u64 = uint64_t; using uint = unsigned int; using usize = std::size_t; using ll = long long; using ull = unsigned long long; using ld = long double; template using arr = T (&)[n]; template using c_arr = const T (&)[n]; template struct simple_cht { private: using line = std::pair; static bool check(line l1, line l2, line l3) { const element_type a1 = l1.first, b1 = l1.second, a2 = l2.first, b2 = l2.second, a3 = l3.first, b3 = l3.second; const value_type da1 = a1 - a2, da2 = a2 - a3, db1 = b1 - b2, db2 = b2 - b3; return (da1 + da2) * db1 > da1 * (db1 + db2) and da2 * (db1 + db2) > (da1 + da2) * db2; } static Val eval(const line& l, const Elem x) { return static_cast(l.first) * static_cast(x) + static_cast(l.second); } usize query_head = 0; Elem px = -inf_v; std::vector lines; public: using element_type = Elem; using value_type = Val; void add(const element_type a, const element_type b) { if (lines.empty()) { return lines.push_back(line{a, b}); } const element_type pa = lines.back().first, pb = lines.back().second; assert(pa >= a); if (pa == a and pb <= b) { return; } if (lines.size() == 1) { return lines.push_back(line{a, b}); } const line l{a, b}; while (lines.size() >= 2 and not check(lines[lines.size() - 2], lines[lines.size() - 1], l)) { lines.pop_back(); } lines.push_back(l), query_head = std::min(query_head, lines.size() - 1); } value_type query(const element_type x) { assert(px <= x), px = x; if (lines.empty()) { return inf_v; } for (; query_head + 1 < lines.size() and eval(lines[query_head], x) > eval(lines[query_head + 1], x); query_head++) {} return eval(lines[query_head], x); } bool empty() const { return lines.empty(); } friend std::ostream& operator<<(std::ostream& os, const simple_cht& cht) { os << "["; for (usize i = 0; i < cht.lines.size(); i++) { os << cht.lines[i].first << "x+" << cht.lines[i].second << (i + 1 == cht.lines.size() ? "" : ","); } return (os << "]\n"); } }; #pragma GCC diagnostic ignored "-Wsign-compare" #pragma GCC diagnostic ignored "-Wsign-conversion" template constexpr T popcount(const T u) { return u ? static_cast(__builtin_popcountll(static_cast(u))) : static_cast(0); } template constexpr T log2p1(const T u) { return u ? static_cast(64 - __builtin_clzll(static_cast(u))) : static_cast(0); } template constexpr T msbp1(const T u) { return log2p1(u); } template constexpr T lsbp1(const T u) { return __builtin_ffsll(u); } template constexpr T clog(const T u) { return u ? log2p1(u - 1) : static_cast(u); } template constexpr bool ispow2(const T u) { return u and (static_cast(u) & static_cast(u - 1)) == 0; } template constexpr T ceil2(const T u) { return static_cast(1) << clog(u); } template constexpr T floor2(const T u) { return u == 0 ? static_cast(0) : static_cast(1) << (log2p1(u) - 1); } template constexpr bool btest(const T mask, const usize ind) { return static_cast((static_cast(mask) >> ind) & static_cast(1)); } template void bset(T& mask, const usize ind) { mask |= (static_cast(1) << ind); } template void breset(T& mask, const usize ind) { mask &= ~(static_cast(1) << ind); } template void bflip(T& mask, const usize ind) { mask ^= (static_cast(1) << ind); } template void bset(T& mask, const usize ind, const bool b) { (b ? bset(mask, ind) : breset(mask, ind)); } template constexpr T bcut(const T mask, const usize ind) { return ind == 0 ? static_cast(0) : static_cast((static_cast(mask) << (64 - ind)) >> (64 - ind)); } template bool chmin(T& a, const T& b) { return (a > b ? a = b, true : false); } template bool chmax(T& a, const T& b) { return (a < b ? a = b, true : false); } auto fix_point = [](auto f) { return [=](auto&&... args) { return f(f, std::forward(args)...); }; }; template T in() { T v; return std::cin >> v, v; } template T in_v(typename std::enable_if<(i == n), c_arr>::type) { return in(); } template auto in_v(typename std::enable_if<(i < n), c_arr>::type& szs) { const usize s = (usize)szs[i]; std::vector(szs))> ans(s); for (usize j = 0; j < s; j++) { ans[j] = in_v(szs); } return ans; } template auto in_v(c_arr szs) { return in_v(szs); } template auto in_t() { return std::tuple...>{in()...}; } struct io_init { io_init() { std::cin.tie(nullptr), std::ios::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(20); } } hogechan; template void out(const T& v) { std::cout << v; } template void out(const T& v, const Args... args) { std::cout << v << ' ', out(args...); } template void outln(const T& v, const Args... args) { std::cout << v << ' ', out(args...), std::cout << '\n'; } template void outel(const T& v, const Args... args) { std::cout << v << ' ', out(args...), std::cout << std::endl; } # define SHOW(...) static_cast(0) constexpr ull TEN(const usize n) { return n == 0 ? 1ULL : TEN(n - 1) * 10ULL; } template auto make_v(typename std::enable_if<(i == n), c_arr>::type, const T& v = T{}) { return v; } template auto make_v(typename std::enable_if<(i < n), c_arr>::type szs, const T& v = T{}) { const usize s = (usize)szs[i]; return std::vector(szs, v))>(s, make_v(szs, v)); } template auto make_v(c_arr szs, const T& t = T{}) { return make_v(szs, t); } int main() { const auto N = in(); std::vector S(N + 2, 0); for (int i = 1; i <= N; i++) { S[i] = in(); } for (int i = 1; i <= N + 1; i++) { S[i] += S[i - 1]; } std::vector> chts(N + 2); // i個閉めるときのCHT auto dp = make_v({N + 2, N + 2}, inf_v); // ドアiを閉めて/合計j個のドアを閉める dp[0][0] = 0; chts[0].add(-2LL * S[0], S[0] * S[0] + dp[0][0]); for (int i = 1; i <= N + 1; i++) { for (int j = 1; j <= N + 1; j++) { if (chts[j - 1].empty()) { continue; } chmin(dp[i][j], S[i - 1] * S[i - 1] + (ll)chts[j - 1].query(S[i - 1])); chts[j].add(-2LL * S[i], dp[i][j] + S[i] * S[i]); } } for (int j = N; j >= 1; j--) { std::cout << dp[N + 1][j] << std::endl; } return 0; }