bzoj 3932: [CQOI2015]任务查询系统 可持久化纤段树

题意

最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。超级计算机中的
任务用三元组(Si,Ei,Pi)描述，(Si,Ei,Pi)表示任务从第Si秒开始，在第Ei秒后结束（第Si秒和Ei秒任务也在运行
），其优先级为Pi。同一时间可能有多个任务同时执行，它们的优先级可能相同，也可能不同。调度系统会经常向
查询系统询问，第Xi秒正在运行的任务中，优先级最小的Ki个任务（即将任务按照优先级从小到大排序后取前Ki个
）的优先级之和是多少。特别的，如果Ki大于第Xi秒正在运行的任务总数，则直接回答第Xi秒正在运行的任务优先
级之和。上述所有参数均为整数，时间的范围在1到n之间（包含1和n）。

分析

对操作查分，然后按照时间建主席树即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define N 100005
#define LL long long
using namespace std;

int n,m,cnt,tot,last[N],root[N];
struct edge{int val,p,next;}e[N*2];
struct tree{int l,r,num;LL val;}t[N*50];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void ins(int &d,int p,int l,int r,int x,int y)
{
    d=++tot;
    t[d].l=t[p].l;t[d].r=t[p].r;t[d].num=t[p].num;t[d].val=t[p].val;
    t[d].val+=x*y;t[d].num+=y;
    if (l==r) return;
    int mid=(l+r)/2;
    if (x<=mid) ins(t[d].l,t[p].l,l,mid,x,y);
    else ins(t[d].r,t[p].r,mid+1,r,x,y);
}

LL query(int d,int l,int r,int k)
{
    if (l==r) return (LL)k*l;
    if (t[d].num<=k) return t[d].val;
    int mid=(l+r)/2;
    if (t[t[d].l].num<k) return t[t[d].l].val+query(t[d].r,mid+1,r,k-t[t[d].l].num);
    else return query(t[d].l,l,mid,k);
}

int main()
{
    m=read();n=read();
    for (int i=1;i<=m;i++)
    {
        int s=read(),E=read(),p=read();
        e[++cnt].val=1;e[cnt].p=p;e[cnt].next=last[s];last[s]=cnt;
        e[++cnt].val=-1;e[cnt].p=p;e[cnt].next=last[E+1];last[E+1]=cnt;
    }
    for (int i=1;i<=n;i++)
    {
        root[i]=root[i-1];
        for (int j=last[i];j;j=e[j].next)
            ins(root[i],root[i],1,10000000,e[j].p,e[j].val);
    }
    LL pre=1;
    for (int i=1;i<=n;i++)
    {
        int k=read(),a=read(),b=read(),c=read();
        pre=query(root[k],1,10000000,1+(LL)((LL)a*pre+(LL)b)%c);
        printf("%lld\n",pre);
    }
    return 0;
}