带修莫队 的 小优化 (针对yxc版本)

带修莫队模板题

请大佬出门右转

y总进阶课有一道类似的题，我把进阶课的代码交上去发现t了好几个点，我以为我代码常数大，我把y总的代码交上去之后，一直被卡一个点没过。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

const int N = 10010, S = 1000010;

int n, m, mq, mc, len;
int w[N], cnt[S], ans[N];
struct Query
{
    int id, l, r, t;
}q[N];
struct Modify
{
    int p, c;
}c[N];

int get(int x)
{
    return x / len;
}

bool cmp(const Query& a, const Query& b)
{
    int al = get(a.l), ar = get(a.r);
    int bl = get(b.l), br = get(b.r);
    if (al != bl) return al < bl;
    if (ar != br) return ar < br;
    return a.t < b.t;
}

void add(int x, int& res)
{
    if (!cnt[x]) res ++ ;
    cnt[x] ++ ;
}

void del(int x, int& res)
{
    cnt[x] -- ;
    if (!cnt[x]) res -- ;
}

int main()
{
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= n; i ++ ) scanf("%d", &w[i]);
    for (int i = 0; i < m; i ++ )
    {
        char op[2];
        int a, b;
        scanf("%s%d%d", op, &a, &b);
        if (*op == 'Q') mq ++, q[mq] = {mq, a, b, mc};
        else c[ ++ mc] = {a, b};
    }

    len = cbrt((double)n * mc) + 1;
    sort(q + 1, q + mq + 1, cmp);

    for (int i = 0, j = 1, t = 0, k = 1, res = 0; k <= mq; k ++ )
    {
        int id = q[k].id, l = q[k].l, r = q[k].r, tm = q[k].t;
        while (i < r) add(w[ ++ i], res);
        while (i > r) del(w[i -- ], res);
        while (j < l) del(w[j ++ ], res);
        while (j > l) add(w[ -- j], res);
        while (t < tm)
        {
            t ++ ;
            if (c[t].p >= j && c[t].p <= i)
            {
                del(w[c[t].p], res);
                add(c[t].c, res);
            }
            swap(w[c[t].p], c[t].c);
        }
        while (t > tm)
        {
            if (c[t].p >= j && c[t].p <= i)
            {
                del(w[c[t].p], res);
                add(c[t].c, res);
            }
            swap(w[c[t].p], c[t].c);
            t -- ;
        }
        ans[id] = res;
    }

    for (int i = 1; i <= mq; i ++ ) printf("%d\n", ans[i]);
    return 0;
}

作者：yxc
链接：https://www.acwing.com/activity/content/code/content/500054/
来源：AcWing
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

跑一遍代码。

一直被卡，当我看到评论说当t很小的时候，几乎是和普通莫队一样的时候，采取块长为(n * mc)^1/3的大小去做的时候，我们观察到mc为0的时候，块长其实已经变成了1 。这个时候其实做法已经不能保证比较好的时间复杂度了，这时候需要我们进行一个优化，我们把len长度的赋值改为以下版本才能保证比较好的时间复杂度了，如果你采取块长为n^2/3 当我没说，本文只针对yxc版本的时间复杂度。

跑一遍代码。

y总的代码在t为0的时候有比较大的缺陷，基本上退化成暴力的时间复杂度，所以会超时，t为0的时候，我们可以考虑把他变成普通莫队，然后时间复杂度可以得以保证在O((n⁴t)^1/3)的时间复杂度。

蒟蒻说了一堆废话 自嗨

~~ 已订正~~

原代码已经会被卡掉了