BZOJ 3932

最新推荐文章于 2019-05-13 08:18:00 发布

weixin_30765319

最新推荐文章于 2019-05-13 08:18:00 发布

阅读量41

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/Stump/p/7944492.html

http://www.lydsy.com/JudgeOnline/problem.php?id=3932

主席树棵题

按时间每次建一棵线段树

然后合并成主席树

考虑差分的性质

对于三元组(S,E,P)

每次修改变为

在S时间时，P处+1

在T+1时间时，P处-1

然后对应的区间也做相同的操作

根据修改操作建树

#include<cstdio>
#include<algorithm>
#define FOR(i,s,t) for(register int i=s;i<=t;++i)
using namespace std;
typedef long long ll;
const int NlogN=8000011;
const int N=400011;
struct node{
    int ls,rs;
    int sum;
    ll res;
}tr[NlogN];
int n,m;
int a[N>>1];
int root[N>>1];
struct cg{
    int t,p,v;
    inline bool operator<(cg A)const{
        if(t!=A.t)return t<A.t;
        if(p!=A.p)return p<A.p;
        return v>A.v;
    }
}c[N];
int s,t,p,tot,cnt,cpy;
inline void disc_init(){
    sort(a+1,a+a[0]+1);
    a[0]=unique(a+1,a+a[0]+1)-a-1;
    FOR(i,1,(n<<1))c[i].p=lower_bound(a+1,a+a[0]+1,c[i].p)-a;
}
inline void change(int p,int &now,int l,int r,int pos,int v){
    tr[now=++cnt]=tr[p];
    tr[now].sum+=v;tr[now].res+=1ll*v*a[pos];
    if(l==r)return;
    int mid=(l+r)>>1;
    pos<=mid?change(tr[p].ls,tr[now].ls,l,mid,pos,v):change(tr[p].rs,tr[now].rs,mid+1,r,pos,v);
}
inline ll query(int now,int l,int r,int k){
    if(!now)return 0;
    if(l==r)return !tr[now].sum?0:1ll*tr[now].res/tr[now].sum*k;
    int mid=(l+r)>>1;
    int data=tr[tr[now].ls].sum;
    if(k<=data)return query(tr[now].ls,l,mid,k);
    return tr[tr[now].ls].res+query(tr[now].rs,mid+1,r,k-data);
}
ll pre=1;
int A,B,C,x,k;
int main(){
    scanf("%d%d",&n,&m);
    FOR(i,1,n){
        scanf("%d%d%d",&s,&t,&p);
        c[++tot]=(cg){s,p,1};c[++tot]=(cg){t+1,p,-1};
        a[++a[0]]=p;
    }
    disc_init();
    sort(c+1,c+(n<<1|1));
    tot=1;
    FOR(i,1,100000){
        root[i]=root[i-1];
        while(c[tot].t==i&&tot<=(n<<1)){
            cpy=root[i];
            change(cpy,root[i],1,a[0],c[tot].p,c[tot].v);
            ++tot;
        }
    }
    while(m--){
        scanf("%d%d%d%d",&x,&A,&B,&C);
        A%=C;B%=C;
        k=(A*pre%C+B)%C+1;
        k=min(k,tr[root[x]].sum);
        pre=query(root[x],1,a[0],k);
        printf("%lld\n",pre);
    }
    return 0;
}

转载于:https://www.cnblogs.com/Stump/p/7944492.html