bzoj3932: [CQOI2015]任务查询系统 主席树

niabby讲过，不过niabby之后就退役了。。。。。。这个题离线后就是裸的主席树了(连删除操作都没有。。。)。然后怎么搞还要我说？空间要开这么大我也是醉了

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long sint;
#define ss printf("orz\n");
#define s2 printf("*******************************************\n");
#define maxn 5001000
#define maxm 101000
int Min(int a,int b)
{
    return a>b?b:a;
}
int Max(int a,int b)
{
    return a>b?a:b;
}
int getint()
{
    char c;int res;
    while(c=getchar(),c<'0'||c>'9');
    res=c-'0';
    while(c=getchar(),c>='0'&&c<='9')
        res=res*10+c-'0';
    return res;
}
int tot,ls[maxn],rs[maxn],siz[maxn];
int root[maxm],n,m,maxl;
sint sum[maxn],lastans=1;
struct node
{
    int f,tim,val;
}sa[maxm*2];
bool cmp(node aa,node bb)
{
    return aa.tim<bb.tim;
}
void update(int l,int r,int x,int &y,int flag,int val)
{
    //printf("%d %d %d %d %d\n",l,r,x,flag,val);
    //system("pause");
    if(!y)
    {
        y=++tot;
        sum[y]=sum[x];
        siz[y]=siz[x];
        ls[y]=ls[x];rs[y]=rs[x];
    }
    sum[y]+=flag*val;
    siz[y]+=flag;
    if(l==r) return;
    int mid=(l+r)>>1;
    if(val<=mid)
    {
        if(ls[y]==ls[x]) ls[y]=0;
        update(l,mid,ls[x],ls[y],flag,val);
    }
    else
    {
        if(rs[y]==rs[x]) rs[y]=0;
        update(mid+1,r,rs[x],rs[y],flag,val);
    }
}
sint getans(int l,int r,int rt,int k)
{
    //printf("%d %d %d %d\n",l,r,k,siz[rt]);
    //system("pause");
    int mid=l+r>>1;
    if(l==r)
        return (long long)mid*Min(k,siz[rt]);
    if(k<=siz[ls[rt]])
        return getans(l,mid,ls[rt],k);
    else
        return getans(mid+1,r,rs[rt],k-siz[ls[rt]])+sum[ls[rt]] ;
}
 
void print(int l,int r,int rt)
{
    printf("%d %d %d %d %d",l,r,rt,siz[rt],sum[rt]);
    //system("pause");
    int mid=(l+r)>>1;
    if(ls[rt])
    print(l,mid,ls[rt]);
    if(rs[rt])
    print(mid+1,r,rs[rt]);
}
int main()
{
    int a,b,c,d,x,y,z;
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++)
    {
        x=getint();y=getint();z=getint();
        sa[i].f=1;sa[i+m].f=-1;
        sa[i].tim=x;sa[i+m].tim=y+1;
        sa[i].val=z;sa[i+m].val=z;
        maxl=Max(maxl,z);
    }
    sort(sa+1,sa+m+m+1,cmp);
 
    /*s2
    for(int i=1;i<=m*2;i++)
    {
        printf("%d %d %d\n",sa[i].tim,sa[i].f,sa[i].val);
    }
    ss
    */
    for(int j=1,i=1;i<=n;i++)
    {
        //s2
        int flag=0;
        root[i]=root[i-1];
        for(;j<=m<<1&&sa[j].tim==i;j++)
        {
         //   s2
            if(flag==0&&sa[j].tim==i){flag=1;root[i]=0;}
            update(1,maxl,root[i-1],root[i],sa[j].f,sa[j].val);
        }
    }
 
    for(int i=1;i<=n;i++)
    {
        a=getint();b=getint();
        c=getint();d=getint();
        int k=(b*lastans+c)%d+1;
        printf("%lld\n",lastans=getans(1,maxl,root[a],k));
    }
    return 0;
}