bzoj 3932: [CQOI2015]任务查询系统 可持久化纤段树

题意

最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。超级计算机中的
任务用三元组(Si,Ei,Pi)描述，(Si,Ei,Pi)表示任务从第Si秒开始，在第Ei秒后结束（第Si秒和Ei秒任务也在运行
），其优先级为Pi。同一时间可能有多个任务同时执行，它们的优先级可能相同，也可能不同。调度系统会经常向
查询系统询问，第Xi秒正在运行的任务中，优先级最小的Ki个任务（即将任务按照优先级从小到大排序后取前Ki个
）的优先级之和是多少。特别的，如果Ki大于第Xi秒正在运行的任务总数，则直接回答第Xi秒正在运行的任务优先
级之和。上述所有参数均为整数，时间的范围在1到n之间（包含1和n）。

分析

对操作查分，然后按照时间建主席树即可。

代码

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#define N 100005#define LL long longusing namespace std;int n,m,cnt,tot,last[N],root[N];struct edge{int val,p,next;}e[N*2];struct tree{int l,r,num;LL val;}t[N*50];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 ins(int &d,int p,int l,int r,int x,int y){    d=++tot;    t[d].l=t[p].l;t[d].r=t[p].r;t[d].num=t[p].num;t[d].val=t[p].val;    t[d].val+=x*y;t[d].num+=y;    if (l==r) return;    int mid=(l+r)/2;    if (x<=mid) ins(t[d].l,t[p].l,l,mid,x,y);    else ins(t[d].r,t[p].r,mid+1,r,x,y);}LL query(int d,int l,int r,int k){    if (l==r) return (LL)k*l;    if (t[d].num<=k) return t[d].val;    int mid=(l+r)/2;    if (t[t[d].l].num<k) return t[t[d].l].val+query(t[d].r,mid+1,r,k-t[t[d].l].num);    else return query(t[d].l,l,mid,k);}int main(){    m=read();n=read();    for (int i=1;i<=m;i++)    {        int s=read(),E=read(),p=read();        e[++cnt].val=1;e[cnt].p=p;e[cnt].next=last[s];last[s]=cnt;        e[++cnt].val=-1;e[cnt].p=p;e[cnt].next=last[E+1];last[E+1]=cnt;    }    for (int i=1;i<=n;i++)    {        root[i]=root[i-1];        for (int j=last[i];j;j=e[j].next)            ins(root[i],root[i],1,10000000,e[j].p,e[j].val);    }    LL pre=1;    for (int i=1;i<=n;i++)    {        int k=read(),a=read(),b=read(),c=read();        pre=query(root[k],1,10000000,1+(LL)((LL)a*pre+(LL)b)%c);        printf("%lld\n",pre);    }    return 0;}