[CQOI2015]任务查询系统

题目描述
最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。
超级计算机中的任务用三元组(Si,Ei,Pi)描述，(Si,Ei,Pi)表示任务从第Si秒开始，在第Ei秒后结束（第Si秒和Ei秒任务也在
运行），其优先级为Pi。同一时间可能有多个任务同时执行，它们的优先级可能相同，也可能不同。
调度系统会经常向查询系统询问，第Xi秒正在运行的任务中，优先级最小的Ki个任务（即将任务按照优先级从小到大
排序后取前Ki个）的优先级之和是多少。特别的，如果Ki大于第Xi秒正在运行的任务总数，则直接回答第Xi秒正在运行
的任务优先级之和。上述所有参数均为整数，时间的范围在1到n之间（包含1和n）。
输入格式
输入文件第一行包含两个空格分开的正整数m和n，分别表示任务总数和时间范围。
接下来m行，每行包含三个空格分开的正整数Si、Ei和Pi(Si≤Ei)，描述一个任务。
接下来n行，每行包含四个空格分开的整数Xi、Ai、Bi和Ci，描述一次查询。查询的参数Ki需要由公式 Ki=1 +
(Ai*Pre+Bi) mod Ci 计算得到。其中Pre表示上一次查询的结果，对于第一次查询，Pre=1。
输出格式
输出共n行，每行一个整数，表示查询结果。
输入样例
4 3
1 2 6
2 3 3
1 3 2
3 3 4
3 1 3 2
1 1 3 4
2 2 4 3
输出样例
281
1
样例解释
K1 = (1*1+3)%2+1 = 1
K2 = (1*2+3)%4+1 = 2
K3 = (2*8+4)%3+1 = 3
数据范围
对于50%的数据，Ai=0

对于100%的数据，1≤m,n,Si,Ei,Ci≤100000，0≤Ai,Bi≤100000，1≤Pi≤10000000，Xi为1到n的一个排列

题解：

题上说了要求区间第k大 可以考虑将一条线段变成两个点的差分，那么1-i的状态之和就是第i个时刻的状态 那么就可以以时间为版本建立主席树 每一棵线段树维护区间个数与优先度之和 下标为优先度离散后的值 查询优先度之和就行了

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define ls(x) (a[x].lc)
#define rs(x) (a[x].rc)
#define LL long long
const int maxn=150000*2;
struct Node{
    int lc,rc,sum;LL date;
}a[maxn*22];
struct Task{
    int x,p,type;
    bool operator < (const Task & a)const{
        return x<a.x;
    }
}q[maxn];
int n,m,cnt,num;
int hash[maxn],root[maxn];
void Insert(int x,int add,int & rt,int l,int r){
    a[++cnt]=a[rt];rt=cnt;
    a[rt].sum+=add;
    a[rt].date+=hash[x]*add;
    if(l==r)return;
    int mid=(l+r)>>1;
    if(x<=mid)Insert(x,add,ls(rt),l,mid);
    else Insert(x,add,rs(rt),mid+1,r);
}
int Get_pos(int k,int rt,int l,int r){
    if(l==r)return l;
    int mid=(l+r)>>1;
    int temp=a[ls(rt)].sum;
    if(temp>=k)return Get_pos(k,ls(rt),l,mid);
    else return Get_pos(k-temp,rs(rt),mid+1,r);
}
LL Query(int pos,int k,int rt,int l,int r){
    if(r<pos)return a[rt].date;
    if(l==r && r==pos)return a[rt].date/a[rt].sum*k;
    int mid=(l+r)>>1;
    LL res=Query(pos,k,ls(rt),l,mid);
    if(pos>mid)res+=Query(pos,k-a[ls(rt)].sum,rs(rt),mid+1,r);
    return res;
}
void Init();
int main(){
    freopen("cqoi15_query.in","r",stdin);freopen("cqoi15_query.out","w",stdout);
    Init();
    getchar();getchar();
    return 0;
}
void Init(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        int s,t,p;scanf("%d%d%d",&s,&t,&p);
        q[i].x=s;q[i].p=p;q[i].type=1;
        q[n+i].x=t+1;q[n+i].p=p;q[n+i].type=-1;
        hash[i]=p;
    }
    sort(hash+1,hash+n+1);
    num=unique(hash+1,hash+n+1)-hash-1;
    sort(q+1,q+2*n+1);
    for(int i=1;i<=n*2;i++){
        int x=lower_bound(hash+1,hash+num+1,q[i].p)-hash;
        for(int j=q[i-1].x+1;j<=q[i].x;j++)root[j]=root[j-1];
        Insert(x,q[i].type,root[q[i].x],1,num);
    }
    for(int i=q[2*n].x+1;i<=m;i++)root[i]=root[i-1];
    LL ans=1;
    for(int i=1;i<=m;i++){
        int x,aa,b,c;scanf("%d%d%d%d",&x,&aa,&b,&c);
        int k=(1ll*ans%c*1ll*aa%c+b)%c+1;
        if(k>a[root[x]].sum)ans=a[root[x]].date;
        else if(k==0)ans=0;
        else {
            int pos=Get_pos(k,root[x],1,num);
            ans=Query(pos,k,root[x],1,num);
        }
        printf("%lld\n",ans);
    }
}