5308: [Zjoi2018]胖

5308: [Zjoi2018]胖

链接

分析：

　　题目转化为一个点可以更新多少个点，一个点可以更新的点一定是一个区间，考虑二分左右端点确定这个区间。

　　设当前点是x，向右二分一个点y，如果x可以更新到y，那么在x~y之间的所有关键点（存在宫殿往这边的点）到y的距离小于x到y的距离，以及y~2*y-x之间的关键点到y的距离小于x到y的距离。

　　当然一个点可能同时被两个点更新，那么让最靠左的点更新它。

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cmath>
#include<cctype>
#include<set>
#include<queue>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;

inline int read() {
    int x=,f=;char ch=getchar();for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-;
    for(;isdigit(ch);ch=getchar())x=x*+ch-'';return x*f;
}

const int N = ;
int Log[N], n, m, k;
LL dis[N];
struct Node {
    int p; LL l;
    Node() {}
    Node(int _p,LL _l) { p = _p, l = _l; }
    bool operator < (const Node &A) const { return p < A.p; }
}a[N];
struct ST{
    LL f[][N];
    void init() {
        for (int i = ; i <= k; ++i) f[][i] = a[i].l;
        for (int j = ; j <= Log[k]; ++j)
            for (int i = ; i + ( << j) -  <= k; ++i)
                f[j][i] = min(f[j - ][i], f[j - ][i + ( << (j - ))]);
    }
    LL query(int l,int r) {
        l = max(, l), r = min(r, n);
        l = lower_bound(a + , a + k + , Node(l, )) - a;
        r = upper_bound(a + , a + k + , Node(r, )) - a - ;
        if (l > r) return 1e18;
        int k = Log[r - l + ];
        return min(f[k][l], f[k][r - ( << k) + ]);
    }
}L, R;

bool check1(int x,int y) { // [2 * y - x + 1, y, x]
    if (x == y) return true;
    LL a = L.query( * y - x + , y) + dis[y];
    LL b = R.query(y, x - ) - dis[y];
    LL c = R.query(x, x) - dis[y];
    if (a <= c || b <= c) return false;
    if ( * y - x >= ) return L.query( * y - x,  * y - x) + dis[y] > c; // 如果距离相同，优先给左边的点
    return true;
}
int findL(int x) {
    int l = , r = x, ans = ;
    while (l <= r) {
        int mid = (l + r) >> ;
        if (check1(x, mid)) r = mid - , ans = mid;
        else l = mid + ;
    }
    return ans;
}
bool check2(int x,int y) { // [x, y, 2 * y - x - 1]
    if (x == y) return true;
    LL a = L.query(x + , y) + dis[y];
    LL b = R.query(y,  * y - x - ) - dis[y];
    LL c = L.query(x, x) + dis[y];
    if (a <= c || b <= c) return false;
    if ( * y - x <= n) return R.query( * y - x,  * y - x) - dis[y] >= c;
    return true;
}
int findR(int x) {
    int l = x, r = n, ans = ;
    while (l <= r) {
        int mid = (l + r) >> ;
        if (check2(x, mid)) l = mid + , ans = mid;
        else r = mid - ;
    }
    return ans;
}
int main() {
    n = read(), m = read();
    for (int i = ; i <= n; ++i) Log[i] = Log[i >> ] + ;
    for (int i = ; i <= n; ++i) dis[i] = dis[i - ] + read();
    while (m --) {
        LL ans = ;
        k = read();
        for (int i = ; i <= k; ++i) a[i].p = read(), a[i].l = read();
        sort(a + , a + k + );
        for (int i = ; i <= k; ++i) a[i].l -= dis[a[i].p]; L.init();
        for (int i = ; i <= k; ++i) a[i].l += dis[a[i].p] * ; R.init();
        for (int i = ; i <= k; ++i) ans += findR(a[i].p) - findL(a[i].p) + ;
        printf("%lld\n", ans);
    }
    return ;
}