// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/843 #include #include #include #include /** * @brief ダブリング * * @tparam L * @tparam Monoid モノイド */ template struct Doubling { private: using T = typename Monoid::value_type; public: template Doubling(const std::vector &to, const std::vector &v) : Doubling(to.size()) { build(to, v); } std::pair jump(int f, std::uint64_t k) { return solve(f, k); } std::pair solve(int f, std::uint64_t k) { assert(-1 <= f && f < _size); T res = Monoid::id(); for (int cnt = 0; k > 0; k >>= 1, ++cnt) { if ((k & 1) && f != -1) { res = Monoid::op(res, data[cnt][f]); f = table[cnt][f]; } } return std::make_pair(f, res); } private: int _size; std::vector> table; std::vector> data; Doubling(int n) : _size(n), table(L, std::vector(n)), data(L, std::vector(n)) {} template void build(const std::vector &to, const std::vector &v) { assert(to.size() == v.size()); for (int i = 0; i < _size; ++i) { assert(-1 <= to[i] && to[i] < _size); table[0][i] = to[i]; data[0][i] = v[i]; } for (int i = 0; i < L - 1; ++i) { for (int j = 0; j < _size; ++j) { int k = table[i][j]; if (k != -1) { table[i + 1][j] = table[i][k]; data[i + 1][j] = Monoid::op(data[i][j], data[i][k]); } else { table[i + 1][j] = table[i][j]; data[i + 1][j] = data[i][j]; } } } } }; /** * @brief ダブリング * * @tparam L */ template struct Doubling { Doubling(const std::vector &v) : Doubling(v.size()) { build(v); } int jump(int f, std::uint64_t k) { return solve(f, k); } int solve(int f, std::uint64_t k) { assert(-1 <= f && f < _size); for (int cnt = 0; k > 0; k >>= 1, ++cnt) { if ((k & 1) && f != -1) f = table[cnt][f]; } return f; } int lca(int u, int v) { for (int i = L - 1; i >= 0; --i) { if (table[u][i] != table[v][i]) u = table[u][i], v = table[v][i]; } return table[u][0]; } private: int _size; std::vector> table; Doubling(int n) : _size(n), table(L, std::vector(n)) {} void build(const std::vector &v) { for (int i = 0; i < _size; ++i) { assert(-1 <= v[i] && v[i] < _size); table[0][i] = v[i]; } for (int i = 0; i < L - 1; ++i) { for (int j = 0; j < _size; ++j) { if (table[i][j] != -1) table[i + 1][j] = table[i][table[i][j]]; } } } }; #ifdef ATCODER #pragma GCC target("sse4.2,avx512f,avx512dq,avx512ifma,avx512cd,avx512bw,avx512vl,bmi2") #endif #pragma GCC optimize("Ofast,fast-math,unroll-all-loops") #include #ifndef ATCODER #pragma GCC target("sse4.2,avx2,bmi2") #endif template constexpr bool chmax(T &a, const U &b) { return a < (T)b ? a = (T)b, true : false; } template constexpr bool chmin(T &a, const U &b) { return (T)b < a ? a = (T)b, true : false; } constexpr std::int64_t INF = 1000000000000000003; constexpr int Inf = 1000000003; constexpr double EPS = 1e-7; constexpr double PI = 3.14159265358979323846; #define FOR(i, m, n) for (int i = (m); i < int(n); ++i) #define FORR(i, m, n) for (int i = (m)-1; i >= int(n); --i) #define FORL(i, m, n) for (int64_t i = (m); i < int64_t(n); ++i) #define rep(i, n) FOR (i, 0, n) #define repn(i, n) FOR (i, 1, n + 1) #define repr(i, n) FORR (i, n, 0) #define repnr(i, n) FORR (i, n + 1, 1) #define all(s) (s).begin(), (s).end() struct Sonic { Sonic() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); std::cout << std::fixed << std::setprecision(20); } constexpr void operator()() const {} } sonic; using namespace std; using ll = std::int64_t; using ld = long double; template std::istream &operator>>(std::istream &is, std::pair &p) { return is >> p.first >> p.second; } template std::istream &operator>>(std::istream &is, std::vector &v) { for (T &i : v) is >> i; return is; } template std::ostream &operator<<(std::ostream &os, const std::pair &p) { return os << '(' << p.first << ',' << p.second << ')'; } template std::ostream &operator<<(std::ostream &os, const std::vector &v) { for (auto it = v.begin(); it != v.end(); ++it) os << (it == v.begin() ? "" : " ") << *it; return os; } template void co(Head &&head, Tail &&...tail) { if constexpr (sizeof...(tail) == 0) std::cout << head << '\n'; else std::cout << head << ' ', co(std::forward(tail)...); } template void ce(Head &&head, Tail &&...tail) { if constexpr (sizeof...(tail) == 0) std::cerr << head << '\n'; else std::cerr << head << ' ', ce(std::forward(tail)...); } void Yes(bool is_correct = true) { std::cout << (is_correct ? "Yes\n" : "No\n"); } void No(bool is_not_correct = true) { Yes(!is_not_correct); } void YES(bool is_correct = true) { std::cout << (is_correct ? "YES\n" : "NO\n"); } void NO(bool is_not_correct = true) { YES(!is_not_correct); } void Takahashi(bool is_correct = true) { std::cout << (is_correct ? "Takahashi" : "Aoki") << '\n'; } void Aoki(bool is_not_correct = true) { Takahashi(!is_not_correct); } struct Monoid { using value_type = ll; static constexpr ll id() { return 0; } static constexpr ll op(const ll &lhs, const ll &rhs) { return lhs + rhs; } }; int main(void) { ll n, m, l; cin >> n >> m >> l; vector a(n); cin >> a; vector to(n); rep (i, n) to[i] = (i + l) % n; Doubling<40, Monoid> db(to, a); ll r = (n * m) % l; ll k = (n * m) / l; ll ans = -INF; rep (i, min(r, n)) chmax(ans, db.solve(i, k + 1).second); rep (i, min(l - r, n)) chmax(ans, db.solve((r + i) % n, k).second); co(ans); return 0; }