题目大意:
最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分。超级计算机中的
任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si秒开始,在第Ei秒后结束(第Si秒和Ei秒任务也在运行
),其优先级为Pi。同一时间可能有多个任务同时执行,它们的优先级可能相同,也可能不同。调度系统会经常向
查询系统询问,第Xi秒正在运行的任务中,优先级最小的Ki个任务(即将任务按照优先级从小到大排序后取前Ki个
)的优先级之和是多少。特别的,如果Ki大于第Xi秒正在运行的任务总数,则直接回答第Xi秒正在运行的任务优先
级之和。上述所有参数均为整数,时间的范围在1到n之间(包含1和n)。
思路:
先将优先值离散,然后对每个时间建一棵记录当前时间正在运行的优先值为l~r的任务的个数、优先值之和的线段树,显然可以用主席树优化。一个从s到t的任务可以看成在s时间发生,在t+1时间结束,于是一个任务只需要在s和t+1更新。
对于查找,判断当前要找的k和左子树任务个数的大小关系,如果k大,查询右子树,否则查询左子树。注意当l==r时返回要注意k的值。
代码:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define N 100001
#define ll long long
inline void Read(long long& x){
char c=getchar();
for(;c<''||c>'';c=getchar());
for(x=;c>=''&&c<='';x=(x<<)+(x<<)+c-,c=getchar());
}
inline void Read(int& x){
char c=getchar();
for(;c<''||c>'';c=getchar());
for(x=;c>=''&&c<='';x=(x<<)+(x<<)+c-,c=getchar());
}
vector<int>a1[N],a2[N];
struct node{
int l,r,s;
ll w;
}c[N*];
struct gj{
int x,y,z;
}b[N];
struct ls{
int w,f;
}a[N];
int i,j,Rt[N],k,Num,n,m,M;
ll Last=,w2[N],x,y,z;
inline void Update(int& Node,int l,int r,int x,int y){
c[++Num]=c[Node];Node=Num;
c[Node].w+=w2[x]*y;c[Node].s+=y;
if(l==r)return;
int Mid=(l+r)>>;
if(x<=Mid)Update(c[Node].l,l,Mid,x,y);else Update(c[Node].r,Mid+,r,x,y);
}
inline ll Query(int Node,int l,int r,int x){
if(l==r)return c[Node].w/c[Node].s*x;
int Mid=(l+r)>>;
if(c[c[Node].l].s>x)return Query(c[Node].l,l,Mid,x);else if(c[c[Node].l].s==x)return c[c[Node].l].w;else return Query(c[Node].r,Mid+,r,x-c[c[Node].l].s)+c[c[Node].l].w;
}
inline bool Cmp(ls a,ls b){
return a.w<b.w;
}
int main()
{
Read(m);Read(n);
for(i=;i<=m;i++)Read(b[i].x),Read(b[i].y),Read(b[i].z),a[i].w=b[i].z,a[i].f=i;
sort(a+,a+m+,Cmp);
for(i=;i<=m;i++)if(a[i].w==a[i-].w)b[a[i].f].z=M;else b[a[i].f].z=++M,w2[M]=a[i].w;
for(i=;i<=m;i++){
a1[b[i].x].push_back(b[i].z);a2[b[i].y+].push_back(b[i].z);
}
for(i=;i<=n;i++){
Rt[i]=Rt[i-];
for(j=;j<a1[i].size();j++)Update(Rt[i],,M,a1[i][j],);
for(j=;j<a2[i].size();j++)Update(Rt[i],,M,a2[i][j],-);
}
for(i=;i<=n;i++){
Read(k);Read(x);Read(y);Read(z);
x=((x*Last+y)%z)+;
if(c[Rt[k]].s<=x)Last=c[Rt[k]].w;else Last=Query(Rt[k],,M,x);
printf("%lld\n",Last);
}
return ;
}
bzoj3932