[LibreOJ β Round]ZQC的手办

题目大意

区间对一个数取max。
区间求最小的x个比k小的数。

线段树

第一个操作很好搞。
第二个操作有个很显然的常数大做法。
实际上可以用堆把这个区间的笛卡尔树按优先级广搜。

#include<cstdio>
#include<algorithm>
#include<queue>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=500000+10,inf=1000000000;
struct dong{
    int x,y;
    friend bool operator <(dong a,dong b){
        return a.x<b.x||a.x==b.x&&a.y<b.y;
    }
    friend dong operator +(dong a,dong b){
        if (a<b) return a;else return b;
    }
} tmp;
struct suan{
    int l,r,p,v;
    friend bool operator <(suan a,suan b){
        return a.v>b.v||a.v==b.v&&a.p<b.p;
    }
} zlt;
priority_queue<suan> dl;
dong tree[maxn*4];
int mx[maxn*4],a[maxn],ans[maxn],sta[80];
bool bz[maxn*4];
int i,j,k,l,r,p,t,n,m,tot,top;
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;
}
void build(int p,int l,int r){
    if (l==r){
        tree[p].x=a[l];
        tree[p].y=l;
        return;
    }
    int mid=(l+r)/2;
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    tree[p]=tree[p*2]+tree[p*2+1];
}
void mark(int p,int v){
    bz[p]=1;
    mx[p]=max(mx[p],v);
    if (tree[p].x<v) tree[p].x=v;
}
void down(int p){
    if (bz[p]){
        mark(p*2,mx[p]);
        mark(p*2+1,mx[p]);
        bz[p]=0;
    }
}
void change(int p,int l,int r,int a,int b,int v){
    if (l==a&&r==b){
        mark(p,v);
        return;
    }
    down(p);
    int mid=(l+r)/2;
    if (b<=mid) change(p*2,l,mid,a,b,v);
    else if (a>mid) change(p*2+1,mid+1,r,a,b,v);
    else change(p*2,l,mid,a,mid,v),change(p*2+1,mid+1,r,mid+1,b,v);
    tree[p]=tree[p*2]+tree[p*2+1];
}
void query(int p,int l,int r,int a,int b){
    if (l==a&&r==b){
        tmp=tmp+tree[p];
        return;
    }
    down(p);
    int mid=(l+r)/2;
    if (b<=mid) query(p*2,l,mid,a,b);
    else if (a>mid) query(p*2+1,mid+1,r,a,b);
    else query(p*2,l,mid,a,mid),query(p*2+1,mid+1,r,mid+1,b);
}
int main(){
    n=read();
    fo(i,1,n) a[i]=read();
    build(1,1,n);
    m=read();
    fo(i,1,m){
        t=read();
        if (t==1){
            l=read();r=read();k=read();
            change(1,1,n,l,r,k);
        }
        else{
            l=read();r=read();k=read();j=read();
            if (r-l+1<j){
                printf("-1\n");
                continue;
            }
            while (!dl.empty()) dl.pop();
            zlt.l=l;zlt.r=r;
            tmp.x=inf;
            query(1,1,n,l,r);
            zlt.p=tmp.y;zlt.v=tmp.x;
            dl.push(zlt);
            tot=0;
            while (j--){
                zlt=dl.top();
                dl.pop();
                if (zlt.v>=k) break;
                ans[++tot]=zlt.v;
                l=zlt.l;r=zlt.r;p=zlt.p;
                if (l<p){
                    zlt.l=l;
                    zlt.r=p-1;
                    tmp.x=inf;
                    query(1,1,n,l,p-1);
                    zlt.p=tmp.y;zlt.v=tmp.x;
                    dl.push(zlt);
                }
                if (p<r){
                    zlt.l=p+1;
                    zlt.r=r;
                    tmp.x=inf;
                    query(1,1,n,p+1,r);
                    zlt.p=tmp.y;zlt.v=tmp.x;
                    dl.push(zlt);
                }
            }
            if (j>=0) printf("-1\n");
            else{
                fo(j,1,tot) printf("%d ",ans[j]);
                printf("\n");
            }
        }
    }
}