#include #include #include #include #include #define _USE_MATH_DEFINES #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() const int INF = 0x3f3f3f3f; const long long LINF = 0x3f3f3f3f3f3f3f3fLL; const double EPS = 1e-8; const int MOD = 1000000007; // const int MOD = 998244353; const int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1}; // const int dy[] = {1, 1, 0, -1, -1, -1, 0, 1}, // dx[] = {0, -1, -1, -1, 0, 1, 1, 1}; struct IOSetup { IOSetup() { cin.tie(nullptr); ios_base::sync_with_stdio(false); cout << fixed << setprecision(20); cerr << fixed << setprecision(10); } } iosetup; /*-------------------------------------------------*/ template struct CHT { using Line = pair; CHT(bool is_minimized = true) : is_minimized(is_minimized) {} void add(T a, T b) { Line now(a, b); if (deq.empty()) { deq.emplace_back(now); } else if (deq.back().first <= a) { if (deq.back().first == a) { if (is_minimized) { if (b >= deq.back().second) return; deq.pop_back(); } else { if (b <= deq.back().second) return; deq.pop_back(); } } while (deq.size() >= 2 && !is_necessary(deq[deq.size() - 2], deq.back(), now)) deq.pop_back(); deq.emplace_back(now); } else { if (deq.front().first == a) { if (is_minimized) { if (b >= deq.front().second) return; deq.pop_front(); } else { if (b <= deq.front().second) return; deq.pop_front(); } } while (deq.size() >= 2 && !is_necessary(now, deq[0], deq[1])) deq.pop_front(); deq.emplace_front(now); } } T query(T x) { int lb = -1, ub = deq.size() - 1; while (ub - lb > 1) { int mid = (lb + ub) >> 1; if (is_minimized) { (f(deq[mid], x) < f(deq[mid + 1], x) ? ub : lb) = mid; } else { (f(deq[mid], x) > f(deq[mid + 1], x) ? ub : lb) = mid; } } return f(deq[ub], x); } T monotone_inc_query(T x) { if (is_minimized) { while (deq.size() >= 2 && f(deq[deq.size() - 2], x) <= f(deq.back(), x)) deq.pop_back(); return f(deq.back(), x); } else { while (deq.size() >= 2 && f(deq[0], x) <= f(deq[1], x)) deq.pop_front(); return f(deq.front(), x); } } T monotone_dec_query(T x) { if (is_minimized) { while (deq.size() >= 2 && f(deq[0], x) >= f(deq[1], x)) deq.pop_front(); return f(deq.front(), x); } else { while (deq.size() >= 2 && f(deq[deq.size() - 2], x) >= f(deq.back(), x)) deq.pop_back(); return f(deq.back(), x); } } private: bool is_minimized; deque deq; using Real = long double; bool is_necessary(const Line &l1, const Line &l2, const Line &l3) { Real lhs = static_cast(l3.second - l2.second) / (l2.first - l3.first); Real rhs = static_cast(l2.second - l1.second) / (l1.first - l2.first); // T lhs = (l1.first - l2.first) * (l3.second - l2.second); // T rhs = (l2.first - l3.first) * (l2.second - l1.second); return is_minimized ? lhs < rhs : lhs > rhs; } T f(const Line &l, T x) { return l.first * x + l.second; } }; int main() { int n; cin >> n; vector a(n + 1, 0); FOR(i, 1, n + 1) cin >> a[i]; FOR(i, 2, n + 1) a[i] += a[i - 1]; vector > dp(n + 2, vector(n + 1, -LINF)); dp[0][0] = 0; FOR(i, 1, n + 2) dp[i][i - 1] = (a[i - 1] - a[0]) * (a[i - 1] - a[0]); FOR(i, 1, n + 1) { CHT cht; for (int k = 0; i + k <= n; ++k) { cht.add(-a[i + k] * 2, dp[i + k][k] + a[i + k] * a[i + k]); dp[i + k + 1][k] = cht.monotone_inc_query(a[i + k]) + a[i + k] * a[i + k]; } } FOR(j, 1, n + 1) cout << dp[n + 1][j] << '\n'; return 0; }