Portal

Description

給出\(n(n\leq10^5)\)個任務，和總時間範圍\(m(m\leq10^5)\)。每個任務有開始/結束時間\(s_i,e_i(1\leq s_i \leq e_i \leq m)\)和優先級\(p_i(p_i\leq10^9)\)。接下來\(m\)個詢問，每次詢問在時刻\(t_i\)時優先級前\(k\)大的任務的優先級之和，若\(k\)大於此時正在進行的任務總數則輸出此時優先級之和。其中\(\{t_m\}\)是\(1..m\)的一個排列。

Solution

可持久化線段樹。
將優先級\(p\)離散化，然後對於每個時刻建立線段樹，記錄此時正在進行的任務情況，線段樹上第\(i\)

時間復雜度\(O(nlogn)\)。

Code

//「CQOI2015」任務查詢系統
#include <algorithm>
#include <cstdio>
using namespace std;
typedef long long lint;
inline char gc()
{
    static char now[1<<16],*s,*t;
    if(s==t) {t=(s=now)+fread(now,1,1<<16,stdin); if
 
(s==t) return EOF;}
    return *s++;
}
inline int read()
{
    int x=0; char ch=gc();
    while(ch<'0'||'9'<ch) ch=gc();
    while('0'<=ch&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x;
}
const int N=2e5+10;
int n,m;
struct optn{int opt,t,x;} seq[N];
bool
 
 cmpT(optn a,optn b) {return a.t<b.t;}
int n0,map[N];
void discrete()
{
    for(int i=1;i<=n+n;i++) map[i]=seq[i].x;
    sort(map+1,map+n+n+1); n0=unique(map+1,map+n+n+1)-map-1;
    for(int i=1;i<=n+n;i++) seq[i].x=lower_bound(map+1,map+n0+1,seq[i].x)-map;
}
const int N1=22*N;
int ndCnt,rt[N],ch[N1][2],cnt[N1] ;lint sum[N1];
void update(int p) {cnt[p]=cnt[ch[p][0]]+cnt[ch[p][1]],sum[p]=sum[ch[p][0]]+sum[ch[p][1]];}
void ndCopy(int p,int q) {ch[q][0]=ch[p][0],ch[q][1]=ch[p][1],cnt[q]=cnt[p],sum[q]=sum[p];}
void ins(int &p,int L0,int R0,int x,int v)
{
    ndCopy(p,++ndCnt); p=ndCnt;
    if(L0==R0) {cnt[p]+=v,sum[p]+=v*map[L0]; return;}
    int mid=L0+R0>>1;
    if(x<=mid) ins(ch[p][0],L0,mid,x,v);
    else ins(ch[p][1],mid+1,R0,x,v);
    update(p);
}
lint query(int p,int L0,int R0,int k)
{
    if(L0==R0) return 1LL*k*map[L0];
    if(cnt[p]<=k) return sum[p];
    int mid=L0+R0>>1;
    if(cnt[ch[p][0]]>=k) return query(ch[p][0],L0,mid,k);
    else return sum[ch[p][0]]+query(ch[p][1],mid+1,R0,k-cnt[ch[p][0]]);
}
int main()
{
    n=read(),m=read();
    for(int i=1;i<=n;i++)
    {
        int t1=read(),t2=read(),x=read();
        int i2=i<<1,i1=i2-1;
        seq[i1].opt=1,seq[i1].t=t1,seq[i1].x=x;
        seq[i2].opt=-1,seq[i2].t=t2+1,seq[i2].x=x;
    }
    discrete();
    sort(seq+1,seq+n+n+1,cmpT);
    for(int i=1,k=1;i<=1e5+10;i++)
    {
        rt[i]=rt[i-1];
        for(k;seq[k].t==i;k++) ins(rt[i],1,n0,seq[k].x,seq[k].opt);
    }
    lint pre=1;
    for(int i=1;i<=m;i++)
    {
        int t=read(); lint a=read(),b=read(),c=read();
        pre=query(rt[t],1,n0,1+(a*pre+b)%c);
        printf("%lld\n",pre);
    }
    return 0;
}

LibreOJ2097 - 「CQOI2015」任務查詢系統