noip模拟赛 街灯

分析：对于前30%的数据直接暴力模拟即可,对于另外30%的数据，因为每次的p是一样的，所以可以用莫队来维护，先离散化一下，再用一个桶统计次数.

100%的做法和之前做过的一道模拟赛题很像,当模数很小的时候分块，否则就暴力修改.其实看到区间操作第一感觉是线段树，但是线段树并不能维护这个,分块维护的信息多一些，所以分块.在模数较小的时候记录一下第i个块，模p等于v的有多少个，即g[i][p][v],利用前缀和统计1~i个块的个数.在模数较大的时候因为只有v,v+p,v+2p对答案有影响，所以记录第i个块值为v的有多少个,即f[i][v],同样也可以用前缀和处理一下.查询的时候还是分模数的大小来进行，先统计完整包含在块里面的，在暴力统计在块外面的就可以了.

当有些信息线段树维护不了的时候可以考虑一下分块，分块的时候可以根据范围来决定什么时候分块，什么时候用其它的方法,前缀和可以加速区间查询操作.

60分暴力：

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

using namespace std;

int n, m, a[], b[], cnt, tot, c[], sizee, L = , R = , tong[], ans[];
bool flag = true;

struct node
{
    int l, r, p, v,id;
}e[];

bool cmp(node a, node b)
{
    if (a.l / sizee == b.l / sizee)
        return a.r < b.r;
    return a.l < b.l;
}

void add(int x)
{
    tong[a[x]]++;
}

void del(int x)
{
    tong[a[x]]--;
}

int main()
{
    scanf("%d%d", &n, &m);
    for (int i = ; i <= n; i++)
        scanf("%d", &a[i]);
    cnt = ;
    for (int i = ; i <= m; i++)
    {
        scanf("%d%d%d%d", &e[i].l, &e[i].r, &e[i].p, &e[i].v);
        e[i].id = i;
        if (i !=  && e[i].p != e[i - ].p)
            flag = false;
        b[++cnt] = e[i].v;
    }
    if (flag)
    {
        for (int i = ; i <= n; i++)
        {
            a[i] %= e[].p;
            b[++cnt] = a[i];
        }
        memcpy(c, b, sizeof(b));
        sort(c + , c +  + cnt);
        tot = unique(c + , c +  + cnt) - c - ;
        for (int i = ; i <= n; i++)
            a[i] = lower_bound(c + , c +  + tot, a[i]) - c - ;
        for (int i = ; i <= m; i++)
            e[i].v = lower_bound(c + , c +  + tot, e[i].v) - c - ;
        sizee = (int)sqrt(n);
        sort(e + , e +  + m, cmp);
        for (int i = ; i <= m; i++)
        {
            int l = e[i].l, r = e[i].r;
            while (R < r)
                add(++R);
            while (R > r)
                del(R--);
            while (L > l)
                add(--L);
            while (L < l)
                del(L++);
            ans[e[i].id] = tong[e[i].v];
        }
        for (int i = ; i <= m; i++)
            printf("%d\n", ans[i]);
    }
    else
    {
        for (int i = ; i <= m; i++)
        {
            int cnt = ;
            for (int j = e[i].l; j <= e[i].r; j++)
                if (a[j] % e[i].p == e[i].v)
                    cnt++;
            printf("%d\n", cnt);
        }
    }

    return ;
}

100分正解：

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

using namespace std;

int n, m, a[], maxn, block, cnt, maxx, l[], r[], ans;
int f[][], g[][][], kuai[];

int main()
{
    scanf("%d%d", &n, &m);
    for (int i = ; i <= n; i++)
    {
        scanf("%d", &a[i]);
        maxn = max(maxn, a[i]);
    }
    block = ;
    cnt = (n - ) / block + ;
    maxx = ;
    for (int i = ; i <= cnt; i++)
    {
        l[i] = (i - ) * block + , r[i] = min(n, i * block);
        for (int j = ; j <= maxn; j++)
            f[i][j] = f[i - ][j];
        for (int j = ; j <= maxx; j++)
            for (int k = ; k < maxx; k++)
                g[i][j][k] = g[i - ][j][k];
        for (int j = l[i]; j <= r[i]; j++)
        {
            kuai[j] = i;
            f[i][a[j]]++;
            for (int k = ; k <= maxx; k++)
                g[i][k][a[j] % k]++;
        }
    }
    for (int i = ; i <= m; i++)
    {
        int ll, rr, p, v;
        scanf("%d%d%d%d", &ll, &rr, &p, &v);
        ans = ;
        v %= p;
        int L = kuai[ll], R = kuai[rr];
        if (L < R)
        {
            if (p <= )
                ans += g[R - ][p][v] - g[L][p][v];
            else
                for (int i = v; i <= maxn; i += p)
                    ans += f[R - ][i] - f[L][i];
            for (int j = ll; j <= r[L]; j++)
                if (a[j] % p == v)
                    ans++;
            for (int j = l[R]; j <= rr; j++)
                if (a[j] % p == v)
                    ans++;
        }
        else
            for (int j = ll; j <= rr; j++)
                if (a[j] % p == v)
                    ans++;
        printf("%d\n", ans);
    }

    return ;
}