【BZOJ5249】IIIDX（贪心，线段树）

题意：

思路：赛季结束之前余总推荐的一道好题，不愧是余总

From https://www.cnblogs.com/suika/p/8748115.html

简略的说就是在预留足够多的位置的前提下贪心取最大的数字

剩余可以使用的数字可以使用线段树维护，每次查询可以使用的最大的数字也可以在线段树上二分

维护相同的数字未被使用的的最大的位置的技巧可以学习一下

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef long double ld;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 typedef vector<PII> VII;
 typedef pair<ll,ll>P;
 #define N  500010
 #define M  1000000
 #define INF 1e9
 #define fi first
 #define se second
 #define MP make_pair
 #define pb push_back
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1
 #define fors(i) for(auto i:e[x]) if(i!=p)

 const int MOD=1e8+,inv2=(MOD+)/;
       int p=1e4+;
       double eps=1e-;
       int dx[]={-,,,};
       int dy[]={,,-,};

 struct node
 {
     int s,tag;
 }t[N<<];

 int a[N],b[N],c[N],fa[N],ans[N],sz[N];

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 ll readll()
 {
    ll v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 bool cmp(int a,int b)
 {
     return a>b;
 }

 void pushup(int p)
 {
     t[p].s=min(t[ls].s,t[rs].s);
 }

 void pushdown(int p)
 {
     if(t[p].tag)
     {
         t[ls].s+=t[p].tag;
         t[rs].s+=t[p].tag;
         t[ls].tag+=t[p].tag;
         t[rs].tag+=t[p].tag;
         t[p].tag=;
     }
 }

 void build(int l,int r,int p)
 {
     if(l==r)
     {
         t[p].s=l;
         t[p].tag=;
         return;
     }
     int mid=(l+r)>>;
     build(l,mid,ls);
     build(mid+,r,rs);
     pushup(p);
 }

 void update(int l,int r,int x,int y,int d,int p)
 {
     if(x<=l&&r<=y)
     {
         t[p].s+=d;
         t[p].tag+=d;
         return;
     }
     pushdown(p);
     int mid=(l+r)>>;
     if(x<=mid) update(l,mid,x,y,d,ls);
     if(y>mid) update(mid+,r,x,y,d,rs);
     pushup(p);
 }

 int query(int l,int r,int x,int p)
 {
     if(l==r)
     {
         if(t[p].s>=x) return l;
          else return (l+);
     }
     pushdown(p);
     int mid=(l+r)>>;
     if(t[rs].s>=x) return query(l,mid,x,ls);
     else return query(mid+,r,x,rs);
 }

 int main()
 {
     int n=read();
     double K;
     scanf("%lf",&K);
     rep(i,,n) a[i]=read(),sz[i]=;
     sort(a+,a+n+,cmp);
     per(i,n-,)
      if(a[i]==a[i+]) c[i]=c[i+]+;
       else c[i]=;
     rep(i,,n) fa[i]=(int)((double)i/K+eps);
     per(i,n,) sz[fa[i]]+=sz[i];
     build(,n,);
     rep(i,,n)
     {
         if(fa[i]&&!b[fa[i]])
         {
             update(,n,ans[fa[i]],n,sz[fa[i]]-,);
             b[fa[i]]=;
         }
         int x=query(,n,sz[i],);
         x+=c[x]; c[x]++; x-=(c[x]-); ans[i]=x;
         update(,n,x,n,-sz[i],);
     }
     rep(i,,n-) printf("%d ",a[ans[i]]);
     printf("%d\n",a[ans[n]]);
     return ;
 }