| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 1000719 | shenfe1 | 수열 (APIO14_sequence) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
#define ll long long
#define lb lower_bound
#define pii pair<int,int>
#define pll pair<ll,ll>
#define F first
#define S second
#define ld long double
#define pb push_back
#define all(v) v.begin(),v.end()
#define in insert
#define sz(s) (int)s.size()
// #define int ll
#define ppb pop_back
using namespace std;
const int MAX=1e5+10;
const int inf=2e9;
int n,k;
int a[MAX];
vector<int> dp[2];
int pr[MAX][211];
ll p[MAX];
struct CHT{
vector<pll> vec;
vector<int> num;
int pos;
void init(){
pos=0;
vec.clear();
num.clear();
}
ld inter(pll a,pll b){
// cout<<b.F<<" "<<a.F<<"\n";
return (0.0+a.S-b.S)/(b.F-a.F);
}
void add(ll k,ll x,int i){
// cout<<k<<" "<<x<<"\n";
while(!vec.empty()&&vec.back().F==k){
if(vec.back().S<x){
vec.ppb();
num.ppb();
}
else{
return;
}
}
while(sz(vec)>=2&&inter(vec.back(),vec[sz(vec)-2])>=inter(vec[sz(vec)-2],{k,x})){
vec.ppb();
num.ppb();
}
vec.pb({k,x});
num.pb(i);
pos=min(pos,sz(vec)-1);
}
ll get(pll a,int x){
return a.F*x+a.S;
}
ll get(ll X){
while(pos+1<sz(vec)&&get(vec[pos],X)<=get(vec[pos+1],X))pos++;
return get(vec[pos],X);
}
}C;
void solve(){
cin>>n>>k;
for(int i=1;i<=n;i++)cin>>a[i];
for(int i=1;i<=n;i++){
p[i]=p[i-1]+a[i];
}
k++;
dp[0]=vector<ll>(MAX,0);
dp[1]=vector<ll>(MAX,0);
for(int K=2;K<=k;K++){
C.init();
for(int i=K;i<=n;i++){
C.add(p[i-1],dp[0][i-1]-p[i-1]*p[i-1],i-1);
dp[1][i]=C.get(p[i]);
pr[i][K]=C.num[C.pos];
// cout<<i<<" "<<K<<" "<<pr[i][K]<<"\n";
}
swap(dp[0],dp[1]);
}
cout<<dp[0][n]<<"\n";
vector<int> ans;
int cur=n;
while(k>0){
ans.pb(cur);
cur=pr[cur][k];
k--;
}
sort(all(ans));
ans.ppb();
for(int x:ans){
cout<<x;
if(x!=ans.back())cout<<" ";
}
}
// #ifdef LOCAL
signed main(){
ios_base::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
int t=1;
// cin>>t;
while(t--)solve();
}
// #endif