Submission #732634

#TimeUsernameProblemLanguageResultExecution timeMemory
732634ac2huSplit the sequence (APIO14_sequence)C++14
100 / 100
793 ms86584 KiB
#include<bits/stdc++.h> using namespace std; using ll = long long; const int N = 1e5 + 10; const int K = 200 + 10; int n, k, b[N][K]; ll dp1[N], dp2[N], a[N]; struct line { ll m, c;int id; ll eval(ll x) { return (ll)(m * x) + c; } long double intersectX(line l) { return (long double) (c - l.c) / (l.m - m); } }; signed main(){ iostream::sync_with_stdio(false);cin.tie(nullptr); cin >> n >> k; for(int i = 1;i<=n;i++)cin >> a[i]; for(int i = 1;i<=n;i++)a[i] += a[i - 1]; for(int l = 1;l<=k;l++){ deque<line> dq;dq.push_front({0, 0, 0}); for(int i = 1;i<=n;i++){ if(i > l){ while(dq.size() >= 2 && dq.back().eval(a[i]) <= dq[dq.size() -2].eval(a[i]))dq.pop_back(); dp2[i] = dq.back().eval(a[i]); b[i][l] = dq.back().id; } if(i > l - 1){ line cur = {a[i], -a[i]*a[i] + dp1[i], i}; while(dq.size() >= 2 && cur.intersectX(dq[0]) <= cur.intersectX(dq[1])) dq.pop_front(); dq.push_front(cur); } } for(int i = 0;i<=n;i++)dp1[i] = dp2[i]; } cout << dp1[n] << "\n"; vector<int> v = {n}; for(int l = k;l>=1;l--){ v.push_back(b[v.back()][l]); } reverse(v.begin(), v.end());v.pop_back();reverse(v.begin(), v.end()); for(int e : v)cout << e << " "; cout << "\n"; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...