BZOJ 2002: [Hnoi2010]Bounce 弹飞绵羊 分块

2002: [Hnoi2010]Bounce 弹飞绵羊

Time Limit:10 Sec  Memory Limit:259 MB
Submit: 10110  Solved: 5226
[Submit][Status][Discuss]

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

第一眼看到就认定是动态树

然而。。。我不会写

查了题解，   分块！！  思路超精巧

分块算法是一种很常见的根号算法，一般它的时间复杂度会带根号。 
分块和线段树的区别在于，分块算法可以维护一些线段树维护不了的东西，例如单调队列等，线段树能维护的东西必须能够进行信息合并，而分块则不需要。不过，它们也有共同点，分块和线段树一样，分块需要支持类似标记合并的东西。 
简单来说，分块算法就是优化过后的暴力。

#include<cmath>
#include<ctime>
#include<cstdio> 
#include<climits>
#include<cstring> 
#include<cstdlib> 
#include<iostream> 
#include<algorithm> 
#include<iomanip> 
#include<vector> 
#include<string>
#include<queue>  
#include<map> 
#include<set>
using namespace std; 
inline int read() 
{ 
	int x=0,f=1;char ch=getchar(); 
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();} 
	return f*x; 
}
const int N=200100; 
int n,m,k[N],belong[N],l[N],r[N],st[N],pt[N];
inline int cal(int x)
{
	int ans=0;
	while(1)
	{
		ans+=st[x];
		if(!pt[x])break;
		x=pt[x];
	}
	return ans;
}
int main()
{
	n=read();int block,num;
	for(int i=1;i<=n;i++)k[i]=read();
	block=sqrt(n);
	if(n%block)num=block+1;
	else num=block;
	for(int i=1;i<=n;i++)
		belong[i]=(i-1)/block+1;
	for(int i=1;i<=num;i++)
		l[i]=(i-1)*block+1,r[i]=block*i;
	r[num]=n;
	for(int i=n;i>0;i--)
	{
		if(i+k[i]>n)st[i]=1;
		else if(belong[i]==belong[i+k[i]])
		{st[i]=st[i+k[i]]+1;pt[i]=pt[i+k[i]];}
		else st[i]=1,pt[i]=i+k[i];
	}
	m=read();int opt,x,y;
	for(int i=1;i<=m;i++)
	{
		opt=read();x=read()+1;
		if(opt==1)printf("%d\n",cal(x));
		else 
		{
			y=read();k[x]=y;
			for(int i=x;i>=l[belong[x]];i--)
			{
				if(belong[i]==belong[i+k[i]])
				{st[i]=st[i+k[i]]+1;pt[i]=pt[i+k[i]];}
				else st[i]=1,pt[i]=i+k[i];
			}
		}
	}
}