節點 pen tps real 系統 最小 CQ line pre

Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 4509 Solved: 1497
[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

很水但是很惡心

首先不難想到對兩個端點進行差分

然後把主席樹建出來就可以了

註意要開ll

另外，葉子節點的總權值可能大於查詢的個數，需要特判

#include<cstdio>
#include<algorithm>
#include<vector>
#define int long long 
using namespace std;
const int MAXN = 1e6 + 10;
inline int read() {
    char c = getchar(); int x = 0, f = 1; 
    while(c < ‘0‘ || c > ‘9‘) {if(c == ‘-‘) f = -1; c = getchar();}
    while(c >= ‘0‘ && c <= ‘9‘) x = x * 10 + c - ‘0‘, c = getchar();
    return x * f;
}
int N, M, PRE = 1;
struct node {
    int ls, rs, sum, siz;
}T[MAXN * 6];
int root[MAXN], tot = 0, realrt[MAXN];
int S[MAXN], E[MAXN], tmp[MAXN], P[MAXN];
vector<int>v[MAXN];
#define rs(x) T[x].rs
#define ls(x) T[x].ls
int insert(int &now, int pre, int l, int r, int pos, int val) {
    now = ++tot;
    T[now] = T[pre]; T[now].sum += val;
    val > 0 ? T[now].siz++ : T[now].siz--;
    if(l == r) return now;
    int mid = l + r >> 1;
    if(pos <= mid) ls(now) = insert(ls(now), ls(pre), l, mid, pos, val);
    else            rs(now) = insert(rs(now), rs(pre), mid + 1, r, pos, val);
    return now;
}
int Query(int now, int l, int r, int val) {
    if(l == r)  return min(val * tmp[l], T[now].sum);
    int mid = l + r >> 1;
    if(val <= T[ls(now)].siz) return Query(ls(now), l, mid, val);
    else return Query(rs(now), mid + 1, r, val - T[ls(now)].siz) + T[ls(now)].sum;
}
main() {
    #ifdef WIN32
    freopen("a.in", "r", stdin);
    freopen("a.out", "w", stdout);
    #endif
    N = read(); M = read();
    for(int i = 1; i <= N; i++) {
        S[i] = read(), E[i] = read(), tmp[i] = P[i] = read();
    }
    sort(tmp + 1, tmp + N + 1);
    int num = unique(tmp + 1, tmp + N + 1) - tmp - 1;
    for(int i = 1; i <= N; i++) {
        P[i] = lower_bound(tmp + 1, tmp + num + 1, P[i]) - tmp;
        v[S[i]].push_back(P[i]); v[E[i] + 1].push_back(-P[i]);
    }
    int cur = 0;
    for(int i = 1; i <= 100000; i++) {
        for(int j = 0; j < v[i].size(); j++) 
            cur++, insert(root[cur], root[cur - 1], 1, num, abs(v[i][j]), v[i][j] < 0 ? -tmp[-v[i][j]] : tmp[v[i][j]]);
        realrt[i] = cur;
    }
    for(int i = 1; i <= M; i++) {
        int Xi = read(), Ai = read(), Bi = read(), Ci = read();
        int Ki = 1 + (Ai * PRE + Bi) % Ci;
        PRE = Query(root[realrt[Xi]], 1, num, Ki);
        printf("%lld\n", PRE);
    }
}

BZOJ 3932: [CQOI2015]任務查詢系統