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\n",pre); } return 0; }