bzoj 3932: [CQOI2015]任务查询系统

Description

最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。超级计算机中的

任务用三元组(Si,Ei,Pi)描述，(Si,Ei,Pi)表示任务从第Si秒开始，在第Ei秒后结束（第Si秒和Ei秒任务也在运行

），其优先级为Pi。同一时间可能有多个任务同时执行，它们的优先级可能相同，也可能不同。调度系统会经常向

查询系统询问，第Xi秒正在运行的任务中，优先级最小的Ki个任务（即将任务按照优先级从小到大排序后取前Ki个

）的优先级之和是多少。特别的，如果Ki大于第Xi秒正在运行的任务总数，则直接回答第Xi秒正在运行的任务优先

级之和。上述所有参数均为整数，时间的范围在1到n之间（包含1和n）。

Input

输入文件第一行包含两个空格分开的正整数m和n，分别表示任务总数和时间范围。接下来m行，每行包含三个空格

分开的正整数Si、Ei和Pi(Si≤Ei)，描述一个任务。接下来n行，每行包含四个空格分开的整数Xi、Ai、Bi和Ci，

描述一次查询。查询的参数Ki需要由公式 Ki=1+(Ai*Pre+Bi) mod Ci计算得到。其中Pre表示上一次查询的结果，

对于第一次查询，Pre=1。

  

Output

输出共n行，每行一个整数，表示查询结果。

 

Sample Input

4 3
1 2 6
2 3 3
1 3 2
3 3 4
3 1 3 2
1 1 3 4
2 2 4 3

Sample Output

2
8
11

HINT

样例解释

K1 = (1*1+3)%2+1 = 1

K2 = (1*2+3)%4+1 = 2

K3 = (2*8+4)%3+1 = 3

对于100%的数据，1≤m,n,Si,Ei,Ci≤100000，0≤Ai,Bi≤100000，1≤Pi≤10000000，Xi为1到n的一个排列

 

Source

唉，一道主席树板子题竟然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
",pre);
   }
  return 0;
}