省选九省联考T2 IIIDX(线段树)

题目传送门：https://www.luogu.org/problemnew/show/P4364

期中考后记：期中考刚考完，感觉不咋滴，年排第3。我抗压力太差了。。期末得把rank1抢回来。

本来感觉是个瞎贪心的（对每个歌曲的先决的歌曲为父子关系建树，然后每次找尽量大的孩子）。。其实不是。

贪心仅能骗60分。（省选题不会水的）

前面那个贪心只是瞎YY的。YY的贪心有一定几率出错。

为啥会出错，我一会儿讲。

我们正经得想一个贪心。其实问题变为：

for(i=1;i<=n;++i){

找一个难度为p的歌曲并删去，使在i后面的节点能依照没有p歌曲的歌单下不会没有难度足够大的歌曲的情况下，p尽量大，第i个歌曲难度就是p

}

其实这个贪心的思路非常好想（我最开始想的就是这个，后来去打暴力了，就用了第一个）

我们把设选了歌曲p后需要Ap个大于等于p的歌曲（包括p）。那么，我们只要计算出所有Ap，不就爽歪歪了吗？

最开始，我们将歌曲难度降序排序记为A，Ap=p。

写个伪代码

for(i=1;i<=n;++i){

　　找到最靠右的节点k使Ak-1<siz[i](i的子树大小)

　　ans[i]=k;

　　将ak+1.ak+2....an都减去siz[i];

}

for(i=1;i<=n;++i)输出d[ans[i]]

好了还记得最开始讲的贪心为啥会错吗？

因为，我刚才写的伪代码是错的（错的地方和第一个贪心一毛一样）

正确的在此

for(i=1;i<=n;++i){

　　找到最靠右的节点k使Ak-1<siz[i](i的子树大小)

　　如果相同的几个歌曲和k的难度一样，那么，将k移到在序列A中最靠左的那个；

　　ans[i]=k;

　　将ak+1.ak+2....an都减去siz[i];

}

为啥要那样干？反正没有后效性，而且贪心吗，就要贪得彻底，预留给自己的子树尽量多的节点。

代码：

#include <bits/stdc++.h>
using namespace std;
;
],dr[N<<],hd[N],pos,lovewxy[N];
double k;
bool cmp(int x,int y){return  x>y;}
bool vis[N];
],wxy[rt<<|]);}
void pushdown(int rt){
    if(dr[rt]){
        dr[rt<<]+=dr[rt];
        dr[rt<<|]+=dr[rt];
        wxy[rt<<]+=dr[rt];
        wxy[rt<<|]+=dr[rt];dr[rt]=;
    }
}
void build(int rt,int l,int r){
    if(l==r){
        wxy[rt]=l;return;
    }
    ;
    build(rt<<,l,mid);build(rt<<|,mid+,r);pushup(rt);
}
void motify(int rt,int l,int r,int ql,int qr,int d){
    if(ql<=l && r<=qr){
        wxy[rt]+=d;dr[rt]+=d;return;
    }pushdown(rt);
    ;
    ,l,mid,ql,qr,d);
    <=qr)motify(rt<<|,mid+,r,ql,qr,d);
    pushup(rt);
}
int query(int rt,int l,int r,int lst){
    ;
    pushdown(rt);;
    |]<lst)|,mid+,r,lst);
    ,l,mid,lst);
}
void dfs(int u){
    if(u>n)return;
    int l,r;double tmp=(double)u*k;
    l=tmp-;
    while(l<tmp)++l;
    tmp=()*k;
    r=tmp-;while(r<tmp)++r;
    --r;siz[u]=;
    for(int i=r;i>=l;--i)
    if(i<=n)
    dfs(i),siz[u]+=siz[i],fa[i]=u;
}
int main(){
    scanf("%d%lf",&n,&k);
    ;i<=n;++i)scanf(,d++n,cmp);d[n+]=-;
    ;--i)hd[i]= d[i]==d[i+] ? pos:++pos;
    ;i<=n;++i)lovewxy[hd[i]]=max(lovewxy[hd[i]],i);
    ;
    while(p>=k)
    --p;
    ;--i){dfs(i);}build(,,n);
    ;i<=n;++i){
        if(fa[i]&&!vis[fa[i]]){
            motify(,,n,ans[fa[i]],n,siz[fa[i]]-);vis[fa[i]]=;
        }
        ,,n,siz[i]);
        tp=lovewxy[hd[tp]];
        --lovewxy[hd[tp]];
        ans[i]=tp;motify(,,n,ans[i],n,-siz[i]);
    }
    printf(]]);
    ;i<=n;++i)printf(" %d",d[ans[i]]);
}