BZOJ2821: 作诗(Poetize)

2821: 作诗(Poetize)

Time Limit: 50 Sec  Memory Limit: 128 MB
Submit: 3300  Solved: 965
[Submit][Status][Discuss]

Description

神犇SJY虐完HEOI之后给傻×LYD出了一题：SHY是T国的公主，平时的一大爱好是作诗。由于时间紧迫，SHY作完诗

之后还要虐OI，于是SHY找来一篇长度为N的文章，阅读M次，每次只阅读其中连续的一段[l,r]，从这一段中选出一

些汉字构成诗。因为SHY喜欢对偶，所以SHY规定最后选出的每个汉字都必须在[l,r]里出现了正偶数次。而且SHY认

为选出的汉字的种类数（两个一样的汉字称为同一种）越多越好（为了拿到更多的素材！）。于是SHY请LYD安排选

法。LYD这种傻×当然不会了，于是向你请教……问题简述：N个数，M组询问，每次问[l,r]中有多少个数出现正偶

数次。

Input

输入第一行三个整数n、c以及m。表示文章字数、汉字的种类数、要选择M次。第二行有n个整数，每个数Ai在[1, c

]间，代表一个编码为Ai的汉字。接下来m行每行两个整数l和r，设上一个询问的答案为ans(第一个询问时ans=0)，

令L=(l+ans)mod n+1, R=(r+ans)mod n+1，若L>R，交换L和R，则本次询问为[L,R]。

Output

输出共m行，每行一个整数，第i个数表示SHY第i次能选出的汉字的最多种类数。

Sample Input

5 3 5
1 2 2 3 1
0 4
1 2
2 2
2 3
3 5

Sample Output

2
0
0
0
1

HINT

对于100%的数据，1<=n,c,m<=10^5

Source

【题解】

不知是什么原因，我用上次做蒲公英的写法写这道题，T了。。。

改用陈立杰论文的写法，又T了。。。

“写这个题以后，我也没干别的，大概三个写法，第一个，就是vector存位置二分它，第二个，就是分块排序二分，第三个，就是预处理每个块中每个数出现次数。如果说还有一点优化，就是加了输出优化，这个对时间秒杀hzw题解意义很重大的，还有取摸优化也是很重要的。很惭愧，就做了一点微小的工作，谢谢大家。” ——姜则岷

发现蒲公英也可以用这个写法，更快。。

设块大小为m预处理p[i][j]表示第i块到第j块有多少出现偶数次的，复杂度n^2/m

预处理pp[i][j]表示第i块，j出现了多少次，然后做前缀和，复杂度m * n（数的个数与n相同）

查询时，中间块直接获得，两边的数依次枚举，查询pp，根据奇偶分情况讨论，复杂度nm（询问次数与n相同）

总复杂度m * n + nm + n^2/m,设m = n^x

2 * n^(x + 1)  + n^(2 - x) <= 4√n^3  当且仅当x + 1 = 2 - x，即x = 0.5时取等号

m = √n

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <vector>
#include <map>
#include <string> 
#include <cmath> 
#include <sstream>
#define min(a, b) ((a) < (b) ? (a) : (b))
#define max(a, b) ((a) > (b) ? (a) : (b))
#define abs(a) ((a) < 0 ? (-1 * (a)) : (a))
template<class T>
inline void swap(T &a, T &b)
{
    T tmp = a;a = b;b = tmp;
}
inline void read(int &x)
{
    x = 0;char ch = getchar(), c = ch;
    while(ch < '0' || ch > '9') c = ch, ch = getchar();
    while(ch <= '9' && ch >= '0') x = x * 10 + ch - '0', ch = getchar();
    if(c == '-') x = -x;
}
const int INF = 0x3f3f3f3f;
const int MAXN = 100000 + 10;
const int SMAXN = 320;
int size,group,l[SMAXN],r[SMAXN],num[MAXN],n,c[MAXN],m,typ,p[SMAXN][SMAXN],sum[SMAXN][MAXN];
inline int find(int ll, int rr, int x)
{
    return sum[rr][x] - sum[ll - 1][x];
}
void make_p()
{
    register int i,j,k,ans,now;
    for(i = 1;i <= group;++ i)
    {
        ans = 0, now = i;
        for(int j = l[i];j <= r[i];++ j) ++ sum[i][num[j]];
        for(int j = l[i];j <= n;++ j)
        {
            ++ c[num[j]];
            if(c[num[j]] != 1) 
            {
                if(c[num[j]] & 1) -- ans;
                else ++ ans;
            }
            if(r[now] == j) p[i][now] = ans, ++ now;
        }
        for(int j = l[i];j <= n;++ j) c[num[j]] = 0;
    }
    for(i = 2;i <= group;++ i)
        for(int j = 1;j <= typ;++ j)
            sum[i][j] += sum[i - 1][j];
}
int query(int ll, int rr)
{
    int bl = (ll - 1) / size + 1, br = (rr - 1) / size + 1, ans = 0;
    if(bl + 2 <= br)
    {
        ans = p[bl + 1][br - 1];
        for(int i = ll;i <= r[bl];++ i) ++ c[num[i]];
        for(int i = l[br];i <= rr;++ i) ++ c[num[i]];
        for(int i = ll;i <= r[bl];++ i) 
            if(c[num[i]])
            {
                int tmp = 0;
                tmp = find(bl + 1, br - 1, num[i]);
                if((tmp & 1) && (c[num[i]] & 1)) ++ ans;
                else if(!tmp && !(c[num[i]] & 1)) ++ ans;
                else if(tmp && !(tmp & 1) && (c[num[i]] & 1)) -- ans; 
                c[num[i]] = 0;
            }
        for(int i = l[br];i <= rr;++ i) 
            if(c[num[i]])
            {
                int tmp = 0;
                tmp = find(bl + 1, br - 1, num[i]);
                if((tmp & 1) && (c[num[i]] & 1)) ++ ans;
                else if(!tmp && !(c[num[i]] & 1)) ++ ans;
                else if(tmp && !(tmp & 1) && (c[num[i]] & 1)) -- ans; 
                c[num[i]] = 0;
            }
    }
    else
    {
        for(int j = ll;j <= rr;++ j)
        {
            ++ c[num[j]];
            if(c[num[j]] == 1) continue;
            if(c[num[j]] & 1) -- ans;
            else ++ ans;
        }
        for(int j = ll;j <= rr;++ j) c[num[j]] = 0;
    }
    return ans;
}
std::ostringstream os;
int main()
{
    read(n), read(typ), read(m);
    for(register int i = 1;i <= n;++ i) read(num[i]);
    int M = 1;for(;(1 << M) <= n;++ M);-- M;
    size = sqrt(n);
    for(register int i = 1;i <= n;i += size)
        l[++ group] = i, r[group] = min(i + size - 1, n);
    make_p();
    register int pre = 0, l, r;
    for(register int i = 1;i <= m;++ i)
    {
        read(l), read(r);
        l = (l + pre), r = (r + pre);
        if(l >= n) l -= n;
        if(r >= n) r -= n;
        l += 1, r += 1;
        if(l > r) swap(l, r);
        pre = query(l, r);/*printf("%d", pre);putchar('
');*/
        os << pre << '
';
    }
    std::cout << os.str();
    return 0;
}