「ZJOI2018」胖（ST表+二分）

「ZJOI2018」胖（ST表+二分）

不开 \(O_2\) 又没卡过去是种怎么体验。。。

这可能是 \(ZJOI2018\) 最简单的一题了。。。我都能 \(A\)。。。

首先我们发现这个奇怪的图每个点扩展的是一个区间 \([L,R]\)，然后我们就可以二分端点了。

一个点 \(x\) 扩展到点 \(y\) 至少要 \(|x-y|\) 的时间，所以我们把 \(a_i\) 排个序，在上面二分一个合法的区间使得 \(|x-a_l|\leq t\) 且 \(|x-a_r|\leq t\)

然后若能扩展到 \(y\)，那么 \(0\) 号点到 \(y\) 号点的距离为 \(|dis_y-dis_x|+l\)。我们用两个 \(ST\) 表把绝对值拆掉，分别维护最小值即可。

时间复杂度 \(O(n\log^2 n)\)

为什么 \(dl\) 出题人会出到 \(2\times 10^5\)。。。两个 \(\log\) 一般只出到 \(10^5\) 的啊。。。

\(Code\ Below:\)

// luogu-judger-enable-o2
#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=200000+10;
int n,m,k,w[maxn],lg[maxn];ll dis[maxn];

struct node{
    int p;ll l;
    node(int p=0,ll l=0):p(p),l(l){}
}a[maxn];
inline bool operator < (const node &a,const node &b){
    return a.p<b.p;
}

struct Sparse_Table{
    ll st[maxn][18];
    inline void init(){
        for(int i=1;i<=k;i++) st[i][0]=a[i].l;
        for(int j=1;j<=lg[k];j++)
            for(int i=1;i+(1<<j)-1<=k;i++) st[i][j]=min(st[i][j-1],st[i+(1<<(j-1))][j-1]);
    }
    inline ll query(int l,int r){
        l=max(l,1);r=min(r,n);
        l=lower_bound(a+1,a+k+1,node(l))-a;
        r=upper_bound(a+1,a+k+1,node(r))-a-1;
        if(l>r) return 1e18;
        int k=lg[r-l+1];
        return min(st[l][k],st[r-(1<<k)+1][k]);
    }
}L,R;

inline int read(){
    register int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return (f==1)?x:-x;
}


inline bool check1(int x,int y){
//	(2*y-x,y,x]
    if(x==y) return 1;
    ll a=L.query(2*y-x+1,y)+dis[y];
    ll b=R.query(y,x-1)-dis[y];
    ll c=R.query(x,x)-dis[y];
    if(a<=c||b<=c) return 0;
    if(2*y-x>=1) return L.query(2*y-x,2*y-x)+dis[y]>=c;
    return 1;
}

inline bool check2(int x,int y){
//	[x,y,2*y-x)
    if(x==y) return 1;
    ll a=L.query(x+1,y)+dis[y];
    ll b=R.query(y,2*y-x-1)-dis[y];
    ll c=L.query(x,x)+dis[y];
    if(a<=c||b<=c) return 0;
    if(2*y-x<=n) return R.query(2*y-x,2*y-x)-dis[y]>c;
    return 1;
}

inline int solve1(int x){
    int l=1,r=x,mid,ans=0;
    while(l<=r){
        mid=(l+r)>>1;
        if(check1(x,mid)) r=mid-1,ans=mid;
        else l=mid+1;
    }
    return ans;
}

inline int solve2(int x){
    int l=x,r=n,mid,ans=0;
    while(l<=r){
        mid=(l+r)>>1;
        if(check2(x,mid)) l=mid+1,ans=mid;
        else r=mid-1;
    }
    return ans;
}

int main()
{
    n=read(),m=read();
    for(int i=2;i<=n;i++) lg[i]=lg[i>>1]+1;
    for(int i=2;i<=n;i++){
        w[i]=read();
        dis[i]=dis[i-1]+w[i];
    }
    ll ans;
    while(m--){
        k=read();
        for(int i=1;i<=k;i++) a[i].p=read(),a[i].l=read();
        sort(a+1,a+k+1);ans=0;
        for(int i=1;i<=k;i++) a[i].l-=dis[a[i].p];L.init();
        for(int i=1;i<=k;i++) a[i].l+=2*dis[a[i].p];R.init();
        for(int i=1;i<=k;i++) ans+=solve2(a[i].p)-solve1(a[i].p)+1;
        printf("%lld\n",ans);
    }
    return 0;
}