CodeForces 13E. Holes 分块处理

正解是动态树，太难了，仅仅好分块处理水之。看了看status大概慢了一倍之多。。

分块算法大体就是在找一个折衷点，使得查询和改动的时间复杂度都不算太高，均为o(sqrt(n)),所以总的时间复杂度为o(m*sqrt(n))。

分块的大体思想就是将整段区间分成sqrt(n)，每次改动影响所在段内的sqrt(n)个点。每次查询遍历sqrt(n)段，然后搞出一些事情。

对于此题，首先说明几个数组的含义：

cnt[i]的意义为从i点開始跳出其所在段所须要的步数。

goal[i]的意义为从i点開始所经历的最后一个点。

jmp[i]的意义为从i点開始跳出所在段之后第一个到达点。

明确了这几个意义之后看代码就好了，文字太无力。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <stack>
#include <map>
#include <ctime>
#include <iomanip>

#pragma comment(linker, "/STACK:1024000000");
#define EPS (1e-6)
#define LL long long
#define ULL unsigned long long
#define _LL __int64
#define INF 0x3f3f3f3f
#define Mod 1000000007

using namespace std;

const int MAXN = 100010,BLOCK = sqrt(100010);

int jmp[MAXN],cnt[MAXN],goal[MAXN],val[MAXN];

void Updata(int x,int y,int n)
{
if(y > n)
{
jmp[x] = MAXN;
cnt[x] = 1;
goal[x] = x;
}
else if(x/BLOCK*BLOCK == y/BLOCK*BLOCK)
{
jmp[x] = jmp[y];
cnt[x] = cnt[y]+1;
goal[x] = goal[y];
}
else
{
jmp[x] = y;
cnt[x] = 1;
goal[x] = y;
}
}

void Cal(int x)
{
int ans = 0,f;

while(x != MAXN)
{
f = goal[x];
ans += cnt[x];
x = jmp[x];
}

printf("%d %d\n",f,ans);

}

int main()
{
int i,j,u,v,n,m,ord,L;

scanf("%d %d",&n,&m);

for(i = 1;i <= n; ++i)
scanf("%d",&val[i]);

for(i = n;i >= 1; --i)
Updata(i,i+val[i],n);

while(m--)
{
scanf("%d",&ord);

if(ord)
{
scanf("%d",&u);
Cal(u);
}
else
{
scanf("%d %d",&u,&v);
val[u] = v;
L = u/BLOCK*BLOCK;
for(i = u;i >= L; --i)
Updata(i,i+val[i],n);
}
}

return 0;
}