BZOJ 2002: [Hnoi2010]Bounce 弹飞绵羊

2002: [Hnoi2010]Bounce 弹飞绵羊

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

——分割线——

弹飞绵羊QAQ这道题目名字很亲和，其实好坑爹！一开始评测还以为要T、结果A了。好开心！

好吧，言归正传，这道题目用分块的做法可以过。看Hzwer的博客貌似用动态树才是更好的做法Orz……

分块的做法就是将n个机器分成根号n块，处理个机器走几步可以走出本块，并且到达哪一块的那一个机器。然后球答案的时候最多只要算根号n次，而每次更新机器也最多只要根号n次，这样把两个的复杂度平衡了。就可以过了。

代码：【中途有一个变量写错了，特意写了个Arrest查了好久、QAQ】

#include<cstdio>
#include<cmath>
using namespace std;
int n;
struct Node{
	int num;
	int index[450];
	int k[450];
	int step[450];
	int to[450];
	int next[450];
};
Node node[450];
int m;
int sn;
int cnt=0;

inline void calc(){
	for (int i=1;i<=cnt;i++){
		for (int j=node[i].num;j>=1;j--){
			if (node[i].index[j]+node[i].k[j]>n){
				node[i].step[j]=1;
				node[i].next[j]=-1;
			}else{
				if (node[i].k[j]+j<=node[i].num){
					node[i].to[j]=node[i].to[j+node[i].k[j]];
					node[i].step[j]=node[i].step[j+node[i].k[j]]+1;
					node[i].next[j]=node[i].next[j+node[i].k[j]];
				}else{
					int tempcnt=i+1;
					int tempk=node[i].k[j]-(node[i].num-j);
					while(tempk>node[tempcnt].num){
						tempk-=node[tempcnt].num;
						tempcnt++;
					}
					node[i].to[j]=tempcnt;
					node[i].next[j]=tempk;
					node[i].step[j]=1;
				}
			}
		}
	}
}

inline void calc(int index,int val){
	node[(index-1)/sn+1].k[(index-1)%sn+1]=val;
	int i=(index-1)/sn+1;
	for (int j=(index-1)%sn+1;j>=1;j--){
		if (node[i].index[j]+node[i].k[j]>n){
				node[i].step[j]=1;
				node[i].next[j]=-1;
			}else{
				if (node[i].k[j]+j<=node[i].num){
					node[i].to[j]=node[i].to[j+node[i].k[j]];
					node[i].step[j]=node[i].step[j+node[i].k[j]]+1;
					node[i].next[j]=node[i].next[j+node[i].k[j]];
				}else{
					int tempcnt=i+1;
					int tempk=node[i].k[j]-(node[i].num-j);
					while(tempk>node[tempcnt].num){
						tempk-=node[tempcnt].num;
						tempcnt++;
					}
					node[i].to[j]=tempcnt;
					node[i].next[j]=tempk;
					node[i].step[j]=1;
				}
			}
	}
}


inline void Arrest(){
	printf("Call Cnts:\n");
	for (int i=1;i<=cnt;i++){
		printf("Cnt %d\n",i);
		for (int j=1;j<=node[i].num;j++){
			printf("num=%d index=%d k=%d step=%d toCnt=%d next=%d\n",j,node[i].index[j],node[i].k[j],node[i].step[j],node[i].to[j],node[i].next[j]);
		}
	}	
}

int main(){
	scanf("%d",&n);
	sn=sqrt(n);
	cnt=n/sn+(n%sn!=0?1:0);
	int index=0;
	for(int i=1;i<=n;i++){
		int temp=0;
		index++;
		scanf("%d",&temp);
		node[(index-1)/sn+1].k[(index-1)%sn+1]=temp;
		node[(index-1)/sn+1].index[(index-1)%sn+1]=index;
		node[(index-1)/sn+1].num++;
	}
	calc();
	scanf("%d",&m);
	for (int i=1;i<=m;i++){
		//Arrest();
		int order,j;
		scanf("%d%d",&order,&j);
		j++;
		if (order==1){
			int icnt=(j-1)/sn+1;
			int iindex=(j-1)%sn+1;
			int ans=0;
			while(iindex!=-1){
				//printf("icnt %d iindex %d -> Add %d\n",icnt,iindex,node[icnt].step[iindex]);
				ans=ans+node[icnt].step[iindex];
				int tindex=iindex;
				iindex=node[icnt].next[iindex];
				icnt=node[icnt].to[tindex];
			}
			printf("%d\n",ans);
		}else{
			int k;
			scanf("%d",&k);
			calc(j,k);
		} 
	}
	
	return 0;
}