#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using pii = pair; using ll=long long; using ld=long double; #define pb push_back #define mp make_pair #define sc second #define fr first #define stpr setprecision #define cYES cout<<"YES"<=0;i--) #define rRep(i,a,b) for(ll i=a;i>=b;i--) #define crep(i) for(char i='a';i<='z';++i) #define psortsecond(A,N) sort(A,A+N,[](const pii &a, const pii &b){return a.second inline bool chmax(T& lhs, const U& rhs) { if (lhs < rhs) { lhs = rhs; return 1; } return 0; } template inline bool chmin(T& lhs, const U& rhs) { if (lhs > rhs) { lhs = rhs; return 1; } return 0; } template istream& operator>>(istream& is,vector& v){for(auto&& x:v)is >> x;return is;} template istream& operator>>(istream& is, pair& p){ is >> p.first; is >> p.second; return is;} template ostream& operator>>(ostream& os, const pair& p){ os << p.first << ' ' << p.second; return os;} template ostream& operator<<(ostream& os, vector& v){ for(auto i=begin(v); i != end(v); ++i){ if(i !=begin(v)) os << ' '; os << *i; } return os; } // modint template class modint { using u64 = std::uint_fast64_t; public: u64 a; constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {} constexpr u64 &value() noexcept { return a; } constexpr const u64 &value() const noexcept { return a; } constexpr modint operator+(const modint rhs) const noexcept { return modint(*this) += rhs; } constexpr modint operator-(const modint rhs) const noexcept { return modint(*this) -= rhs; } constexpr modint operator*(const modint rhs) const noexcept { return modint(*this) *= rhs; } constexpr modint operator/(const modint rhs) const noexcept { return modint(*this) /= rhs; } constexpr modint &operator+=(const modint rhs) noexcept { a += rhs.a; if (a >= Modulus) { a -= Modulus; } return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (a < rhs.a) { a += Modulus; } a -= rhs.a; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { a = a * rhs.a % Modulus; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = Modulus - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; ll A[3007],sum[3007],dp[3007][3007]; bool SLcomp(ll a,ll b){ return a>b; } int main(){ ll N,K,M; cin >> N >> K >> M; rep(i,N){ cin >> A[i]; } rep(i,N){ sum[i+1]=sum[i]+A[i]; } rep(i,N+2){ rep(j,K+2){ dp[i][j]=-LINF; } } vector>> dq1(K + 1), dq2(K + 1); dq1[0].emplace_back(0, 0); dq2[0].emplace_back(0, 0); dp[0][0] = 0; for (int i = 0; i < N; i++) { for (int j = 1; j <= K; j++) { ll mx1 = -INF; ll mx2 = -INF; while (!dq1[j - 1].empty() and dq1[j - 1].front().first < i + 1 - M) dq1[j - 1].pop_front(); while (!dq2[j - 1].empty() and dq2[j - 1].front().first < i + 1 - M) dq2[j - 1].pop_front(); if (!dq1[j - 1].empty()) mx1 = dq1[j - 1].front().second; if (!dq2[j - 1].empty()) mx2 = dq2[j - 1].front().second; chmax(dp[i + 1][j], mx1 - sum[i + 1]); chmax(dp[i + 1][j], mx2 + sum[i + 1]); } for (int j = 0; j <= K; j++) { while (!dq1[j].empty() and dq1[j].back().second < dp[i + 1][j] + sum[i + 1]) dq1[j].pop_back(); while (!dq2[j].empty() and dq2[j].back().second < dp[i + 1][j] - sum[i + 1]) dq2[j].pop_back(); dq1[j].emplace_back(i + 1, dp[i + 1][j] + sum[i + 1]); dq2[j].emplace_back(i + 1, dp[i + 1][j] - sum[i + 1]); } } cout << dp[N][K] << endl; }