bzoj3932

3932: [CQOI2015]任务查询系统

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 1326 Solved: 480
[Submit][Status][Discuss]

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

题解：同bzoj1901，写一个进来时加，出去时删。。。。。

  1 #include<set>
  2 #include<map>
  3 #include<ctime>
  4 #include<queue>
  5 #include<cmath>
  6 #include<cstdio>
  7 #include<vector>
  8 #include<cstring>
  9 #include<cstdlib>
 10 #include<iostream>
 11 #include<algorithm>
 12 #define inf 100000000000000LL
 13 #define pa pair<int,int>
 14 #define ll long long 
 15 using namespace std;
 16 ll read()
 17 {
 18     ll x=0,f=1;char ch=getchar();
 19     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
 20     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
 21     return x*f;
 22 }
 23 int m,n,cnt,mx;
 24 int q[100005],top;
 25 int S[100005],E[100005],P[100005],disc[100005];
 26 int root[100005],size[20000005],a[20000005][2];
 27 ll sum[20000005];
 28 int insert(int x,int val,int f)
 29 {
 30     if(!x)x=++cnt;
 31     size[x]+=f;sum[x]+=f*disc[val];
 32     int tmp=x;
 33     int l=1,r=2*m;
 34     while(l!=r)
 35     {
 36         int mid=(l+r)>>1,t=0;
 37         if(val>mid)t^=1,l=mid+1;
 38         else r=mid;
 39         if(!a[x][t])a[x][t]=++cnt;
 40         x=a[x][t];size[x]+=f;
 41         sum[x]+=f*disc[val];
 42     }
 43     return tmp;
 44 }
 45 void update(int x,int val,int f)
 46 {
 47     for(int i=x;i<=n;i+=i&-i)
 48         root[i]=insert(root[i],val,f);
 49 }
 50 ll query(int x,int K)
 51 {
 52     ll ans=0,tot=0;
 53     top=0;
 54     for(int i=x;i;i-=i&-i)
 55     {
 56         q[++top]=root[i],tot+=size[root[i]];
 57         ans+=sum[root[i]];
 58     }
 59     if(K>tot)return ans;
 60     else ans=0;
 61     int l=1,r=2*m;
 62     while(l!=r)
 63     {
 64         int mid=(l+r)>>1,t,tmp=0;
 65         for(int i=1;i<=top;i++)
 66             tmp+=size[a[q[i]][0]];
 67         if(tmp<K)
 68         {
 69             for(int i=1;i<=top;i++)ans+=sum[a[q[i]][0]];
 70             t=1;
 71             K-=tmp,l=mid+1;
 72         }
 73         else r=mid,t=0;
 74         for(int i=1;i<=top;i++)
 75             q[i]=a[q[i]][t];
 76     }
 77     ans+=1LL*K*disc[l];
 78     return ans;
 79 }
 80 int main()
 81 {
 82     m=read();n=read();
 83     for(int i=1;i<=m;i++)
 84     {
 85         S[i]=read(),E[i]=read(),P[i]=read();
 86         disc[i]=P[i];
 87     }
 88     sort(disc+1,disc+m+1);
 89     for(int i=1;i<=m;i++)
 90     {
 91         P[i]=lower_bound(disc+1,disc+m+1,P[i])-disc;
 92         update(S[i],P[i],1),update(E[i]+1,P[i],-1);
 93     }
 94     ll Pre=1,X,A,B,C,K;
 95     for(int i=1;i<=n;i++)
 96     {
 97         X=read();A=read();B=read();C=read();
 98         K=1+(A*Pre+B)%C;
 99         Pre=query(X,K);
100         printf("%lld
",Pre);
101     }
102     return 0;
103 }

View Code