bzoj2002 [Hnoi2010]Bounce 弹飞绵羊

http://www.elijahqi.win/2018/02/12/bzoj2002/ ‎

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

方法：分块/lct

分块：son为该节点指向的节点step为跳出这个块需要的步数 然后根号n分组暴力模拟即可

#include<cstdio>
#include<cmath>
#define N 220000
struct node{
    int son,step;
}map[N];
int n,m,a[N],q;
int query1(int posi){
    int ans=0;for (int i=posi;i!=0;i=map[i].son) ans+=map[i].step;
    return ans;
}
void rebuild(int p,int x){
    a[p]=x;int op,cl;
    if (p/q*q!=p) op=p/q*q+1,cl=(p/q+1)*q;else op=(p/q-1)*q+1,cl=p;
    for (int i=p;i>=op;--i)
        if (a[i]+i>cl)map[i].son=a[i]+i,map[i].step=1;else
        map[i].son=map[a[i]+i].son,map[i].step=map[a[i]+i].step+1;
}
void build(){
    for (int i=n;i>=1;--i) rebuild(i,a[i]);
}
int main(){
    freopen("bounce.in","r",stdin);
    freopen("bounce.out","w",stdout);
    scanf("%d",&n);
    for (int i=1;i<=n;++i) scanf("%d",&a[i]);
    for (int i=1;i<N;++i) map[i].son=0;
    q=sqrt(n);build();
    scanf("%d",&m);int tmp,tmp1,k;
    for (int i=1;i<=m;++i){
        scanf("%d",&tmp);
        if (tmp==1) scanf("%d",&tmp1),printf("%d\n",query1(tmp1+1));
        if (tmp==2) scanf("%d%d",&tmp1,&k),rebuild(tmp1+1,k);
    }
    return 0;
}

lct

增加一个n+1这个点 每次跳跃的时候算一算 如果超过n了 那么就直接把我这条链接在n+1上

#include<cstdio>
#include<algorithm>
#define N 220000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S)return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int size[N],c[N][2],fa[N],rev[N],s[N],q[N],top,n,m;
inline void update(int x){
    size[x]=size[c[x][0]]+size[c[x][1]]+1;
}
inline bool isroot(int x){
    return c[fa[x]][0]!=x&&c[fa[x]][1]!=x;
}
inline void pushdown(int x){
    if (!rev[x]) return;
    int l=c[x][0],r=c[x][1];
    swap(c[x][0],c[x][1]);
    rev[l]^=1;rev[r]^=1;rev[x]=0;
}
inline void rotate(int x){
    int y=fa[x],z=fa[y];
    if (!isroot(y)) c[z][c[z][1]==y]=x;
    int l=c[y][1]==x,r=l^1;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void print(int x){
    pushdown(x);
    if (c[x][0]) print(c[x][0]);
    printf("%d %d %d %d\n",x,fa[x],c[x][0],c[x][1]);
    if (c[x][1]) print(c[x][1]);
}
inline void splay(int x){
    q[top=1]=x;for (int i=x;!isroot(i);i=fa[i]) q[++top]=fa[i];
    for (int i=top;i;--i) pushdown(q[i]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if (!isroot(y)) if (c[y][0]==x^c[z][0]==y) rotate(x);else rotate(y);//,print(y),puts("fd");
        rotate(x);//print(x);
    }
}
inline void access(int x){
    for (int t=0;x;t=x,x=fa[x]) splay(x),c[x][1]=t,update(x);
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void cut(int x,int y){
    makeroot(y);access(x);splay(x);c[x][0]=fa[c[x][0]]=0;update(x);
}

int main(){
    //freopen("bzoj2002.in","r",stdin);
    n=read();for (int i=1;i<=n;++i) {int x=read();s[i]=x;if (i+x>n) x=n+1;else x=x+i;fa[i]=x;}
    m=read();for (int i=1;i<=n;++i) size[i]=1;
    //for (int i=1;i<=n+1;++i) printf("%d ",fa[i]);
    while(m--){
        int op=read();
        if (op==1){
            int st=read()+1;makeroot(n+1);
            access(st);//print(n+1);puts("sdf");
            splay(st);
            printf("%d\n",size[st]-1);
        }
        if (op==2){
            int x=read()+1,y=read();int to=x+s[x]>n?n+1:x+s[x];cut(x,to);s[x]=y;
            to=x+s[x]>n?n+1:x+s[x];makeroot(x);fa[x]=to;
        }
    }
    return 0;
}