#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> 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; } } dp[0][0]=0; rep(i,N){ Rep(j,1,K+1){ deque> DQ1,DQ2; auto key1=dp[i][j-1]+sum[i]-sum[i+1]; auto key2=dp[i][j-1]-sum[i]+sum[i+1]; while(1){ if(DQ1.size() && !SLcomp(DQ1.back().fr,key1)) DQ1.pop_back(); else break; } while(1){ if(DQ2.size() && !SLcomp(DQ2.back().fr,key2)) DQ2.pop_back(); else break; } DQ1.emplace_back(mp(key1,j)); DQ2.emplace_back(mp(key2,j)); if(DQ1.front().sc==i-M) DQ1.pop_front(); if(DQ2.front().sc==i-M) DQ2.pop_front(); dp[i+1][j]=max(dp[i+1][j],DQ1.front().fr); dp[i+1][j]=max(dp[i+1][j],DQ2.front().fr); Rep(k,max(0ll,i-M+1),i+1){ dp[i+1][j]=max(dp[i+1][j],dp[k][j-1]+sum[k]-sum[i+1]); dp[i+1][j]=max(dp[i+1][j],dp[k][j-1]-sum[k]+sum[i+1]); } } } cout << dp[N][K] << endl; }