BZOJ3932（主席树上二分+差分

按时间作为主席树的版本，每个版本的主席树都是一个权值线段树。

差分消去时间影响

对于当前时间版本的主席树查询前K大即可。

树上二分时结束后切记判定l==r的状态（易错

l==r叶子节点可能存在多个值（值大小为sum/siz ）

用I64dOLE了好久 。。 。。。。 。 。。 。 。 。 。 。 。 。 。。 。

用bit/stdc++.h  CE.....

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N = ;
const int LOGN = ;
const int maxn = N*LOGN;
#define lson son[rt][0]
#define rson son[rt][1]
#define un(x) (x>0?x:-x)
struct node{
    int pos,val;
};
node cf[N<<];
ll sum[maxn],siz[maxn];
int son[maxn][],root[N],id[N];
void pushup(int rt)
{
    siz[rt] = siz[lson]+siz[rson];
    sum[rt] = sum[lson]+sum[rson];
}
int tot = ;
void update(int x,int& y,int l,int r,int p)
{
    y = ++tot;
    if(l==r)
    {
        if(p>) siz[y] = siz[x]+;
        else siz[y] = siz[x]-;
        sum[y] = sum[x]+p;
        return ;
    }
    son[y][] = son[x][];
    son[y][] = son[x][];
    int m = (r+l)>>;
    if(un(p)<=m) update(son[x][],son[y][],l,m,p);
    else update(son[x][],son[y][],m+,r,p);
    pushup(y);

}
ll  bfind(int l,int r,int rt,int k)
{
    //cout<<"size:: "<<siz[rt]<<sum[rt]<<endl;
    if(siz[rt]<=k) return sum[rt];
    ll ans = ;
    while(l<r)
    {
        int m = (r+l)>>;
        if(k<siz[lson])
        {
            r = m;
            rt = lson;
        }
        else{
            l = m+;
            k -= siz[lson];
            ans+=sum[lson];
            rt = rson;
        }
    }
    if(k&&siz[rt]>) ans+=sum[rt]/siz[rt]*min(k,siz[rt]);
    return ans;
}
bool cmp(node a,node b)
{
    return a.pos<b.pos;
}
int main()
{
    int n,m;
    scanf("%d%d",&m,&n);
    int s,e,p;
    int lim = ;
    int cnt = ;
    for(int i=;i<m;i++){
        scanf("%d%d%d",&s,&e,&p);
        cf[++cnt] = (node){s,p};
        cf[++cnt] = (node){e+,-p};
        lim = max(lim,p);
    }
    sort(cf+,cf+cnt+,cmp);
    ll pre = ;
    for(int i=;i<=cnt;i++)update(root[i-],root[i],,lim,cf[i].val);
    for(int i=cnt;i>=;i--){
        if(cf[i].pos!=cf[i+].pos) id[cf[i].pos] = i;
    }
    for(int i=;i<=n;i++){
        if(id[i]==)id[i] = id[i-];
        //cout<<id[i]<<endl;
    }
    int x,a,b,c,k;
    for(int i=;i<n;i++)
    {
        scanf("%d%d%d%d",&x,&a,&b,&c);
        k = (a*pre+b)%c+;
        pre = bfind(,lim,root[id[x]],k);
        //cout<<x<<"  "<<id[x]<<" ";
        printf("%lld\n",pre);
    }
    return ;
}

AC代码