[九省联考 2018] iiidx

题目描述：

给出一个序列
重新安排序列的下标
使得 $d_i<=d_{i/k}$ 且字典序最大

题目分析：

出题人的题解好皮啊
30分：
爆搜啊啊啊啊
55分：
贪心哇.
我们令 fa[i]=i/k 那么就形成了一个森林（或者是以0为根的树，后序遍历填数
100分：
上面的贪心仅对$d_i$不相同的情况成立
比如
4 2.0
1 1 1 2 这组样例 用贪心搞出来的是
1 1 1 2
其实正解应该是1 1 2 1
如果权值有重复的，我们要找和第siz大的权值相等的权值中最靠右的那个点，再给子树预定siz-1个，这样做能够让当前点权值尽可能大的前提下给子树预定的权值尽可能的小，这样，万一下一个要分配的点不在这个点的子树中，也能被分配尽可能大的权值。
先 降序排序
考虑搞一个线段树，维护每个数左边可用的数字个数f[i]以及f[i]的最小值
当我们操作一个树根时，我们只需要从 树根的节点处-n 区间-siz[树根] 表示子树预定了这些就好了，尽量在树上找靠右的
不要忘了搞到儿子时消除父亲的影响
为什么BZOJ卡了精度洛谷却没有？数据不应该是官方数据吗。。。

题目链接：

BZOJ 5249
Luogu 4364

Ac 代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
const int maxm=510000;
struct node{
    int l,r,del,lm;
}st[maxm<<2];
int a[maxm],ans[maxm],siz[maxm],cnt[maxm];
int n;
double k;
int comp(int a,int b)
{
    return a>b;
}
void pushup(int o)
{
    st[o].lm=std::min(st[(o<<1)].lm,st[(o<<1)|1].lm);
}
void col(int o,int num)
{
    st[o].lm+=num,st[o].del+=num;
}
void pushdown(int o)
{
    if(st[o].del)
    {
        col((o<<1),st[o].del);
        col((o<<1)|1,st[o].del);
        st[o].del=0;
    }
}
void build(int o,int l,int r)
{
    st[o].l=l,st[o].r=r;
    if(l>=r)
    {
        st[o].lm=l;
        return;
    }
    int mid=(l+r)>>1;
    build((o<<1),l,mid),build((o<<1)|1,mid+1,r);
    pushup(o);
}
void modify(int o,int l,int r,int ql,int qr,int num)
{
    if(ql<=l&&r<=qr)
    {
        col(o,num);
        return;
    }
    pushdown(o);
    int mid=(l+r)>>1;
    if(ql<=mid) modify((o<<1),l,mid,ql,qr,num);
    if(qr>mid) modify((o<<1)|1,mid+1,r,ql,qr,num);
    pushup(o);
}
int query(int o,int l,int r,int rank)
{
    if(l>=r) return st[o].lm>=rank?l:l+1;
    pushdown(o);
    int mid=(l+r)>>1;
    if(st[(o<<1)|1].lm>=rank) return query((o<<1),l,mid,rank);
    else return query((o<<1)|1,mid+1,r,rank);
}
int fa(int x)
{
    return (int)floor(x/k);
}
int main()
{
    scanf("%d%lf",&n,&k);
    for(int i=1;i<=n;i++)
     scanf("%d",&a[i]),siz[i]=1;
    build(1,1,n);
    std::sort(a+1,a+n+1,comp);
    for(int i=n-1;i;i--) if(a[i]==a[i+1]) cnt[i]=cnt[i+1]+1;
    for(int i=n;i;i--) siz[fa(i)]+=siz[i];
    for(int i=1;i<=n;i++)
    {
        if(fa(i)&&fa(i)!=fa(i-1)) modify(1,1,n,ans[fa(i)],n,siz[fa(i)]-1);
        ans[i]=query(1,1,n,siz[i]);
        ans[i]+=cnt[ans[i]],ans[i]-=cnt[ans[i]]++;
        modify(1,1,n,ans[i],n,-siz[i]);
    }
    for(int i=1;i<=n;i++) printf("%d ",a[ans[i]]);
    return 0;
}