「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;
}