唉，一道主席树板子题竟然WA得怀疑人生。。。
首先任务是[s,t]这一段出现，我们可以选择区间加法也可以选择直接差分，显然差分好打得多。。。
也就是说把询问拆成两个，一个是在s加，另一个是在t+1减。。。把这些事件按时间排序。。。
我们对于每一个时间点建立一颗权值线段树，然后先复制前一个时间点的历史版本。。。
然后一个一个的把该时间点的事件加入，该点的权值线段树中，此时以自己为历史版本。。。
然后query的就是一个很正常的一个权值线段树上的区间求和，然而这样会获得80分。。。
query 的递归边界(l==r)时，一个权值上可能有多个任务，就是说可能大于query传上来的K。。
所以我们要return sum[x]/size[x]*K。。。
哦，忘了说要离散化一下。。。
// MADE BY QT666
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<cstring>
#define int long long
using namespace std;
typedef long long ll;
const int N=100050;
int root[N*2],ls[N*40],rs[N*40],sum[N*40],size[N*40],sz;
int hsh[N],tot,m,n,cnt;
struct data{
int t,val,type;
}q[N*2];
bool cmp(const data &a,const data &b){
return a.t<b.t;
}
void insert(int &y,int x,int l,int r,int v,int type){
y=++sz;size[y]=size[x]+type;
ls[y]=ls[x];rs[y]=rs[x];
sum[y]=sum[x]+hsh[v]*type;
if(l==r) return;
int mid=(l+r)>>1;
if(v<=mid) insert(ls[y],ls[x],l,mid,v,type);
else insert(rs[y],rs[x],mid+1,r,v,type);
}
int query(int x,int l,int r,int K){
if(l==r){
if(size[x]) return sum[x]/size[x]*K;
else return 0;
}
int mid=(l+r)>>1;
if(size[ls[x]]>=K) return query(ls[x],l,mid,K);
else return sum[ls[x]]+query(rs[x],mid+1,r,K-size[ls[x]]);
}
main(){
scanf("%lld%lld",&m,&n);
for(int i=1;i<=m;i++){
int st,ed,p;scanf("%lld%lld%lld",&st,&ed,&p);
q[++cnt]=(data){st,p,1};q[++cnt]=(data){ed+1,p,-1};
hsh[++tot]=p;
}
sort(hsh+1,hsh+1+tot);tot=unique(hsh+1,hsh+tot+1)-hsh-1;
sort(q+1,q+1+cnt,cmp);
for(int i=1;i<=cnt;i++) q[i].val=lower_bound(hsh+1,hsh+1+tot,q[i].val)-hsh;
int j=1;
for(int i=1;i<=n;i++){
root[i]=root[i-1];
while(q[j].t<=i&&j<=cnt) insert(root[i],root[i],1,m,q[j].val,q[j].type),j++;
}
int pre=1;
for(int i=1;i<=n;i++){
int x,a,b,c;scanf("%lld%lld%lld%lld",&x,&a,&b,&c);
int k=1+(a*pre+b)%c;
k=min(k,size[root[x]]);pre=query(root[x],1,m,k);
printf("%lld\n",pre);
}
return 0;
}
原文：http://www.cnblogs.com/qt666/p/7231552.html