bzoj5249 [2018多省省队联测]IIIDX（树+贪心+线段树）

把i的父亲设为i/k，我们显然会得到一座森林。（或者说以0为根的一棵树）

我们要把给定的权值分配给每一个点，满足每个点上的权值都$\leq$ 子树内点的权值，并使这些点编号从小到大的权值字典序最大。

我们贪心的想，把权值从小到大排序，每次给能给的最靠后位置即可。

但是如果有相同的这样可能会有问题，比如：
n=4,k=2,d={1,1,1,2}。
我们贪心会得到:{1,1,1,2}
实际上有更好的:{1,1,2,1}

问题出在哪呢？问题在于你本来后一个数一定比当前这个数大，所以当前这个数最靠后一定最优。但是现在你后一个数可能和当前数一样大，但是你如果先放一个不得不放的靠前的位置，那么可能可以扩展出更靠后的位置来放当前这个数。

因此现在每次做的时候，你要找不得不放的最靠前的位置来访。
比如{1,1,1,2}，第一个1只能放在1号点，考虑下一个1，就算放在3，下一个还是要放在2，不如先放2，看能不能扩展出更靠后的，比如说这个例子中扩展出了更靠后的位置4，这样就把位置3留给了更大的数。

tips：bzoj O2 毒瘤卡精度。需要eps=1e-8

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
#define eps 1e-8
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,a[N],c[N],sz[N],ans[N];double k;
vector<int>son[N];
struct node{
    int x;
}tr[N<<2];
inline void add(int p,int l,int r,int x,int val){
    tr[p].x+=val;if(l==r) return;int mid=l+r>>1;
    if(x<=mid) add(p<<1,l,mid,x,val);
    else add(p<<1|1,mid+1,r,x,val);
}
inline int ask(int p,int l,int r,int k){//求第k大 
    if(l==r) return l;
    int mid=l+r>>1;
    if(k<=tr[p<<1|1].x) return ask(p<<1|1,mid+1,r,k);
    return ask(p<<1,l,mid,k-tr[p<<1|1].x);
}
int main(){
//  freopen("iiidx.in","r",stdin);
    n=read();scanf("%lf",&k);
    for(int x=n;x>=1;--x){
        int y=x/k+eps;son[y].push_back(x);sz[x]++;sz[y]+=sz[x];a[x]=read();
    }sort(a+1,a+n+1);
    for(int i=n;i>=1;--i) if(a[i]==a[i+1]) c[i]=c[i+1]+1;else c[i]=1;
    for(int i=0;i<son[0].size();++i) add(1,1,n,son[0][i],sz[son[0][i]]);
    for(int i=1;i<=n;++i){
        int x=ask(1,1,n,c[i]);ans[x]=a[i];
        add(1,1,n,x,-sz[x]);
        for(int i=0;i<son[x].size();++i) add(1,1,n,son[x][i],sz[son[x][i]]);
    }for(int i=1;i<=n;++i) printf("%d ",ans[i]);puts("");
    return 0;
}