BZOJ:4209: 西瓜王

原题链接：http://www.lydsy.com/JudgeOnline/problem.php?id=4209

（虽然仅仅是看在名字的份上，我们还是得说这题是一道）

绝世好题！西瓜王！西瓜王！西瓜王！[wàng]

说点正经的，一开始以为合成圣堂的时候必须相邻的，然后变成了区间内取k个数不能相邻的数的最大总价值，诶我的瓜妈妈呀，不会做啊……

看一下CA爷blog：贪心，然后替换……

诶，我看错题了？

最后就是主席树维护一下区间前k大的总价值，再查一下这最大的k个里面有多少个奇数多少偶数，如果奇数数量是奇数的话（好拗口）那就再查一下在前k大的数中最小的奇数和偶数，不在前k大中最大奇数和偶数，然后替换一下就吼了啊。。。

卡常卡不过啊T_T

#include<cstdio>
#include<algorithm>
#define MN 300001
#define ll long long
using namespace std;

char ibuf[],*iptr=ibuf-,C;
inline int _(){
    int x=,c=*++iptr;
    while(c<)c=*++iptr;
    while(c>)x=x*+c-,c=*++iptr;
    return x;
}
char obuf[],*optr=obuf;
inline void _(ll x){
    static int stk[],stp=;
    if (x<) *(optr++)='-',x=-x;
    if(!x)stk[stp++]=;
    while(x)stk[stp++]=x%,x/=;
    while(stp)*(optr++)=stk[--stp]+;
    *(optr++)=;
}
struct tree{int l,r,s[],S;ll al;}t[MN*];
int n,m,a[MN],num=,l,r,k,s[],S[MN][],ro[MN],mi,ma;
ll mmh,c;
inline ll min(ll a,ll b){return a<b?a:b;}
inline ll max(ll a,ll b){return a>b?a:b;}
void in(int &P,int x,int v){
    int l=mi,r=ma,mid,*p;p=&P;
    for (;;){
        *p=++num;t[*p]=t[x];t[*p].s[v&]++;t[*p].S++;t[*p].al+=v;
        if (l==r) break;
        int mid=l+r>>;
        if (v<=mid) p=&t[*p].l,x=t[x].l,r=mid;else p=&t[*p].r,x=t[x].r,l=mid+;
    }
}
void ask(int p,int x,int l,int r,int k){
    while (l<r){
        if (k>=t[t[p].r].S-t[t[x].r].S) s[]+=t[t[p].r].s[]-t[t[x].r].s[],s[]+=t[t[p].r].s[]-t[t[x].r].s[],
        mmh+=t[t[p].r].al-t[t[x].r].al,k-=t[t[p].r].S-t[t[x].r].S,p=t[p].l,x=t[x].l,r=l+r>>;else p=t[p].r,x=t[x].r,l=(l+r>>)+;
        if (k==) return;
    }
    if (k) mmh+=1LL*k*l,s[l&]+=k;
}
ll askmax(int p,int x,int k,bool bo){
    if (t[p].s[bo]-t[x].s[bo]<=k) return -1e18;
    int l=mi,r=ma;
    while (l<r)
    if (t[t[p].r].s[bo]-t[t[x].r].s[bo]>k) p=t[p].r,x=t[x].r,l=(l+r>>)+;else k-=t[t[p].r].s[bo]-t[t[x].r].s[bo],p=t[p].l,x=t[x].l,r=l+r>>;
    return l;
}
ll askmin(int p,int x,int k,bool bo){
    if (t[p].s[bo]-t[x].s[bo]<k) return 1e18;
    int l=mi,r=ma;
    while (l<r)
    if (t[t[p].r].s[bo]-t[t[x].r].s[bo]>=k) p=t[p].r,x=t[x].r,l=(l+r>>)+;else k-=t[t[p].r].s[bo]-t[t[x].r].s[bo],p=t[p].l,x=t[x].l,r=l+r>>;
    return l;
}
int main(){
    ibuf[fread(ibuf,,,stdin)]=;
    register int i;
    n=_();t[].l=t[].r=t[].s[]=t[].s[]=t[].al=t[].S=;
    for (i=;i<=n;i++) a[i]=_(),S[i][]=S[i-][],S[i][]=S[i-][],S[i][a[i]&]++,mi=min(mi,a[i]),ma=max(ma,a[i]);
    for (i=;i<=n;i++) in(ro[i],ro[i-],a[i]);
    m=_();
    while (m--){
        l=_();r=_();k=_();
        if ((S[r][]-S[l-][]>>)+(S[r][]-S[l-][]>>)<(k>>)){_(-);continue;}
        mmh=s[]=s[]=;ask(ro[r],ro[l-],mi,ma,k);
        if (s[]&) mmh-=min(askmin(ro[r],ro[l-],s[],)-askmax(ro[r],ro[l-],s[],),askmin(ro[r],ro[l-],s[],)-askmax(ro[r],ro[l-],s[],));
        _(mmh);
    }
    fwrite(obuf,,optr-obuf,stdout);
    return ;
}