排序列表

题目大意

有若干个区间，C(m)表示所有包含m这个点的区间编号排序后的序列。
求本质不同的非空字典序第k小的序列。

做法

先离散化，因为本质不同不会超过2n个序列。
接下来顺序扫，并维护每个位置的hash值。
遇到之前出现过的hash值就叉掉。
然后接下来枚举按字典序枚举，每次看看往字典序末尾加入i会有多少种可能。
对于k,如果它不在答案序列中，不能选择它区间所包含的m。
对于k，如果它在答案序列中，只能选择它区间所包含的m。
所以需要线段树。
还要注意每次要检查当前答案序列后面能不能不放东西了，这个也可以线段树。

#include<cstdio>#include<algorithm>#include<map>#define fo(i,a,b) for(i=a;i<=b;i++)using namespace std;const int maxn=100000+10,inf=10000000,mo1=1000000007,mo2=998244353;typedef long long ll;typedef pair<int,int> pi;pi zlt;map<pi,bool> s;int father[maxn],tree[maxn][2],mic[maxn][2],ha[maxn][2],size[maxn];int h[maxn*2],go[maxn*2],nxt[maxn*2],fx[maxn*2];int sum[maxn*8],mi[maxn*8],ad[maxn*8];bool bz[maxn*8];int a[maxn],b[maxn],d[maxn*2],ans[maxn];int i,j,k,l,r,t,n,m,tot,top,cnt,root;bool czy;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;}void add(int x,int y,int z){    go[++cnt]=y;    fx[cnt]=z;    nxt[cnt]=h[x];    h[x]=cnt;}void build(int p,int l,int r){    if (l==r){        sum[p]=1;        return;    }    int mid=(l+r)/2;    build(p*2,l,mid);build(p*2+1,mid+1,r);    sum[p]=sum[p*2]+sum[p*2+1];}void markcl(int p){    bz[p]=1;    sum[p]=0;    mi[p]=inf;}void markad(int p,int v){    ad[p]+=v;    mi[p]+=v;}void down(int p){    if (bz[p]){        markcl(p*2);        markcl(p*2+1);        bz[p]=0;    }    if (ad[p]){        markad(p*2,ad[p]);        markad(p*2+1,ad[p]);        ad[p]=0;    }}void change(int p,int l,int r,int a,int b){     if (a>b) return;    if (l==a&&r==b){        markcl(p);        return;    }    down(p);    int mid=(l+r)/2;    if (b<=mid) change(p*2,l,mid,a,b);    else if (a>mid) change(p*2+1,mid+1,r,a,b);    else change(p*2,l,mid,a,mid),change(p*2+1,mid+1,r,mid+1,b);    sum[p]=sum[p*2]+sum[p*2+1];}int query(int p,int l,int r,int a,int b){    if (l==a&&r==b) return sum[p];    down(p);    int mid=(l+r)/2;    if (b<=mid) return query(p*2,l,mid,a,b);    else if (a>mid) return query(p*2+1,mid+1,r,a,b);    else return query(p*2,l,mid,a,mid)+query(p*2+1,mid+1,r,mid+1,b);}void change2(int p,int l,int r,int a,int b,int v){    if (l==a&&r==b){            markad(p,v);        return;    }    down(p);    int mid=(l+r)/2;    if (b<=mid) change2(p*2,l,mid,a,b,v);    else if (a>mid) change2(p*2+1,mid+1,r,a,b,v);    else change2(p*2,l,mid,a,mid,v),change2(p*2+1,mid+1,r,mid+1,b,v);    mi[p]=min(mi[p*2],mi[p*2+1]);}int pd(int x){    return tree[father[x]][1]==x;}void update(int x){    size[x]=size[tree[x][0]]+size[tree[x][1]]+1;    int t=((ll)ha[tree[x][0]][0]*(n+1)%mo1+x)%mo1;    t=((ll)t*mic[size[tree[x][1]]][0]%mo1+ha[tree[x][1]][0])%mo1;    ha[x][0]=t;    t=((ll)ha[tree[x][0]][1]*(n+1)%mo2+x)%mo2;    t=((ll)t*mic[size[tree[x][1]]][1]%mo2+ha[tree[x][1]][1])%mo2;    ha[x][1]=t;}void rotate(int x){    int y=father[x],z=pd(x);    father[x]=father[y];    if (father[y]) tree[father[y]][pd(y)]=x;    tree[y][z]=tree[x][1-z];    if (tree[x][1-z]) father[tree[x][1-z]]=y;    tree[x][1-z]=y;    father[y]=x;    update(y);    update(x);}void splay(int x,int y){    while (father[x]!=y){        if (father[father[x]]!=y)            if (pd(x)==pd(father[x])) rotate(father[x]);else rotate(x);        rotate(x);    }}void insert(int &x,int y){    if (!x){        x=y;        update(x);        return;    }    if (y<x){        insert(tree[x][0],y);        father[tree[x][0]]=x;    }    else{        insert(tree[x][1],y);        father[tree[x][1]]=x;    }    update(x);}int merge(int a,int b){    if (!a||!b) return a+b;    while (tree[a][1]) a=tree[a][1];    splay(a,0);    tree[a][1]=b;    father[b]=a;    update(a);    return a;}int main(){    freopen("list.in","r",stdin);freopen("list.out","w",stdout);    n=read();m=read();    fo(i,1,n) d[++top]=a[i]=read(),d[++top]=b[i]=read();    sort(d+1,d+top+1);    top=unique(d+1,d+top+1)-d-1;    fo(i,1,n){        a[i]=lower_bound(d+1,d+top+1,a[i])-d;        a[i]++;        b[i]=lower_bound(d+1,d+top+1,b[i])-d;        change2(1,1,top,a[i],b[i],1);    }    //change2(1,1,top,1,1,inf);    build(1,1,top);    fo(i,1,n){        add(a[i],i,1);        add(b[i]+1,i,-1);    }    mic[0][0]=mic[0][1]=1;    fo(i,1,n){        mic[i][0]=(ll)mic[i-1][0]*(n+1)%mo1;        mic[i][1]=(ll)mic[i-1][1]*(n+1)%mo2;    }    fo(i,1,top){        t=h[i];        while (t){            if (fx[t]==-1){                splay(go[t],0);                father[tree[go[t]][0]]=father[tree[go[t]][1]]=0;                root=merge(tree[go[t]][0],tree[go[t]][1]);            }            t=nxt[t];        }        t=h[i];        while (t){            if (fx[t]==1){                insert(root,go[t]);                splay(go[t],0);                root=go[t];            }            t=nxt[t];        }        if (!root) change(1,1,top,i,i);        else{            zlt=make_pair(ha[root][0],ha[root][1]);            if (s[zlt]){                change(1,1,top,i,i);            }else s[zlt]=1;        }    }    fo(i,1,n){        t=query(1,1,top,a[i],b[i]);        if (t<m){            m-=t;            change(1,1,top,a[i],b[i]);            change2(1,1,top,a[i],b[i],-1);        }        else{            ans[++tot]=i;            change(1,1,top,1,a[i]-1);            change(1,1,top,b[i]+1,top);            change2(1,1,top,a[i],b[i],-1);            if (mi[1]==0) m--;            if (m==0) break;        }    }    printf("%d\n",tot);    fo(i,1,tot) printf("%d ",ans[i]);    printf("\n");    fclose(stdin);fclose(stdout);    return 0;}