제출 #208829

#제출 시각아이디문제언어결과실행 시간메모리
208829baohiep수열 (APIO14_sequence)C++14
0 / 100
65 ms83320 KiB
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 1e5 + 5, K = 205;

struct Line {
  ll a, b;
  int id;
  Line(ll _a, ll _b, int _id): a(_a), b(_b), id(_id) {}

  ll calc(ll x) {
    return a * (x - a) + b;
  }
};

int n, k, a[N], trace[N][K];
ll sum[N], pre[N], dp[N];
vector<Line> dq;

bool overlap(const Line &p1, const Line &p2, const Line &p3) {
  return (p2.a - p1.a) * (p1.b - p3.b) <= (p3.a - p1.a) * (p1.b - p2.b);
}

int main() {
  ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
  cin >> n >> k;
  for (int i = 1; i <= n; ++i) {
    cin >> a[i];
    sum[i] = sum[i - 1] + a[i];
  }

  for (int i = 1; i <= k; ++i) {
    int seen = 0;
    dq.clear();
    dq.emplace_back(0, 0, 0);
    for (int j = 1; j <= n; ++j) {
      while (seen + 1 < (int)dq.size() && dq[seen].calc(sum[j]) <= dq[seen + 1].calc(sum[j]))
        ++seen;
      dp[j] = dq[seen].calc(sum[j]);
      trace[j][i] = dq[seen].id;
      Line newLine = Line(sum[j], pre[j], j);
      while (dq.size() > 2 && overlap(dq[(int)dq.size() - 2], dq.back(), newLine))
        dq.pop_back();
      dq.push_back(newLine);
    }
    for (int j = 1; j <= n; ++j) pre[j] = dp[j];
  }
  cout << dp[n] << "\n";
  int cur = n;
  while (k) {
    cur = trace[cur][k--];
    cout << cur << " ";
  }
  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...