主席树||可持久化线段树||离散化||[CQOI2015]任务查询系统||BZOJ 3932||Luogu P3168

题目： [CQOI2015]任务查询系统

题解：

是一道很经典的题目。大体思路是抓优先级来当下标做主席树，用时刻作为主席树的版本。然而优先级范围到1e7去了，就离散化一遍。然后把每个事件的开始（s）、结束（e）（e记得+1，因为一个事件是第e+1时刻结束的）时间点抓出来排序一遍，按时刻从早到晚维护主席树。差不多就是这样。

代码：

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define ll long long
 #define max(a,b) ((a)>(b)?(a):(b))
 #define min(a,b) ((a)<(b)?(a):(b))
 using namespace std;
 const int maxn=+,maxm=maxn,maxnum_p=maxm,maxnum_b=maxm<<;
 int N,M,cnt_p=,num_b=,num_treenode,root[maxnum_b],X,A,B,C,belong_root[maxn];
 ll Pre=,belong[maxm],K;
 struct _A{int p,s,e,hp;}a[maxm];
 bool cmp1(const _A&c,const _A&d){return(c.p<d.p);}
 struct _B{int p,o,data;}b[maxm<<];//b数组负责记录把每个时间点都拆出来后的数据 ，data记录时间点，o是时间点的类型
 bool cmp2(const _B&c,const _B&d){return(c.data<d.data);}
 struct Tree{int l,r,ls,rs,cnt;ll sum;}t[(maxnum_p<<)+maxnum_b*];
 void Build(int x,int l,int r){
     t[x].l=l;t[x].r=r;int mid=(l+r)>>;
     if(l==r)return;
     Build(t[x].ls=++num_treenode,l,mid);Build(t[x].rs=++num_treenode,mid+,r);
     return;
 }
 void Update(int u,int x,int p,int o){
     int l=t[u].l,r=t[u].r,mid=(l+r)>>;
     t[x].l=l;t[x].r=r;
     if(l==r&&l==p){
         if(o==){t[x].cnt=t[u].cnt+; t[x].sum=t[u].sum+belong[p];}
             else{t[x].cnt=t[u].cnt-; t[x].sum=t[u].sum-belong[p];}
         return;
     }
     if(p<=mid){t[x].rs=t[u].rs; Update(t[u].ls,t[x].ls=++num_treenode,p,o);}
           else{t[x].ls=t[u].ls; Update(t[u].rs,t[x].rs=++num_treenode,p,o);}
     t[x].cnt=t[t[x].ls].cnt+t[t[x].rs].cnt;
     t[x].sum=t[t[x].ls].sum+t[t[x].rs].sum;
     return;
 }
 ll Query(int x,ll k){
     if(k==)return ;
     if(k>=t[x].cnt)return t[x].sum;
     int l=t[x].l,r=t[x].r,ls=t[x].ls,rs=t[x].rs;
     if(l==r)return (k*belong[l]);
     if(t[ls].cnt>=k)return Query(ls,k);
     else return (t[ls].sum+Query(rs,k-t[ls].cnt));
 }
 int main(){
     scanf("%d%d",&M,&N);
     for(int i=;i<=M;i++)scanf("%d%d%d",&a[i].s,&a[i].e,&a[i].p),a[i].e++;
     sort(a+,a+M+,cmp1);
     for(int i=;i<=M;i++)
         if(a[i].p!=a[i-].p)belong[a[i].hp=++cnt_p]=a[i].p;else a[i].hp=cnt_p;
     for(int i=;i<=M;i++){
         b[++num_b].data=a[i].s;b[num_b].o=;b[num_b].p=a[i].hp;
         b[++num_b].data=a[i].e;b[num_b].o=;b[num_b].p=a[i].hp;
     }
     Build(root[]=++num_treenode,,cnt_p);
     sort(b+,b+num_b+,cmp2);
     for(int i=;i<=num_b;i++){
         int v=b[i].data,o=b[i].o,p=b[i].p;
         Update(root[i-],root[i]=++num_treenode,p,o);
         if(v!=b[i+].data)belong_root[v]=i;
     }
     for(int i=;i<=N;i++)if(belong_root[i]==)belong_root[i]=belong_root[i-];
     while(N--){
         scanf("%d%d%d%d",&X,&A,&B,&C);
         K=+(Pre*A+B)%C;
         Pre=Query(root[belong_root[X]],K);
         printf("%lld\n",Pre);
     }
     return ;
 }

---

By:AlenaNuna