显然,每个点跳k步能到达的点一定是一个区间。
考虑倍增,l[i][k]表示点i跳2^k步能到的左端点,r[i][k]同理。
转移:
l[i][k]=min{l[j][k-1]} l[i][k-1]<=j<=r[i][k-1]
r[i][k]=max{r[j][k-1]} l[i][k-1]<=j<=r[i][k-1]
线段树搞搞。
二分答案求出l[ ],r[ ]后,就是问存不存在l[i]
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=100010;
int e[maxn][18][2],n,a[maxn],now[maxn][2],temp[maxn][2],s[maxn];
struct tree
{
tree *ls,*rs;
int mi,mx;
tree()
{
ls=rs=NULL;
}
void update()
{
mx=max(ls->mx,rs->mx);
mi=min(ls->mi,rs->mi);
}
void build(int lx,int rx,int k,int b)
{
if(lx==rx)
{
mi=mx=e[lx][k][b];
return ;
}
int mid=(lx+rx)>>1;
(ls=new tree)->build(lx,mid,k,b);
(rs=new tree)->build(mid+1,rx,k,b);
update();
}
int getmin(int l,int r,int lx,int rx)
{
if(l==lx&&r==rx) return mi;
int mid=(l+r)>>1;
if(rx<=mid) return ls->getmin(l,mid,lx,rx);
else if(lx>mid) return rs->getmin(mid+1,r,lx,rx);
else return min(ls->getmin(l,mid,lx,mid),rs->getmin(mid+1,r,mid+1,rx));
}
int getmax(int l,int r,int lx,int rx)
{
if(l==lx&&r==rx) return mx;
int mid=(l+r)>>1;
if(rx<=mid) return ls->getmax(l,mid,lx,rx);
else if(lx>mid) return rs->getmax(mid+1,r,lx,rx);
else return max(ls->getmax(l,mid,lx,mid),rs->getmax(mid+1,r,mid+1,rx));
}
}*xtr[18][2];
bool check(int x,int lb)
{
if ((x|(1<<lb))>n) return 0;
for(int i=1;i<=n;i++)
{
temp[i][0]=xtr[lb][0]->getmin(1,n,now[i][0],now[i][1]);
temp[i][1]=xtr[lb][1]->getmax(1,n,now[i][0],now[i][1]);
}
bool re=0;
s[0]=n+1;
for(int i=1;i<=n;i++)
{
s[i]=(i==1)?temp[i][1]:min(temp[i][1],s[i-1]);
if(s[temp[i][0]-1]<i){re=1; break;}
}
if(re==1)
for(int i=1;i<=n;i++)
{
now[i][0]=temp[i][0];
now[i][1]=temp[i][1];
}
return re;
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
e[i][0][0]=max(1,i-a[i]);
e[i][0][1]=min(i+a[i],n);
}
(xtr[0][0]=new tree)->build(1,n,0,0);
(xtr[0][1]=new tree)->build(1,n,0,1);
for(int k=1;(1<<k)<=n;k++)
{
for(int i=1;i<=n;i++)
{
e[i][k][0]=xtr[k-1][0]->getmin(1,n,e[i][k-1][0],e[i][k-1][1]);
e[i][k][1]=xtr[k-1][1]->getmax(1,n,e[i][k-1][0],e[i][k-1][1]);
}
(xtr[k][0]=new tree)->build(1,n,k,0);
(xtr[k][1]=new tree)->build(1,n,k,1);
}
int ans=0;
for(int i=1;i<=n;i++)
now[i][0]=now[i][1]=i;
for(int t=17;t>=0;t--)
if(check(ans,t)) ans=ans|(1<<t);
printf("%d\n",ans+1);
return 0;
}
/*8
7 1 1 1 1 1 1 7
10*/
卡内存。。。