标签：复杂度 logs lap bzoj 时间复杂度 long splay ref link
我其实是在codevs上看到它的题号后才去做这道题的。。。2333。。。
题目传送门：codevs:http://codevs.cn/problem/2333/ bzoj:http://www.lydsy.com/JudgeOnline/problem.php?id=2002
题意已经说的很清晰了，蓝鹅……我不会做。
看题解有两种写法：分块和link-cut tree。但是lct我不会，只能用分块了。
但是怎么分块呢？
我们想一下有什么暴力写法。一种是记录每个反弹装置能往后弹多长距离，这样效率是修改O(1)+查询O(n)。另一种是记录从每个反弹装置开始到被弹飞要被弹多少次，这样效率是修改O(n)+查询O(1)。我们可以想一个介于两者之间的做法，比如……把序列平均分成m块，记录每个反弹装置弹到下一个块要弹多少次 和 会弹到下一个块的哪一个装置。这样一来，效率就变成了修改O(n/m)+查询O(m)，总时间复杂度就是O(n*(n/m+m))。显然，如果要使效率最快m就取√n。
接下来就可以打代码+AC了。（代码真心很短）
代码：
var a,ne,x,y:array[0..200010]of longint; n,m,i,j,k,p,t,q,ans:longint; begin read(n); p:=trunc(sqrt(n)); for i:=0 to n-1 do begin read(a[i]); ne[i]:=(i div p+1)*p; end; for i:=0 to n-1 do begin x[i]:=i; y[i]:=0; while(x[i]<ne[i])and(x[i]<n)do begin x[i]:=x[i]+a[x[i]]; inc(y[i]); end; end; read(m); for i:=1 to m do begin read(q); if q=1 then begin read(k); ans:=0; while k<n do begin ans:=ans+y[k]; k:=x[k]; end; writeln(ans); end else begin read(k,t); a[k]:=t; for j:=ne[k]-1 downto ne[k]-p do if j+a[j]>=ne[j] then begin x[j]:=j+a[j]; y[j]:=1; end else begin x[j]:=x[j+a[j]]; y[j]:=y[j+a[j]]+1; end; end; end; end.
【codevs2333】&【BZOJ2002】弹飞绵羊[分块]
标签：复杂度 logs lap bzoj 时间复杂度 long splay ref link
原文地址：http://www.cnblogs.com/quzhizhou/p/6756655.html
