#include using namespace std; // Define using ll = long long; using ull = unsigned long long; using ld = long double; template using pvector = vector>; template using rpriority_queue = priority_queue, greater>; constexpr const ll dx[4] = {1, 0, -1, 0}; constexpr const ll dy[4] = {0, 1, 0, -1}; constexpr const ll MOD = 1e9 + 7; constexpr const ll mod = 998244353; constexpr const ll INF = 1LL << 60; constexpr const ll inf = 1 << 30; constexpr const char rt = '\n'; constexpr const char sp = ' '; #define rt(i, n) (i == (ll)(n) -1 ? rt : sp) #define len(x) ((ll)(x).size()) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() #define mp make_pair #define mt make_tuple #define pb push_back #define eb emplace_back #define ifn(x) if (not(x)) #define elif else if #define elifn else ifn #define fi first #define se second template bool chmax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; } template bool chmin(T &a, const T &b) { if (a > b) { a = b; return 1; } return 0; } // Debug #define debug(...) \ { \ cerr << __LINE__ << ": " << #__VA_ARGS__ << " = "; \ for (auto &&X : {__VA_ARGS__}) cerr << "[" << X << "] "; \ cerr << rt; \ } #define dump(a, h, w) \ { \ cerr << __LINE__ << ": " << #a << " = [" << rt; \ rep(_i, h) { \ rep(_j, w) { \ if (abs(a[_i][_j]) >= INF / 2 and a[_i][_j] <= -INF / 2) \ cerr << '-'; \ if (abs(a[_i][_j]) >= INF / 2) \ cerr << "∞" << sp; \ else \ cerr << a[_i][_j] << sp; \ } \ cerr << rt; \ } \ cerr << "]" << rt; \ } #define vdump(a, n) \ { \ cerr << __LINE__ << ": " << #a << " = ["; \ rep(_i, n) { \ if (_i) cerr << sp; \ if (abs(a[_i]) >= INF / 2 and a[_i] <= -INF / 2) cerr << '-'; \ if (abs(a[_i]) >= INF / 2) \ cerr << "∞" << sp; \ else \ cerr << a[_i]; \ } \ cerr << "]" << rt; \ } // Loop #define inc(i, a, n) for (ll i = (a), _##i = (n); i <= _##i; ++i) #define dec(i, a, n) for (ll i = (a), _##i = (n); i >= _##i; --i) #define rep(i, n) for (ll i = 0, _##i = (n); i < _##i; ++i) #define each(i, a) for (auto &&i : a) // Stream #define fout(n) cout << fixed << setprecision(n) struct io { io() { cin.tie(nullptr), ios::sync_with_stdio(false); } } io; // Speed #pragma GCC optimize("Ofast,unroll-loops") #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native") // Math inline constexpr ll gcd(const ll a, const ll b) { return b ? gcd(b, a % b) : a; } inline constexpr ll lcm(const ll a, const ll b) { return a / gcd(a, b) * b; } inline constexpr ll modulo(const ll n, const ll m = MOD) { ll k = n % m; return k + m * (k < 0); } inline constexpr ll chmod(ll &n, const ll m = MOD) { n %= m; return n += m * (n < 0); } inline constexpr ll mpow(ll a, ll n, const ll m = MOD) { ll r = 1; rep(i, 64) { if (n & (1LL << i)) r *= a; chmod(r, m); a *= a; chmod(a, m); } return r; } inline ll inv(const ll n, const ll m = MOD) { ll a = n, b = m, x = 1, y = 0; while (b) { ll t = a / b; a -= t * b; swap(a, b); x -= t * y; swap(x, y); } return modulo(x, m); } struct UshiTree { using Type = ll; Type unit = -INF; function comb = [](Type left, Type right) { return max(left, right); }; static const ll n = 1 << 12; vector node; UshiTree() : node(2 * n - 1, unit) {} inline void update(ll a, Type x) { ll i = n - 1 + a; node[i] = x; while (i != 0) { i = (i - 1) / 2; node[i] = comb(node[i * 2 + 1], node[i * 2 + 2]); } } inline Type result(ll a, ll b, ll i = 0, ll l = 0, ll r = n) { Type res = unit; if (b <= l or r <= a) return res; if (a <= l and r <= b) { return comb(node[i], res); } else { return comb(result(a, b, i * 2 + 1, l, (l + r) / 2), result(a, b, i * 2 + 2, (l + r) / 2, r)); } } }; struct slidemax { ll N; vector data; ll l, r; deque deq; slidemax(ll n = 3001, ll start = 0) : N(n), data(n), l(start), r(start) {} void setdata(ll itr, ll val) { data[itr] = val; } void pop_front() { while (deq.size() && deq.front() <= l) { deq.pop_front(); } l++; } void push_back() { while (deq.size() && data[deq.back()] < data[r]) { deq.pop_back(); } deq.push_back(r); r++; } ll min() { return data[deq.front()]; } }; UshiTree seg[3001]; slidemax slm[3001]; signed main() { ll N, M, K; cin >> N >> K >> M; ll A[N + 1] = {}; pvector V; rep(i, N) cin >> A[i + 1]; rep(i, N) A[i + 1] += A[i]; rep(i, N + 1) V.pb({A[i], i}); sort(all(V)); ll id[N + 1]; rep(i, N + 1) id[V[i].se] = i; ll DP[N + 1][K + 1]; rep(i, N + 1) rep(j, K + 1) DP[i][j] = -INF; DP[0][0] = 0; seg[0].update(id[0], -A[0]); slm[0].setdata(0, A[0]); slm[0].pb(); rep(i, N) { rep(j, K) { if (i - M >= 0) { seg[j].update(id[i - M], -INF); slm[j].pop_front(); } DP[i + 1][j + 1] = max(seg[j].result(0, id[i + 1] + 1) + A[i + 1], slm[j].min() - A[i + 1]); seg[j + 1].update(id[i + 1], DP[i + 1][j + 1] - A[i + 1]); slm[j + 1].setdata(i + 1, DP[i + 1][j + 1] + A[i + 1]); slm[j + 1].pb(); } } cout << DP[N][K] << rt; }