LOJ #6282. 数列分块入门 6

单点插入，单点查询。

因为要么插入O(1)查询O(sqrt(n))，要么插入O(sqrt(n))，查询O(1)，所以选择代码难度更小的非链表暴力后移。

但是直接这样做虽然在随机数据上表现不错，但是在那些特意构造出的在某个区间疯狂加点的数据上会退化到与暴力差不多。

解决办法是每sqrt（n）次insert后将块重构，增大一点常数但使其保持均衡，从而保证复杂度不会被卡。

#include<iostream>
#include<cmath>
#include<vector>
#include<cstdio>
using namespace std;

void read(int& a){
    a = 0;char c = getchar(),l = c;
    while(c < '0'||c > '9')l = c,c = getchar();
    while('0' <= c&&c <= '9')a = a*10+c-'0',c = getchar();
    if(l == '-')a *= -1;
}

vector<int> ve[350];
int n,blo,cnti;

void insert(int x,int pl){
    int suf,p;
    for(suf = 0,p = 1;suf+ve[p].size() < pl;p++)
        suf += ve[p].size();
    pl -= suf+1;int y = *--ve[p].end();
    ve[p].push_back(y);
    for(register int i = ve[p].size()-2;i > pl;i--)ve[p][i] = ve[p][i-1];
    ve[p][pl] = x;
}

int ask(int r){
    int suf,p;
    for(suf = 0,p = 1;suf+ve[p].size() < r;p++)
        suf += ve[p].size();
    return ve[p][r-suf-1];
}

vector<int> all;
void reset(){
    all.clear();
    for(int i = 1;i <= (n-1)/blo+1;i++){
        for(int j = 0;j < ve[i].size();j++)
            all.push_back(ve[i][j]);
        ve[i].clear();
    }
    n = all.size(),blo = sqrt(n);
    for(int i = 0;i < n;i++)
        ve[i/blo+1].push_back(all[i]);
}

int main(){
    read(n); blo = sqrt(n);
    for(register int i = 1;i <= n;i++){
        int x; read(x);
        ve[(i-1)/blo+1].push_back(x);
    }
    for(register int i = n;i > 0;i--){
        int opt,l,r,c;
        read(opt),read(l),read(r),read(c);
        switch(opt){
            case 0:
                insert(r,l); cnti++;
                if(cnti%blo == 0)reset();
            break;
            case 1:
                printf("%d
",ask(r));
            break;
        }
    }
return 0;
}