BOJ)14577 일기예보

문제:icpc.me/14577

4가지 종류의 질의가 들어오는 문제이다.

지역 수 가 N 질의의 수가 M일 때 N,M이 10만까지 들어오므로 적어도 질의 하나당 log 시간에 해결해줘야 한다.

결국 3,4번 질의를 로그시간에 답해줘야 하는데 구간에 대한 질의이므로 우선 세그먼트 트리를 생각해낼 수 있다.

하지만 좌표 자체에 업데이트를 할 경우 좌표의 범위부터가 감당이 안된다. 

그렇기 때문에 우리는 한 지역에 쌓일 수 있는 눈의 양에 업데이트를 해야한다. 

나올 수 있는 눈의 양은 3번 쿼리의 경우 2개 1,2번 쿼리의 경우 1개 일 뿐이니 최대 2*M개 뿐이기 때문이다. 

말을 어렵게 표현했는데 우선 쿼리를 오프라인으로 처리할거다.

모든 쿼리를 받은 뒤에 해당 질의 때 마다 나올 수 있는 모든 눈의 양을 저장해 준 뒤 

크기 순서대로 정렬하여 번호를 매겨준다.

이제 1,2번의 경우 해당 눈의 양에 업데이트를 시켜준다면 

3번 질의는 구간의 합을 출력하는 쿼리로 4번 질의는 이분 탐색 쿼리로 변형된다.

#include <cstdio>
#include <algorithm>
#include <map>
#include <set>
using namespace std;
typedef long long ll;
ll seg[1000400];
ll update(ll pos, ll val, ll node, ll x, ll y) {
    if (y < pos || pos < x)return seg[node];
    if (x == y)return seg[node] += val;
    ll mid = (x + y) >> 1;
    return seg[node] = update(pos, val, node * 2, x, mid) + update(pos, val, node * 2 + 1, mid + 1, y);
}
ll query(ll lo, ll hi, ll node, ll x, ll y) {
    if (y < lo || hi < x)return 0;
    if (lo<=x&&y<=hi)return seg[node];
    ll mid = (x + y) >> 1;
    return query(lo, hi, node * 2, x, mid) + query(lo, hi, node * 2 + 1, mid + 1, y);
}
ll bquery(ll val, ll node, ll x, ll y) {
    if (x == y)return x;
    ll mid = (x + y) >> 1;
    if (seg[node * 2] >= val) 
        return bquery(val, node * 2, x, mid);
    return bquery(val - seg[node * 2], node * 2 + 1, mid + 1, y);
}
ll q[100010][3], a[100010], b[100010], n, m, rmp[400040];
set<ll> st;
map<ll, ll> mp;
int main() {
    scanf("%lld%lld", &n, &m);
    for (int i = 1; i <= n; i++) {
        scanf("%lld", &a[i]);
        b[i] = a[i];
        st.insert(a[i]);
    }
    for (int i = 0; i < m; i++) {
        scanf("%lld", &q[i][0]);
        if (q[i][0] == 4) 
            scanf("%lld", &q[i][1]);
        else {
            scanf("%lld%lld", &q[i][1], &q[i][2]);
            if (q[i][0] == 1) {
                a[q[i][1]] += q[i][2];
                st.insert(a[q[i][1]]);
            }
            else if (q[i][0] == 2) {
                a[q[i][1]] -= q[i][2];
                st.insert(a[q[i][1]]);
            }
            else {
                st.insert(q[i][1]);
                st.insert(q[i][2]);
            }
        }
    }
    ll cnt = 0;
    for (auto next : st) {
        mp[next] = ++cnt;
        rmp[cnt] = next;
    }
    for (int i = 1; i <= n; i++) 
        update(mp[b[i]], 1, 1, 1, cnt);
    for (int i = 0; i < m; i++) {
        if (q[i][0] == 1) {
            ll idx = q[i][1];
            update(mp[b[idx]], -1, 1, 1, cnt);
            b[idx] += q[i][2];
            update(mp[b[idx]], 1, 1, 1, cnt);
        }
        else if (q[i][0] == 2) {
            ll idx = q[i][1];
            update(mp[b[idx]], -1, 1, 1, cnt);
            b[idx] -= q[i][2];
            update(mp[b[idx]], 1, 1, 1, cnt);
        }
        else if (q[i][0] == 3) {
            printf("%lld\n", query(mp[q[i][1]], mp[q[i][2]], 1, 1, cnt));
        }
        else {
            ll f = n + 1 - q[i][1];
            printf("%lld\n", rmp[bquery(f, 1, 1, cnt)]);
        }
    }
    return 0;
}