bzoj 3932 [CQOI2015]任务查询系统可持久化二维线段树

最新推荐文章于 2019-03-07 08:46:23 发布

原创最新推荐文章于 2019-03-07 08:46:23 发布 · 415 阅读

0 ·

CC 4.0 BY-SA版权

树套树同时被 2 个专栏收录

9 篇文章

订阅专栏

主席树

8 篇文章

订阅专栏

本文介绍了一种使用主席树解决区间查询问题的方法，并通过实际案例分享了如何避免内存溢出和超时错误。作者对比了两种不同的实现思路，一种是采用可持久化二维线段树，另一种则是更简洁高效的主席树方案。

好好的一道主席树被我强行套了一层。。。
然后MLE+TLE，然后我惊奇地发现区间的sum只要开int就行了。。。

我的做法：
一个点需要满足si<=x且ei>=x，然后求满足这个条件的pi前k小的点的和。
那么对si排一遍序，把所有点按ei第一维，pi第二维插入可持久化二维线段树。。。。

正常的做法：
对于一个点，在si时+1，ti+1时-1，然后按这个建主席树。

#include <bits/stdc++.h>
using namespace std;
#define N 110000
#define M 1800000
#define A 31300000
#define ll long long
#define ls l,mid,ch[now][0]
#define rs mid+1,r,ch[now][1]
int n,m,v1,v2,K;
int st[N],root[N],top,pos[N];
int st1[1100],top1;
ll ans;
char getc()
{
    static const int LEN = 4096;
    static char buf[LEN],*S=buf,*T=buf;
    if(S == T)
    {
        T = (S=buf)+fread(buf,1,LEN,stdin);
        if(S == T)return EOF;
    }
    return *S++;
}
int read()
{
    static char ch;
    static int D,type;
    type=0;
    while(!isdigit(ch=getc()))type|=(ch=='-');
    for(D=ch-'0'; isdigit(ch=getc());)
        D=(D<<3)+(D<<1)+(ch-'0');
    return type ? -D:D;
}
struct node
{
    int s,e,p;
    void red(){s=read();e=read();p=read();st[++top]=p;}
    friend bool operator < (const node &r1,const node &r2)
    {return r1.s<r2.s;};
}a[N];
struct seg_tree2
{
    int ch[A][2],num[A],cnt;
    int sum[A];
    void insert(int l,int r,int &now,int pre)
    {
        sum[now=++cnt]=sum[pre]+st[v2];
        num[now]=num[pre]+1;
        if(l==r)return;
        int mid=(l+r)>>1;
        ch[now][0]=ch[pre][0];
        ch[now][1]=ch[pre][1];
        if(mid>=v2)insert(ls,ch[pre][0]);
        else insert(rs,ch[pre][1]);
    }
    void trs(int x)
    {
        for(int i=1;i<=top1;i++)
            st1[i]=ch[st1[i]][x];
    }
    ll query(int l,int r,int rem)
    {
        ll nm=0,sm=0;
        if(l==r)
        {
            for(int i=1;i<=top1;i++)
                nm+=num[st1[i]];        
            return (ll)min((ll)rem,nm)*st[l];
        }
        for(int i=1;i<=top1;i++)
        {
            nm+=num[ch[st1[i]][0]];
            sm+=sum[ch[st1[i]][0]];
        }
        int mid=(l+r)>>1;
        if(nm>=rem){trs(0);return query(l,mid,rem);}
        else {trs(1);return sm+query(mid+1,r,rem-nm);}
    }
}tr2;
struct seg_tree1
{
    int ch[M][2],root[M],cnt;
    void insert(int l,int r,int &now,int pre)
    {
        now=++cnt;
        tr2.insert(1,top,root[now],root[pre]);
        if(l==r)return;
        int mid=(l+r)>>1;
        ch[now][0]=ch[pre][0];
        ch[now][1]=ch[pre][1];
        if(mid>=v1)insert(ls,ch[pre][0]);
        else insert(rs,ch[pre][1]);
    }
    void query(int l,int r,int now)
    {
        if(l>=v1)
            {st1[++top1]=root[now];return;}
        int mid=(l+r)>>1;
        query(rs);
        if(mid>=v1)query(ls);
    }
    ll query(int x)
    {
        top1=0;query(1,m,x);
        return tr2.query(1,top,K);
    }
}tr1;
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)a[i].red();
    sort(st+1,st+1+top);
    top=unique(st+1,st+1+top)-st-1;
    for(int i=1;i<=n;i++)
        a[i].p=lower_bound(st+1,st+1+top,a[i].p)-st;
    sort(a+1,a+1+n);
    for(int i=1;i<=n;i++)
    {
        v1=a[i].e;v2=a[i].p;
        tr1.insert(1,m,root[i],root[i-1]);
    }
    for(int i=1,j=0;i<=m;i++)
    {
        while(j<n&&a[j+1].s<=i)j++;
        pos[i]=j;
    }
    ans=1;
    for(int i=1,x,a,b,c;i<=m;i++)
    {
        x=read();a=read();b=read();c=read();
        K=1+(a*(ans%c)+b)%c;v1=x;
        printf("%lld\n",ans=tr1.query(root[pos[x]]));
    }
    return 0;
}