题目描述
最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。超级计算机中的任务用三元组(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行，每行一个整数，表示查询结果。
输入输出样例
说明
样例解释
K1 = (1*1+3)%2+1 = 1
K2 = (1*2+3)%4+1 = 2
K3 = (2*8+4)%3+1 = 3
对于100%的数据，1<=m,n,Si,Ei,Ci<=100000，0<=Ai,Bi<=100000，1<=Pi<=10000000，Xi为1到n的一个排列
按p排序，按时间建立主席树，维护区间内的数个数，权值和
每一次查询前二分，二分的判断用一个主席树查询x上数的个数
找到了最小的满足x点上数个数>=k的点，主席树上单点查询和
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 using namespace std; 6 typedef long long lol; 7 struct ZYYS 8 { 9 int l,r; 10 lol p; 11 }P[100001]; 12 int pos,ch[5000001][2]; 13 lol sum[5000001],lazy[5000001]; 14 int size[5000001],n,m,root[200001]; 15 lol ans; 16 bool cmp(ZYYS a,ZYYS b) 17 { 18 return a.p<b.p; 19 } 20 void update(int &rt,int l,int r,int L,int R,lol k) 21 { 22 int Rt=rt; 23 rt=++pos; 24 ch[rt][0]=ch[Rt][0];ch[rt][1]=ch[Rt][1]; 25 sum[rt]=sum[Rt];size[rt]=size[Rt]; 26 lazy[rt]=lazy[Rt]; 27 if (l>=L&&r<=R) 28 { 29 sum[rt]+=(r-l+1)*k; 30 lazy[rt]+=k; 31 size[rt]+=1; 32 return; 33 } 34 int mid=(l+r)/2; 35 if (L<=mid) update(ch[rt][0],l,mid,L,R,k); 36 if (R>mid) update(ch[rt][1],mid+1,r,L,R,k); 37 sum[rt]=sum[ch[rt][0]]+sum[ch[rt][1]]+(r-l+1)*lazy[rt]; 38 } 39 int querysize(int rt,int l,int r,int x,int la) 40 { 41 if (l==r) return size[rt]+la; 42 int mid=(l+r)/2; 43 if (x<=mid) return querysize(ch[rt][0],l,mid,x,la+size[rt]); 44 else return querysize(ch[rt][1],mid+1,r,x,la+size[rt]); 45 } 46 lol querysum(int rt,int l,int r,int x,lol laz) 47 { 48 if (l==r) 49 return sum[rt]+laz; 50 int mid=(l+r)/2; 51 if (x<=mid) return querysum(ch[rt][0],l,mid,x,laz+lazy[rt]); 52 else return querysum(ch[rt][1],mid+1,r,x,laz+lazy[rt]); 53 } 54 int find(int x,int k) 55 { 56 int l=1,r=n,as=n; 57 while (l<=r) 58 { 59 int mid=(l+r)/2; 60 if (querysize(root[mid],1,m,x,0)>=k) as=mid,r=mid-1; 61 else l=mid+1; 62 } 63 return as; 64 } 65 int main() 66 {int i,s,e,x; 67 lol p,a,b,c; 68 cin>>n>>m; 69 for (i=1;i<=n;i++) 70 { 71 scanf("%d%d%lld",&s,&e,&p); 72 P[i].l=s;P[i].r=e;P[i].p=p; 73 } 74 sort(P+1,P+n+1,cmp); 75 for (i=1;i<=n;i++) 76 { 77 root[i]=root[i-1]; 78 update(root[i],1,m,P[i].l,P[i].r,P[i].p); 79 } 80 ans=1; 81 for (i=1;i<=m;i++) 82 { 83 scanf("%d%lld%lld%lld",&x,&a,&b,&c); 84 lol k=1+(a*ans%c+b)%c; 85 int t=find(x,k); 86 ans=querysum(root[t],1,m,x,0); 87 printf("%lld ",ans); 88 } 89 }