BZOJ5249: [2018多省省队联测]IIIDX

$n leq 500000$个数字，给实数$k$，问用这些数字填上的，满足$d_i geq d_{left lfloor frac{i}{k} ight floor}$的字典序最大的序列。

如果数字不同的话按后序遍历从大到小填就行了。

如果数字相同的话就不行。比如输入4 2 1 1 1 2。

预定。每选一个数，要预定几个比他大的数。首先从小到大排序。

方法零：线段树上维护每个数选了没选。在同一层的时候，每选一个数，就预定他左边紧贴着他的，他子树大小个数的数字。然后这层选完再把它还原回去。用个可持久化可以还原回去。难写。

方法一：线段树上维护每个数右边有几个数还没预定。每次找数需要找一个x，满足x左边的“未预定数”都大于等于x的子树大小，找权值尽量大的、同权值尽量靠左的数字选中。选完后，把它左边所有数的“未预定数”减去他的子树大小。决定一个点时取消它父亲的影响，注意只用取消一次。

//#include<iostream>
#include<cstring>
#include<cstdio>
//#include<time.h>
//#include<complex>
//#include<queue>
#include<algorithm>
#include<stdlib.h>
using namespace std;

#define LL long long
int qread()
{
    char c; int s=0; while ((c=getchar())<'0' || c>'9');
    do s=s*10+c-'0'; while ((c=getchar())>='0' && c<='9'); return s;
}

//Pay attention to '-' , LL and double of qread!!!!

int n; double K;
#define maxn 500011
int num[maxn];
struct Edge{int to,next;}edge[maxn]; int first[maxn],le=2;
void in(int x,int y) {Edge &e=edge[le]; e.to=y; e.next=first[x]; first[x]=le++;}

int size[maxn];
void dfs(int x)
{
    size[x]=1;
    for (int i=first[x];i;i=edge[i].next)
    {
        Edge &e=edge[i];
        dfs(e.to); size[x]+=size[e.to];
    }
}

struct SMT
{
    struct Node{int ls,rs; int Min,add;}a[maxn<<1];
    int size,n;
    void up(int x) {a[x].Min=min(a[a[x].ls].Min,a[a[x].rs].Min);}
    void build(int &x,int L,int R)
    {
        x=++size;
        if (L==R) {a[x].Min=n-L; return;}
        int mid=(L+R)>>1;
        build(a[x].ls,L,mid); build(a[x].rs,mid+1,R);
        up(x);
    }
    void build(int m) {n=m; size=0; int x; build(x,1,n);}
    
    void addsingle(int x,int v) {a[x].add+=v; a[x].Min+=v;}
    void down(int x) {if (a[x].add) {addsingle(a[x].ls,a[x].add); addsingle(a[x].rs,a[x].add); a[x].add=0;}}
    
    int ql,qr,v;
    void Add(int x,int L,int R)
    {
        if (ql<=L && R<=qr) {addsingle(x,v); return;}
        down(x);
        int mid=(L+R)>>1;
        if (ql<=mid) Add(a[x].ls,L,mid);
        if (qr>mid) Add(a[x].rs,mid+1,R);
        up(x);
    }
    void add(int L,int R,int v) {if (L>R) return; ql=L; qr=R; this->v=v; Add(1,1,n);}
    
    int Query(int x,int L,int R)
    {
        if (L==R) return L;
        down(x);
        int mid=(L+R)>>1;
        if (a[a[x].ls].Min>=v) return Query(a[x].rs,mid+1,R);
        return Query(a[x].ls,L,mid);
    }
    int query(int val) {v=val; return Query(1,1,n);}
    
    int Find(int x,int L,int R)
    {
        if (L==R) return L;
        down(x);
        int mid=(L+R)>>1;
        if (num[mid]>=v) return Find(a[x].ls,L,mid);
        return Find(a[x].rs,mid+1,R);
    }
    int find(int val) {v=val; return Find(1,1,n);}
    
}t;

int ans[maxn],fa[maxn];
int main()
{
    n=qread(); scanf("%lf",&K);
    for (int i=1;i<=n;i++) num[i]=qread();
    
    fa[0]=-1;
    for (int i=1;i<=n;i++) in(fa[i]=(int)(i/K+1e-9),i);
    dfs(0);
    
    sort(num+1,num+1+n); num[++n]=0x3f3f3f3f;
    t.build(n);
    for (int i=1;i<n;i++)
    {
        if (fa[i] && fa[i]!=fa[i-1]) t.add(1,ans[fa[i]],size[fa[i]]-1);
        int val=num[t.query(size[i])-1]; int p=t.find(val);
        ans[i]=p; t.add(1,ans[i],-size[i]);
    }
    
    for (int i=1;i<n;i++) printf("%d ",num[ans[i]]);
    return 0;
}

方法二：看不太懂