/** * @FileName a.cpp * @Author kanpurin * @Created 2020.07.19 21:32:15 **/ #include "bits/stdc++.h" using namespace std; typedef long long ll; template < class Monoid > struct SegmentTree { private: using Func = std::function< Monoid(Monoid, Monoid) >; Func F; Monoid UNITY; int n; std::vector< Monoid > node; public: SegmentTree() {} SegmentTree(int m, const Monoid& val, const Func f, const Monoid& unity) { F = f; UNITY = unity; n = 1; while (n < m) n <<= 1; node.resize(n * 2 - 1, UNITY); if (val != UNITY) { for (int i = 0; i < m; i++) node[i] = val; build(); } } SegmentTree(const std::vector< Monoid >& v, const Func f, const Monoid& unity) { F = f; UNITY = unity; int sz = v.size(); n = 1; while (n < sz) n <<= 1; node.resize(n * 2 - 1, UNITY); for (int i = 0; i < sz; i++) node[i + n - 1] = v[i]; build(); } void set(int k, const Monoid& x) { node[n + k - 1] = x; } void build() { for (int i = n - 2; i >= 0; i--) node[i] = F(node[2 * i + 1], node[2 * i + 2]); } void update_query(int x, const Monoid& val) { if (x >= n || x < 0) return; x += n - 1; node[x] = val; while (x > 0) { x = (x - 1) >> 1; node[x] = F(node[2 * x + 1], node[2 * x + 2]); } } Monoid get_query(int a, int b) { Monoid L = UNITY, R = UNITY; if (a < 0) a = 0; if (b < 0) return UNITY; if (a >= n) return UNITY; if (b >= n) b = n; for (a += n, b += n; a < b; a >>= 1, b >>= 1) { if (a & 1) L = F(L, node[a++ - 1]); if (b & 1) R = F(node[--b - 1], R); } return F(L, R); } Monoid operator[](int x) const { return node[n + x - 1]; } int size() { return n; } void print() { for (int i = 0; i < n; i++) { std::cout << i << "\t: " << node[n + i - 1] << std::endl; } } }; int main() { int n, k, m; cin >> n >> k >> m; vector< int > a(n); vector< ll > sum(n + 1); for (int i = 0; i < n; i++) { cin >> a[i]; sum[i + 1] = sum[i] + a[i]; } long long LLINF = 1e18 + 1; vector< vector< ll > > dp(k + 1, vector< ll >(n + 1, -LLINF)); dp[0][0] = 0; SegmentTree< ll > seg1(n + 1, -LLINF, [](ll a, ll b) { return max(a, b); }, -LLINF); SegmentTree< ll > seg2(n + 1, -LLINF, [](ll a, ll b) { return max(a, b); }, -LLINF); seg1.update_query(0, 0); seg2.update_query(0, 0); for (int j = 0; j < k; j++) { for (int i = n - 1; i >= 0; i--) { dp[j + 1][i + 1] = max(dp[j + 1][i + 1], seg1.get_query(i + 1 - m, i + 1) + sum[i + 1]); dp[j + 1][i + 1] = max(dp[j + 1][i + 1], seg2.get_query(i + 1 - m, i + 1) - sum[i + 1]); seg1.update_query(i + 1, max(seg1[i + 1],dp[j + 1][i + 1] - sum[i + 1])); seg2.update_query(i + 1, max(seg2[i + 1],dp[j + 1][i + 1] + sum[i + 1])); } } cout << dp[k][n] << endl; return 0; }