//#define _GLIBCXX_DEBUG #include using namespace std; #define endl '\n' #define lfs cout<= (ll)(n); i--) using ll = long long; using ld = long double; const ll MOD1 = 1e9+7; const ll MOD9 = 998244353; const ll INF = 1e18; using P = pair; templatebool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;} templatebool chmax(T1 &a,T2 b){if(avoid ans(bool x,T1 y,T2 z){if(x)cout<void debug(vector>&v,ll h,ll w){for(ll i=0;i&v,ll h,ll w){for(ll i=0;ivoid debug(vector&v,ll n){if(n!=0)cout<vector>vec(ll x, ll y, T w){vector>v(x,vector(y,w));return v;} ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;} vectordx={1,-1,0,0,1,1,-1,-1};vectordy={0,0,1,-1,1,-1,1,-1}; templatevector make_v(size_t a,T b){return vector(a,b);} templateauto make_v(size_t a,Ts... ts){return vector(a,make_v(ts...));} templateostream &operator<<(ostream &os, const pair&p){return os << p.first << " " << p.second;} templateostream &operator<<(ostream &os, const vector &v){for(auto &z:v)os << z << " ";cout<<"|"; return os;} //mt19937 mt(chrono::steady_clock::now().time_since_epoch().count()); int popcount(ll x){return __builtin_popcountll(x);}; int poplow(ll x){return __builtin_ctzll(x);}; int pophigh(ll x){return 63 - __builtin_clzll(x);}; template< typename T, bool isMin > struct ConvexHullTrickAddMonotone { #define F first #define S second using P = pair< T, T >; deque< P > H; ConvexHullTrickAddMonotone() = default; bool empty() const { return H.empty(); } void clear() { H.clear(); } inline ll sgn(T x) { return x == 0 ? 0 : (x < 0 ? -1 : 1); } using D = long double; inline bool check(const P &a, const P &b, const P &c) { if(b.S == a.S || c.S == b.S) return sgn(b.F - a.F) * sgn(c.S - b.S) >= sgn(c.F - b.F) * sgn(b.S - a.S); //return (b.F-a.F)*(c.S-b.S) >= (b.S-a.S)*(c.F-b.F); return D(b.F - a.F) * sgn(c.S - b.S) / D(abs(b.S - a.S)) >= D(c.F - b.F) * sgn(b.S - a.S) / D(abs(c.S - b.S)); } void add(T a, T b) { if(!isMin) a *= -1, b *= -1; P line(a, b); if(empty()) { H.emplace_front(line); return; } if(H.front().F <= a) { if(H.front().F == a) { if(H.front().S <= b) return; H.pop_front(); } while(H.size() >= 2 && check(line, H.front(), H[1])) H.pop_front(); H.emplace_front(line); } else { assert(a <= H.back().F); if(H.back().F == a) { if(H.back().S <= b) return; H.pop_back(); } while(H.size() >= 2 && check(H[H.size() - 2], H.back(), line)) H.pop_back(); H.emplace_back(line); } } inline T get_y(const P &a, const T &x) { return a.F * x + a.S; } T query(T x) { assert(!empty()); ll l = -1, r = H.size() - 1; while(l + 1 < r) { ll m = (l + r) >> 1; if(get_y(H[m], x) >= get_y(H[m + 1], x)) l = m; else r = m; } if(isMin) return get_y(H[r], x); return -get_y(H[r], x); } T query_monotone_inc(T x) { assert(!empty()); while(H.size() >= 2 && get_y(H.front(), x) >= get_y(H[1], x)) H.pop_front(); if(isMin) return get_y(H.front(), x); return -get_y(H.front(), x); } T query_monotone_dec(T x) { assert(!empty()); while(H.size() >= 2 && get_y(H.back(), x) >= get_y(H[H.size() - 2], x)) H.pop_back(); if(isMin) return get_y(H.back(), x); return -get_y(H.back(), x); } #undef F #undef S }; int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); ll res=0,buf=0; bool judge = true; ll n;cin>>n; vectora(n); rep(i,0,n)cin>>a[i]; vector>cht(n+1); vectorb(n+1); rep(i,0,n)b[i+1]=b[i]+a[i]; cht[0].add(-2*b[0],b[0]*b[0]); auto dp=make_v(n+1,n+1,INF); dp[0][0]=0; rep(i,0,n){ rep(j,0,i+1){ ll mi=cht[j].query_monotone_inc(b[i+1])+b[i+1]*b[i+1]; //このタイミングでi-j+1回処理を行っている chmin(dp[i+1][i-j+1],mi); chmin(dp[i+1][j],dp[i][j]); } rep(j,0,i+1){ cht[i-j+1].add(-2*b[i+1],b[i+1]*b[i+1]+dp[i][j]); } } //debug(dp,n+1,n+1); vectorret(n,INF); rep(i,0,n+1)rep(j,1,n+1)chmin(ret[j-1],dp[i][j]); rep(i,0,n)cout<