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;

}