【BZOJ 2724】蒲公英 &【BZOJ 2821】作诗【分块】

今天全是奇妙的分块大法。

【BZOJ 2714】蒲公英

给一个序列，求一个区间的众数。（强制在线）
【题解】
先分块。
　　设f[i][j]表示第i块到第j块总的众数。
　　然后每次对于一个询问[l,r]，x到y已经知晓，答案只可能是l~x，y~r，和块x~y的众数之一，所以接下来要做的就是判断l~x和y~r中存不存在数字出现次数比块x到y众数出现更多的数。
　　不过怎么判断呢？可以先用vector来保存每个数字出现的位置（因为本题的值域较大，所以要离散化一下，就是map来搞一搞），在用二分答案查找这个数字在[l,r]内出现的次数即可。
详见代码。
【代码】

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>

#define N 40010
#define B 200
using namespace std;

map<int,int> mp;
vector<int> pos[N];

int cnt,val[N],a[N];
int belong[N],tot[N],f[B+5][B+5];
int n;

void pre(int x)
{
    int now = 0,maxnsum = 0;
    memset(tot,0,sizeof(tot));
    for(int i = (x-1)*B + 1;i <= n;i++) {
        tot[a[i]] ++;
        if(tot[a[i]] > maxnsum || (tot[a[i]] == maxnsum && val[now] > val[a[i]]) )
            now = a[i],maxnsum = tot[a[i]];
        f[x][belong[i]] = now;
    }
}

int sigma(int l,int r,int x) {
    return upper_bound(pos[x].begin(),pos[x].end(),r)-lower_bound(pos[x].begin(),pos[x].end(),l);
}

int query(int l,int r)
{
    int now = f[belong[l]+1][belong[r]-1],maxnsum = sigma(l,r,now);
    for(int i = l;i <= min(r,belong[l]*B);i++) {
        int tmp = sigma(l,r,a[i]);
        if(tmp > maxnsum || (tmp == maxnsum && val[now] > val[a[i]]) )
            now = a[i],maxnsum = tmp;
    }
    if(belong[l] != belong[r]) {
        for(int i = (belong[r]-1)*B+1;i <= r;i++) {
            int tmp = sigma(l,r,a[i]);
            if(tmp > maxnsum || (tmp == maxnsum && val[now] > val[a[i]]) )
                now = a[i],maxnsum = tmp;
        }
    }
    return now;
}

int main()
{
    int m,l,r;
    int ans = cnt = 0;
    scanf("%d%d",&n,&m);
    for(int i = 1;i <= n;i++) {
        scanf("%d",&a[i]);
        if(!mp[a[i]]) {
            mp[a[i]] = ++cnt; val[cnt] = a[i];
        }
        a[i] = mp[a[i]];
        pos[a[i]].push_back(i);
    }
    for(int i = 1;i <= n;i++) belong[i] = (i-1)/B + 1;
    for(int i = 1;i <= belong[n];i++) pre(i);

    for(int i = 1;i <= m;i++) {
        scanf("%d%d",&l,&r);
        l = (l+ans-1) % n + 1; r = (r+ans-1) % n + 1;
        if(l > r) swap(l,r);
        ans = val[query(l,r)];
        printf("%d\n",ans);
    }
    return 0;
}

【BZOJ 2821】作诗

给一个序列，求一个区间内出现正偶次的数的个数。（强制在线）
【题解】
本题和上面一题的套路是一样的。
　　用f[i][j]表示第i块到第j块中符合条件的数的个数。
　　对于一个询问l~r，中间大块x~y已经知晓，要做的就是看l-x和y-r之间的数对答案的影响。
　　用vector保存每个数的位置，分别统计除了大块外的那些数（不超过2*sqrt(n)个）在x~y中出现的次数t1，以及在l~r之间出现的次数t2，根据t1和t2的奇偶性判断对答案的影响。同样，统计t1、t2也需要在vector中二分。
注：对f[i][j]的预处理可以在O(n * sqrt(n))的时间完成。
　　同时，块的大小取sqrt(n/log n)时最优。
【代码】

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

#define N 100010
#define B 1505
using namespace std;
vector<int> pos[N];
int a[N],belong[N],L[B],R[B];
int n,bs;
int f[B+5][B+5],tmp[N];
bool vis[N];

void pre() {
    int tot;
    for(int i = 1;i <= bs;i++) {
        for(int j = L[i];j <= n;j++) tmp[a[j]] = 0;
        tot = 0;
        for(int j = L[i];j <= n;j++) {
            if(!(tmp[a[j]]&1) && tmp[a[j]] != 0) tot--;
            tmp[a[j]]++;
            if(!(tmp[a[j]]&1) )tot++;
            f[i][belong[j]] = tot;
        }
    }
}

int sigma(int l,int r,int x) {
    return max(upper_bound(pos[x].begin(),pos[x].end(),r)-lower_bound(pos[x].begin(),pos[x].end(),l),0);
}

int query(int l,int r)
{
    int x = belong[l],y = belong[r],ans = 0,t1,t2;
    memset(vis,0,sizeof(vis));
    if(x == y || x+1 == y) {
        for(int i = l;i <= r;i++)
            if(!vis[a[i]]) {
                if(!(sigma(l,r,a[i])&1)) ans++;
                vis[a[i]] = 1;
            }
    }
    else {
        x++; y--;
        ans = f[x][y];
        for(int i = l;i <= L[x]-1;i++)
            if(!vis[a[i]]) {
                t2 = sigma(l,L[x]-1,a[i])+sigma(R[y]+1,r,a[i]);
                t1 = sigma(L[x],R[y],a[i]);
                if(!(t2&1) && t2 != 0 && t1 == 0) ans++;
                if((t2&1) && (t1&1)) ans++;
                if((t2&1) && !(t1&1) && t1 != 0) ans--;
                vis[a[i]] = 1;
            }
        for(int i = R[y]+1;i <= r;i++)
            if(!vis[a[i]]) {
                t2 = sigma(l,L[x]-1,a[i])+sigma(R[y]+1,r,a[i]);
                t1 = sigma(L[x],R[y],a[i]);
                if(!(t2&1) && t2 != 0 && t1 == 0) ans++;
                if((t2&1) && (t1&1)) ans++;
                if((t2&1) && !(t1&1) && t1 != 0) ans--;
                vis[a[i]] = 1;
            }
    }
    return ans;
}

int main()
{
    int c,m,l,r;
    scanf("%d%d%d",&n,&c,&m);
    for(int i = 1;i <= n;i++) {
        scanf("%d",&a[i]);
        pos[a[i]].push_back(i);
    }
    int block=sqrt((double)n/log((double)n)*log(2));
    if(n % block) bs = n/block+1; else bs = n/block;
    for(int i = 1;i <= n;i++) belong[i] = (i-1)/block+1;
    for(int i = 1;i <= bs;i++) {
        L[i] = (i-1)*block+1; R[i] = i*block;
    }
    R[bs] = n;
    pre();
    int ans = 0;
    while(m--) {
        scanf("%d%d",&l,&r);
        l = (l+ans) % n + 1;r = (r+ans) % n + 1;
        if(l > r) swap(l,r);
        ans = query(l,r);
        printf("%d\n",ans);
    }
    return 0;
}