제출 #7653

#제출 시각아이디문제언어결과실행 시간메모리
7653imsifile수열 (APIO14_sequence)C++98
33 / 100
8 ms83356 KiB
#include<stdio.h>
#include<deque>
#define sq(a) ((a)*(a))
#define sik(l,x) (-(l.a)*(x)+(l.b))
using namespace std;

typedef long long lld;

struct lin {
	int ix;
	lld a, b; // -ax+b
}im;

deque<lin> dq;

int n, k, di[100100][202];
lld ba[100100], dy[100100][2], sum[100100];

int main(){
	int i, j, t=1;
	scanf("%d%d", &n, &k);
	for(i=1; i<=n; i++){
		scanf("%lld", &ba[i]);
		sum[i]=sum[i-1]+ba[i], dy[i][0]=sq(sum[i]);
	}
	for(j=1; j<=k; j++){
		im.a=2*sum[j], im.b=dy[j][1-t]+sq(sum[j]), im.ix=j;
		dq.push_back(im);
		for(i=1; i<=j; i++)dy[i][t]=0;
		for(i=j+1; i<=n; i++){
			for(; dq.size()>=2;){
				if(sik(dq[0],sum[i])<sik(dq[1],sum[i]))break;
				dq.pop_front();
			}
			di[i][j]=dq[0].ix, dy[i][t]=sq(sum[i])+sik(dq[0],sum[i]);
			im.ix=i, im.a=2*sum[i], im.b=dy[i][1-t]+sq(sum[i]);
			for(; dq.size();){
				if(dq[dq.size()-1].b<im.b)break;
				dq.pop_back();
			}
			dq.push_back(im);
		}
		dq.clear();
		t=1-t;
	}
	printf("%lld\n", (sq(sum[n])-dy[n][1-t])/2);
	for(i=n, j=k; j>0; j--)printf("%d ", di[i][j]), i=di[i][j];
	return 0;
}
#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...