#include <bits/stdc++.h>
#include "holiday.h"
#define ll long long
#define pb push_back
#define mp make_pair
using namespace std;
const int maxn = 150100;
ll value[maxn], d;
ll arr[maxn];
class node {
public:
ll sum, cnt, lb, rb;
node *leftchild, *rightchild;
node(ll _sum, ll _cnt) {
sum = _sum;
cnt = _cnt;
}
node(node *a, node *b) {
leftchild = a;
rightchild = b;
sum = a->sum + b->sum;
cnt = a->cnt + b->cnt;
}
};
node* build(int li=1, int ri=d) {
if(li == ri) return new node(0LL, 0LL);
else {
int mid = (li + ri) / 2;
return new node(build(li, mid), build(mid+1, ri));
}
}
node* update(node *curr, ll uind, ll uval, int li=1, int ri=d) {
if(li == ri) {
return new node(curr->sum + uval, curr->cnt + 1);
}
else {
int mid = (li + ri) / 2;
if(uind <= mid) return new node(update(curr->leftchild, uind, uval, li, mid), curr->rightchild);
else return new node(curr->leftchild, update(curr->rightchild, uind, uval, mid+1, ri));
}
}
ll solve(node *p, node *pm, ll k, int li=1, int ri=d) {
ll total = p->cnt - pm->cnt;
//cout<<li<<" "<<ri<<" -> "<<total<<", "<<k<<"\n";
if(k == 0) return 0LL;
if(total < k) return (p->sum - pm->sum);
if(li == ri) {
return k * value[li];
}
else {
int mid = (li + ri) / 2;
ll total_right = p->rightchild->cnt - pm->rightchild->cnt;
if(total_right >= k) return solve(p->rightchild, pm->rightchild, k, mid+1, ri);
else return (p->rightchild->sum - pm->rightchild->sum) + solve(p->leftchild, pm->leftchild, k - total_right, li, mid);
}
}
node *prefix[maxn];
set<int>st;
map<int,int>ind;
ll n, start, c;
ll attraction[maxn];
ll dp[maxn];
void solvedp(int l=1, int r=start, int optl=start, int optr=n) {
if(l > r) return;
int mid = (l + r) / 2;
// calculate dp[mid]
dp[mid] = -1;
int optmid = -1;
for(int i=optl;i<=optr;i++) {
ll walk = i - mid + min(start - mid, i - start);
ll curr = solve(prefix[i], prefix[mid-1], max(0LL, c - walk));
if(curr > dp[mid]) {
dp[mid] = curr;
optmid = i;
}
}
solvedp(l, mid-1, optl, optmid);
solvedp(mid+1, r, optmid, optr);
}
long long int findMaxAttraction(int _n, int _start, int _c, int _attraction[]) {
//ifstream fin("input.txt");
n = _n;
start = _start;
c = _c;
start++;
for(int i=0;i<n;i++) {
attraction[i] = _attraction[i];
st.insert(attraction[i]);
}
int br = 1;
for(int i:st) {
value[br] = i;
ind[i] = br++;
}
d = br - 1;
for(int i=n-1;i>=0;i--) {
arr[i+1] = ind[attraction[i]];
}
node *curr = build();
prefix[0] = curr;
for(int i=1;i<=n;i++) {
curr = update(curr, arr[i], value[arr[i]]);
prefix[i] = curr;
}
solvedp();
ll result = 0LL;
for(int i=1;i<=n;i++) {
result = max(result, dp[i]);
}
return result;
}
