【BZOJ】3922: Karin的弹幕

题意

给定一个长度为\(n(1 \le n \le 70000)\)序列，\(m(1 \le m \le 70000)\)次操作：1. 对一段下标是等差数列的子序列求最大值；2. 单点修改。

分析

如果公差很大，那么速度是很快的。所以我们考虑阈值。

题解

设阈值\(k\)表示如果询问的公差小于等于\(k\)则用线段树求，否则暴力。
则我们对阈值每一个公差\(i(1 \le i \le k)\)建\(i\)棵线段树。
复杂度\(O(kn + m(\frac{n}{k} + log n))\)

#include <bits/stdc++.h>
using namespace std;
const int N=70005, MXD=10, oo=(~0u>>1)+1;
int a[N], D, sz[MXD+1][MXD+1], b[N];
struct node *null;
struct node {
    node *c[2];
    int mx;
    node() {
        c[0]=c[1]=0;
        mx=oo;
    }
    void up() {
        mx=max(c[0]->mx, c[1]->mx);
    }
}Po[10000005], *iT=Po, *root[MXD+1][MXD+1];
node *newnode() {
    return iT++;
}
node* build(int l, int r) {
    node *x=newnode();
    if(l==r) {
        x->mx=b[l];
        return x;
    }
    int mid=(l+r)>>1;
    x->c[0]=build(l, mid);
    x->c[1]=build(mid+1, r);
    x->up();
    return x;
}
int query(int L, int R, int l, int r, node *x) {
    if(L<=l && r<=R) {
        return x->mx;
    }
    int mid=(l+r)>>1, ret=oo;
    if(L<=mid) {
        ret=query(L, R, l, mid, x->c[0]);
    }
    if(mid<R) {
        ret=max(ret, query(L, R, mid+1, r, x->c[1]));
    }
    return ret;
}
void update(int p, int go, int l, int r, node *x) {
    if(l==r) {
        x->mx+=go;
        return;
    }
    int mid=(l+r)>>1;
    if(p<=mid) {
        update(p, go, l, mid, x->c[0]);
    }
    else {
        update(p, go, mid+1, r, x->c[1]);
    }
    x->up();
}
void init(int n) {
    D=min(MXD, n);
    for(int d=1; d<=D; ++d) {
        for(int i=1; i<=d; ++i) {
            int &s=sz[d][i];
            for(int j=i; j<=n; j+=d) {
                b[++s]=a[j];
            }
            root[d][i]=build(1, s);
        }
    }
}
int query(int x, int d, int n) {
    if(d>D) {
        int mx=oo;
        for(int i=x; i<=n; i+=d) {
            mx=max(mx, a[i]);
        }
        return mx;
    }
    int begin=(x-1)%d+1, pos=(x-1)/d+1, len=sz[d][begin];
    return query(pos, len, 1, len, root[d][begin]);
}
void update(int x, int y) {
    a[x]+=y;
    for(int d=1; d<=D; ++d) {
        int begin=(x-1)%d+1, pos=(x-1)/d+1, len=sz[d][begin];
        update(pos, y, 1, len, root[d][begin]);
    }
}
int main() {
    int n, m;
    scanf("%d", &n);
    for(int i=1; i<=n; ++i) {
        scanf("%d", &a[i]);
    }
    init(n);
    scanf("%d", &m);
    while(m--) {
        int op, x, y;
        scanf("%d%d%d", &op, &x, &y);
        if(op) {
            printf("%d\n", query(x, y, n));
        }
        else {
            update(x, y);
        }
    }
    return 0;
}

来源：https://www.cnblogs.com/iwtwiioi/p/4986313.html