제출 #365168

#제출 시각아이디문제언어결과실행 시간메모리
365168YJU수열 (APIO14_sequence)C++14
100 / 100
998 ms83832 KiB
#include<bits/stdc++.h>
#pragma GCC optimize("unroll-loops,no-stack-protector,Ofast")
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll,ll> pll;
const ll MOD=1e9+7;
const ll MOD2=998244353;
const ll N=1e5+5;
const ll K=2e2+5;
const ld pi=acos(-1);
const ll INF=(1LL<<60);
#define SQ(i) ((i)*(i))
#define REP(i,n) for(ll i=0;i<n;i++)
#define REP1(i,n) for(ll i=1;i<=n;i++)
#define pb push_back
#define mp make_pair
#define X first
#define Y second
#define setp setprecision
#define lwb lower_bound
#define SZ(_a) (ll)_a.size()

ll dp[2][N],n,k,u[N];
int from[K][N];
vector<ll> ans;

deque<pair<pll,ll> > dq;

int main(){
	ios_base::sync_with_stdio(0);cin.tie(0);
	cin>>n>>k;
	REP1(i,n)cin>>u[i],u[i]+=u[i-1];
	REP1(i,n)dp[0][i]=-INF;
	dp[0][0]=0;
	REP1(T,k+1){
		dq.clear();
		dq.pb(mp(mp(0,dp[(T-1)&1][0]),0));
		REP1(i,n){
			while(SZ(dq)>1&&dq[0].X.X*u[i]+dq[0].X.Y<=dq[1].X.X*u[i]+dq[1].X.Y)dq.pop_front();
			
			dp[T&1][i]=dq[0].X.X*u[i]+dq[0].X.Y;
			from[T][i]=dq[0].Y;
			
			ll a=u[i],b=dp[(T-1)&1][i]-u[i]*u[i];
			while(SZ(dq)>1){
				ll m=SZ(dq)-1;
				if((b-dq[m].X.Y)*(dq[m-1].X.X-dq[m].X.X)<=(dq[m].X.Y-dq[m-1].X.Y)*(dq[m].X.X-a))dq.pop_back();
				else break;
			}
			dq.pb(mp(mp(a,b),i));
		}
	}
	cout<<dp[(k+1)&1][n]<<"\n";
	for(int T=k+1,i=n;T>=1&&from[T][i];T--){
		ans.pb(i=from[T][i]);
	}
	reverse(ans.begin(),ans.end());
	REP(i,k)cout<<ans[i]<<(i==k-1?"\n":" ");
	return 0;
}
/*
7 3
4 1 3 4 0 2 3
*/
#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...