Description
最近实验室正在为其管理的超级计算机编制一套任务管理系统,而你被安排完成其中的查询部分。超级计算机中的
任务用三元组(Si,Ei,Pi)描述,(Si,Ei,Pi)表示任务从第Si秒开始,在第Ei秒后结束(第Si秒和Ei秒任务也在运行
),其优先级为Pi。同一时间可能有多个任务同时执行,它们的优先级可能相同,也可能不同。调度系统会经常向
查询系统询问,第Xi秒正在运行的任务中,优先级最小的Ki个任务(即将任务按照优先级从小到大排序后取前Ki个
)的优先级之和是多少。特别的,如果Ki大于第Xi秒正在运行的任务总数,则直接回答第Xi秒正在运行的任务优先
级之和。上述所有参数均为整数,时间的范围在1到n之间(包含1和n)。
Input
输入文件第一行包含两个空格分开的正整数m和n,分别表示任务总数和时间范围。接下来m行,每行包含三个空格
分开的正整数Si、Ei和Pi(Si≤Ei),描述一个任务。接下来n行,每行包含四个空格分开的整数Xi、Ai、Bi和Ci,
描述一次查询。查询的参数Ki需要由公式 Ki=1+(Ai*Pre+Bi) mod Ci计算得到。其中Pre表示上一次查询的结果,
对于第一次查询,Pre=1。
Output
输出共n行,每行一个整数,表示查询结果。
Sample Input
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
Sample Output
2 8 11
HINT
样例解释
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的一个排列
/* 做完这道题,感觉自己对主席树的理解还是蒙蔽的 很显然,我们需要对于每个时间点维护一个线段树,修改时运用差分的原理,对于在x~y上修改就在x上+1,在y+1上-1, 这样就完美契合了主席树。 对于每棵线段树,从小到大维护优先级即可。 */#include#include #include #include #include #define N 300010#define lon long longusing namespace std;lon n,m,lim,cnt,cn,to[N*2],root[N*2],son[N*40][2],w[N*40],sum[N*40];struct node{ lon pos,v;};node e[N*2];lon read(){ lon num=0,flag=1;char c=getchar(); while(c<'0'||c>'9'){ if(c=='-')flag=-1;c=getchar();} while(c>='0'&&c<='9'){num=num*10+c-'0';c=getchar();} return num*flag;}bool cmp(const node&s1,const node&s2){ return s1.pos 0) w[y]=w[x]+1; else w[y]=w[x]-1; sum[y]=sum[x]+v; return; } son[y][0]=son[x][0]; son[y][1]=son[x][1]; lon mid=l+r>>1; if(abs(v)<=mid) insert(son[x][0],son[y][0],l,mid,v); else insert(son[x][1],son[y][1],mid+1,r,v); push_up(y);}lon query(lon R,lon k){ lon x=root[R]; if(w[x]<=k)return sum[x]; lon l=1,r=lim,ans=0; while(l >1; if(w[son[x][0]]>=k){ r=mid; x=son[x][0]; } else { ans+=sum[son[x][0]]; k-=w[son[x][0]]; l=mid+1; x=son[x][1]; } if(l>=r) return ans+sum[x]/w[x]*k; } return ans;}int main(){ n=read();m=read(); for(lon i=1;i<=n;i++){ lon x=read(),y=read(),p=read(); e[++cnt].pos=x;e[cnt].v=p; e[++cnt].pos=y+1;e[cnt].v=-p; lim=max(lim,p); } sort(e+1,e+cnt+1,cmp); for(lon i=1;i<=cnt;i++) insert(root[i-1],root[i],1,lim,e[i].v); for(lon i=cnt;i>=1;i--) if(e[i].pos!=e[i+1].pos) to[e[i].pos]=i; for(lon i=1;i<=n;i++) if(!to[i])to[i]=to[i-1]; lon pre=1; for(lon i=1;i<=m;i++){ lon x=read(),a=read(),b=read(),c=read(); lon k=(a*pre+b)%c+1; pre=query(to[x],k); printf("%lld\n",pre); } return 0;}