#ifndef LOCAL #define FAST_IO #endif // ===== template.hpp ===== #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define OVERRIDE(a, b, c, d, ...) d #define REP2(i, n) for (i32 i = 0; i < (i32) (n); ++i) #define REP3(i, m, n) for (i32 i = (i32) (m); i < (i32) (n); ++i) #define REP(...) OVERRIDE(__VA_ARGS__, REP3, REP2)(__VA_ARGS__) #define PER(i, n) for (i32 i = (i32) (n) - 1; i >= 0; --i) #define ALL(x) begin(x), end(x) using namespace std; using u32 = unsigned int; using u64 = unsigned long long; using u128 = __uint128_t; using i32 = signed int; using i64 = signed long long; using i128 = __int128_t; using f64 = double; using f80 = long double; template using Vec = vector; template bool chmin(T &x, const T &y) { if (x > y) { x = y; return true; } return false; } template bool chmax(T &x, const T &y) { if (x < y) { x = y; return true; } return false; } istream &operator>>(istream &is, i128 &x) { i64 v; is >> v; x = v; return is; } ostream &operator<<(ostream &os, i128 x) { os << (i64) x; return os; } istream &operator>>(istream &is, u128 &x) { u64 v; is >> v; x = v; return is; } ostream &operator<<(ostream &os, u128 x) { os << (u64) x; return os; } template > Vec sort_index(i32 n, F f, Comp comp = Comp()) { Vec idx(n); iota(ALL(idx), 0); sort(ALL(idx), [&](i32 i, i32 j) -> bool { return comp(f(i), f(j)); }); return idx; } [[maybe_unused]] constexpr i32 INF = 1000000100; [[maybe_unused]] constexpr i64 INF64 = 3000000000000000100; #ifdef FAST_IO __attribute__((constructor)) void fast_io() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(10); } #endif // ===== template.hpp ===== #ifdef DEBUGF #include "cpl/template/debug.hpp" #else #define DBG(x) (void) 0 #endif int main() { i32 n; i64 m; i64 l; cin >> n >> m >> l; Vec a(n); REP(i, n) { cin >> a[i]; } i32 lg = 0; while (n * m > (1LL << lg)) { ++lg; } i32 step = l % n; Vec steps(lg + 1); steps[0] = step; REP(i, lg) { steps[i + 1] = 2 * steps[i] % n; } Vec> gain(lg + 1, Vec(n)); gain[0] = a; REP(i, lg) { REP(j, n) { gain[i + 1][j] = gain[i][j] + gain[i][(j + steps[i]) % n]; } } DBG(gain); const auto calc_cx = [&](i64 x) -> i64 { i64 cnt = n * m / l; if (x < n * m % l) { ++cnt; } DBG(x); DBG(cnt); i64 sum = 0; i32 cur = x % n; REP(i, lg + 1) { if ((cnt >> i) & 1) { sum += gain[i][cur]; cur = (cur + steps[i]) % n; } } DBG(sum); return sum; }; i64 ans = -INF64; REP(i, n) { chmax(ans, calc_cx(i % l)); } if (l > n) { for (i64 i = l - n; i < l; ++i) { chmax(ans, calc_cx(i)); } } for (i64 i = n * m - n; i < n * m; ++i) { chmax(ans, calc_cx(i % l)); } cout << ans << '\n'; }