弹飞绵羊（bzoj 2002）

Description

某天，Lostmonkey发明了一种超级弹力装置，为了在他的绵羊朋友面前显摆，他邀请小绵羊一起玩个游戏。游戏一开始，Lostmonkey在地上沿着一条直线摆上n个装置，每个装置设定初始弹力系数ki，当绵羊达到第i个装置时，它会往后弹ki步，达到第i+ki个装置，若不存在第i+ki个装置，则绵羊被弹飞。绵羊想知道当它从第i个装置起步时，被弹几次后会被弹飞。为了使得游戏更有趣，Lostmonkey可以修改某个弹力装置的弹力系数，任何时候弹力系数均为正整数。

Input

第一行包含一个整数n，表示地上有n个装置，装置的编号从0到n-1,接下来一行有n个正整数，依次为那n个装置的初始弹力系数。第三行有一个正整数m，接下来m行每行至少有两个数i、j，若i=1，你要输出从j出发被弹几次后被弹飞，若i=2则还会再输入一个正整数k，表示第j个弹力装置的系数被修改成k。对于20%的数据n,m<=10000，对于100%的数据n<=200000,m<=100000

Output

对于每个i=1的情况，你都要输出一个需要的步数，占一行。

Sample Input

4
1 2 1 1
3
1 1
2 1 1
1 1

Sample Output

2
3

/*
  求每棵子树的节点数。
  这个题很巧妙的是建立了一个n+1节点，但是不太清楚这个n+1节点在程序中起了什么作用。。。
*/
#include<iostream>
#include<cstdio>
#define N 200010
using namespace std;
int n,m;
int next[N],c[N][],fa[N],size[N],st[N];
bool rev[N];
bool isroot(int k){
    return c[fa[k]][]!=k&&c[fa[k]][]!=k;
}
void pushup(int x){
    size[x]=size[c[x][]]+size[c[x][]]+;
}
void pushdown(int k){
    int l=c[k][],r=c[k][];
    if(rev[k]){
        rev[k]^=;rev[l]^=;rev[r]^=;
        swap(c[k][],c[k][]);
    }
}
void rotate(int x){
    int y=fa[x],z=fa[y],l,r;
    if(c[y][]==x)l=;else l=;r=l^;
    if(!isroot(y)){
        if(c[z][]==y)c[z][]=x;else c[z][]=x;
    }
    fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
    c[y][l]=c[x][r];c[x][r]=y;
    pushup(y);pushup(x);
}
void splay(int x){
    int top=;st[++top]=x;
    for(int i=x;!isroot(i);i=fa[i])
        st[++top]=fa[i];
    for(int i=top;i;i--)pushdown(st[i]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(c[y][]==x^c[z][]==y)rotate(x);
            else rotate(y);
        }
        rotate(x);
    }
}
void access(int x){
    int t=;
    while(x){
        splay(x);
        c[x][]=t;
        t=x;x=fa[x];}
}
void rever(int x){
    access(x);splay(x);rev[x]^=;
}
void link(int x,int y){
    rever(x);fa[x]=y;splay(x);
}
void cut(int x,int y){
    rever(x);access(y);splay(y);c[y][]=fa[x]=;
}
int main(){
    scanf("%d",&n);
    for(int i=;i<=n;i++){
        int x;scanf("%d",&x);
        fa[i]=x+i;size[i]=;
        if(fa[i]>n+)fa[i]=n+;
        next[i]=fa[i];
    }
    size[n+]=;
    scanf("%d",&m);
    for(int i=;i<=m;i++){
        int f;scanf("%d",&f);
        if(f==){
            rever(n+);
            int x;scanf("%d",&x);x++;
            access(x);splay(x);printf("%d\n",size[c[x][]]);
        }
        else {
            int x,y;scanf("%d%d",&x,&y);x++;
            int t=min(n+,x+y);
            cut(x,next[x]);link(x,t);next[x]=t;
        }
    }
    return ;
}

/*
  如果没有修改操作，可能会想预处理出在每个点被弹飞的步数，这样的查询是O(1)的，但由于有修改操作，
  使每次操作的复杂度可能变成O(n)，所以可以考虑一种折中的方法，分块！
  对于每一块，处理出每个点跳出这个块的步数以及会落在哪一个点上，这样使查询和修改都变成了O(√n)。
*/
#include<cstdio>
#include<iostream>
#include<cmath>
#define N 200010
using namespace std;
int k[N],bl[N],stp[N],next[N],n,m,len;
void init(int x){
    for(int i=(x-)*len+;i<=min(x*len,n);i++){
        int tot=,pos=i;
        while(pos<=n&&bl[pos]==x){
            tot++;
            pos+=k[pos];
        }
        stp[i]=tot;next[i]=pos;
    }
}
int query(int x){
    int ans=;
    for(int i=x;i<=n;i=next[i])
        ans+=stp[i];
    return ans;
}
void modify(int x,int y){
    k[x]=y;
    for(int i=x;i>=(bl[x]-)*len+;i--){
        if(i+k[i]>bl[x]*len) stp[i]=,next[i]=i+k[i];
        else stp[i]=stp[i+k[i]]+,next[i]=next[i+k[i]];
    }
}
int main(){
    scanf("%d",&n);len=sqrt(n);
    for(int i=;i<=n;i++){
        scanf("%d",&k[i]);
        bl[i]=(i-)/len+;
    }
    for(int i=;i<=bl[n];i++)
        init(i);
    scanf("%d",&m);
    for(int i=;i<=m;i++){
        int opt,x,y;scanf("%d%d",&opt,&x);x++;
        if(opt==) printf("%d\n",query(x));
        else scanf("%d",&y),modify(x,y);
    }
    return ;
}