bzoj3932
3932: [CQOI2015]任务查询系统
Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 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
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
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\n",Pre); 101 } 102 return 0; 103 }
我太蒟蒻了,所以神犇们留下意见让我跪膜