This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include<bits/stdc++.h>
using namespace std;
typedef long long lint;
typedef long double ldb;
typedef unsigned long long uli;
#define X first
#define Y second
#define F(i, l, r) for(int i = l; i != r; i++)
#define Df(i, l, r) for(int i = l; i != r; i--)
#define I(i, a) for(auto i : a)
#define pb push_back
#define rs resize
#define mk make_pair
#define asg assign
#define all(x) x.begin(),x.end()
#define ret return
#define cont continue
#define brk break
#define ins insert
#define era erase
#define fi0(x) memset(x, 0, sizeof(x))
#define finf(x) memset(x, 127, sizeof(x))
#define acpy(y, x) memcpy(y, x, sizeof(y))
#define y1 adjf
#define tm dhgdg
const int MAXN = 1e3 + 3;
const int MAXK = 2e2 + 2;
lint dp[MAXN][MAXK];
int pr[MAXN][MAXK];
int a[MAXN];
int pref[MAXN];
int n, k;
struct line{
lint k, b;
int id;
line(lint pk, lint pb, int pid):k(pk),b(pb),id(pid){}
lint f(int x){
ret k * x + b;
}
};
struct convex_hull{
vector<line> l;
int ptr = 0;
void add_line(line curr){
while(l.size() >= 2){
line l1 = l[(int)l.size() - 1], l2 = l[(int)l.size() - 2];
if(l1.k == curr.k && l1.b <= curr.b){l.pop_back(); cont;}
lint f1 = (l1.k - l2.k) * (l1.b - curr.b), f2 = (l2.b - l1.b) * (curr.k - l1.k);
if(f2 < f1)brk;
l.pop_back();
}
l.pb(curr);
ptr = min(ptr, (int)l.size() - 1);
}
int get(int x){
if(l.empty())ret -1;
while(ptr < (int)l.size() - 1 && l[ptr].f(x) <= l[ptr + 1].f(x))ptr++;
ret l[ptr].id;
}
};
convex_hull ch[MAXK];
void restore(int i, int j){
if(i < 0)ret;
printf("%d ", i + 1);
restore(pr[i][j], j - 1);
}
int main(){
ios_base::sync_with_stdio(0);
cin.tie(0);
scanf("%d%d", &n, &k);
F(i, 0, n)scanf("%d", &a[i]);
F(i, 0, n)pref[i] = a[i] + ((i == 0) ? 0 : pref[i - 1]);
k++;
F(i, 0, n)
F(j, 1, k + 1)
pr[i][j] = -2;
dp[0][1] = 0;
pr[0][1] = -1;
ch[1].add_line(line(pref[0], -pref[0] * 1ll * pref[0], 0));
F(i, 1, n){
dp[i][1] = 0;
pr[i][1] = -1;
ch[1].add_line(line(pref[i], -pref[i] * 1ll * pref[i], i));
F(j, 2, k + 1){
int pi = ch[j - 1].get(pref[i]);
if(pi == -1)cont;
pr[i][j] = pi;
dp[i][j] = dp[pi][j - 1] + pref[pi] * 1ll * (pref[i] - pref[pi]);
ch[j].add_line(line(pref[i], dp[i][j] - pref[i] * 1ll * pref[i], i));
}
}
printf("%lld\n", dp[n - 1][k]);
restore(pr[n - 1][k], k - 1);
ret 0;
}
Compilation message (stderr)
sequence.cpp: In function 'int main()':
sequence.cpp:79:23: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d%d", &n, &k);
^
sequence.cpp:80:30: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
F(i, 0, n)scanf("%d", &a[i]);
^
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |