// >>> TEMPLATES #include using namespace std; #define rep(i,n) for (int i = 0; i < (int)(n); i++) #define rep1(i,n) for (int i = 1; i <= (int)(n); i++) #define repR(i,n) for (int i = (int)(n)-1; i >= 0; i--) #define rep1R(i,n) for (int i = (int)(n); i >= 1; i--) #define loop(i,a,B) for (int i = a; i B; i++) #define loopR(i,a,B) for (int i = a; i B; i--) #define all(x) (x).begin(), (x).end() #define allR(x) (x).rbegin(), (x).rend() #define pb push_back #define eb emplace_back #define mp make_pair #define fst first #define snd second #ifdef LOCAL #define dump(...) cerr << "[" << __LINE__ << ":" << __FUNCTION__ << "] ", dump_impl(#__VA_ARGS__, __VA_ARGS__) #define say(x) cerr << "[" << __LINE__ << ":" << __FUNCTION__ << "] " << x << endl; #define debug if (1) void dump_impl(const char*) { cerr << endl; } template void dump_impl(const char *s, T const& x, U const& ...y) { const char *o = "({[", *e = "]})"; for (int i = 0; *s != '\0'; cerr << *s++) { if (count(o,o+3,*s)) i++; if (count(e,e+3,*s)) i--; if (!i && *s == ',') break; } cerr << " = " << x; if (*s == ',') cerr << ", ", s++; dump_impl(s, y...); } #else #define dump(...) #define say(x) #define debug if (0) #endif using ll = long long; using ld = long double; template using pque_max = priority_queue; template using pque_min = priority_queue, greater >; template ::value>::type> ostream& operator<<(ostream& os, T const& v) { os << "{"; for (auto const& x : v) os << " " << x; return os << " }"; } template ::value>::type> istream& operator>>(istream& is, T &v) { for (auto& x : v) is >> x; return is; } template ostream& operator<<(ostream& os, pair const& p) { return os << "(" << p.first << ", " << p.second << ")"; } template istream& operator>>(istream& is, pair& p) { return is >> p.first >> p.second; } template typename enable_if= tuple_size::value>::type output_tuple(ostream&, T const&) { } template typename enable_if::value>::type output_tuple(ostream& os, T const& t) { os << (i ? " " : "") << get(t); output_tuple(os,t); } template ostream& operator<<(ostream& os, tuple const& t) { return output_tuple(os,t), os; } struct IOSetup { IOSetup() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); } } iosetup; template struct FixPoint : private F { constexpr FixPoint(F&& f) : F(forward(f)) {} template constexpr auto operator()(T&&... x) const { return F::operator()(*this, forward(x)...); } }; struct MakeFixPoint { template constexpr auto operator|(F&& f) const { return FixPoint(forward(f)); } }; #define MFP MakeFixPoint()| #define def(name, ...) auto name = MFP [&](auto &&name, __VA_ARGS__) template struct vec_impl { using type = vector::type>; template static type make_v(size_t n, U&&... x) { return type(n, vec_impl::make_v(forward(x)...)); } }; template struct vec_impl { using type = T; static type make_v(T const& x = {}) { return x; } }; template using vec = typename vec_impl::type; template auto make_v(Args&&... args) { return vec_impl::make_v(forward(args)...); } template void quit(T const& x) { cout << x << endl; exit(0); } template constexpr bool chmin(T& x, T const& y) { if (x > y) { x = y; return true; } return false; } template constexpr bool chmax(T& x, T const& y) { if (x < y) { x = y; return true; } return false; } template constexpr auto sumof(It b, It e) { return accumulate(b,e,typename iterator_traits::value_type{}); } #define int ll #define double ld template int sz(T const& x) { return x.size(); } const ll INF = (1LL<<62)-1; // ~ 4.6e18 // <<< int32_t main() { int n; cin >> n; vector a(n); cin >> a; vector<__int128> s(n+1); rep (i,n) s[i+1] = s[i] + a[i]; auto dp = make_v(n+1,n+1,INF); rep (i,n+1) dp[0][i] = s[i]*s[i]; rep (j,n) { deque > q; auto add = [&](int a, int b) -> void { while (1) { if (sz(q) < 2) { q.eb(a,b); return; } auto q1 = q.back(); auto q0 = *prev(q.end(),2); int x1 = q1.fst - q0.fst, y1 = q1.snd - q0.snd; int x2 = a - q0.fst, y2 = b - q0.snd; if (x1*y2-y1*x2 >= 0) q.pop_back(); else { q.eb(a,b); return; } } }; auto query = [&](int x) -> int { while (1) { assert(sz(q)); if (sz(q) == 1) return q[0].fst * x + q[0].snd; if (q[0].fst * x + q[0].snd <= q[1].fst * x + q[1].snd) { q.pop_front(); } else { return q[0].fst * x + q[0].snd; } } }; rep (i,n) { add(2*s[i+1], -(dp[j][i]+s[i+1]*s[i+1])); dp[j+1][i+1] = s[i+1]*s[i+1] - query(s[i+1]); } } rep1 (k,n) cout << dp[n-k][n] << "\n"; }