luoguP4168: [Violet]蒲公英

题目背景

亲爱的哥哥:

你在那个城市里面过得好吗?

我在家里面最近很开心呢。昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多小朋友也被它杀掉了。我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢。不过奶奶说他是很难受的时候才做出这样的事的……

最近村子里长出了一大片一大片的蒲公英。一刮风,这些蒲公英就能飘到好远的地方了呢。我觉得要是它们能飘到那个城市里面,让哥哥看看就好了呢!

哥哥你要快点回来哦!

爱你的妹妹 Violet

Azure 读完这封信之后微笑了一下。

“蒲公英吗……”

题目描述

在乡下的小路旁种着许多蒲公英,而我们的问题正是与这些蒲公英有关。

为了简化起见,我们把所有的蒲公英看成一个长度为n的序列 (a_1,a_2..a_n)(a1,a2..an) ,其中 a_iai 为一个正整数,表示第i棵蒲公英的种类编号。

而每次询问一个区间 [l,r],你需要回答区间里出现次数最多的是哪种蒲公英,如果有若干种蒲公英出现次数相同,则输出种类编号最小的那个。

注意,你的算法必须是在线的

输入输出格式

输入格式:

第一行两个整数 n,m ,表示有n株蒲公英,m 次询问。

接下来一行n个空格分隔的整数 a_iai ,表示蒲公英的种类

再接下来m 行每行两个整数 l_0,r_0l0,r0 ,我们令上次询问的结果为 x(如果这是第一次询问, 则 x=0)。

令 l=(l_0+x-1)\bmod n + 1,r=(r_0+x-1) \bmod n + 1l=(l0+x1)modn+1,r=(r0+x1)modn+1 ,如果 l>r,则交换 l,r 。

最终的询问区间为[l,r]。

输出格式:

输出m 行。每行一个整数,表示每次询问的结果。

输入输出样例

输入样例#1:  复制
6 3 
1 2 3 2 1 2 
1 5 
3 6 
1 5
输出样例#1:  复制
1 
2 
1

说明

对于 20% 的数据,保证 1\le n,m \le 30001n,m3000 。

对于 100% 的数据,保证 1\le n \le 40000,1\le m \le 50000,1\le a_i \le 10^91n40000,1m50000,1ai109 。

题解:

考虑分块,每次询问的答案为询问区间内所有完整的块的众数或者散块中出现过的数

对于每一个数,先离散化,开个vector存下每个数在序列中出现的位置,然后预处理出f[i][j],表示第i个块到第j个块的众数,对于每次询问,考察散块中的每个数及相应的f[i][j],即在vector中二分查找对应的数在询问区间内最近和最远的出现位置,其在vetcor中的位置相减加1即为其在询问区间中的出现次数

另外这题不吸氧气过不去......

#include<bits/stdc++.h>
using namespace std;
int n,m,sz,tot;
int arr[40005],tp[40005],f[205][205],idx[40005],cnt[40005],val[40005];
vector<int>pos[40005];
inline int read()
{
    int X=0,w=0; char ch=0;
    while(!isdigit(ch)) {w|=ch=='-';ch=getchar();}
    while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
inline int disc(int a)
{
    return lower_bound(tp+1,tp+1+n,a)-tp;
}
inline void pre_work(int k)
{
    int mx=0,num;
    memset(cnt,0,sizeof(cnt));
    for(register int i=(k-1)*sz+1;i<=n;i++)
        {
            cnt[arr[i]]++;
            if(mx<=cnt[arr[i]])
               {
                      if(mx==cnt[arr[i]])
                         {
                            if(num>arr[i])
                               num=arr[i];
                         }
                      else
                         {
                               mx=cnt[arr[i]];
                               num=arr[i];
                         }   
               }
            f[k][idx[i]]=num;   
        }
}
inline int find(int l,int r,int p)
{
    return upper_bound(pos[p].begin(),pos[p].end(),r)-lower_bound(pos[p].begin(),pos[p].end(),l);
}
inline int query(int l,int r)
{
    int ret,mx=0;
    if(idx[l]==idx[r])
       {
          for(register int i=l;i<=r;i++)
              {
                    int t=find(l,r,arr[i]);
                    if(mx<=t)
                       {
                            if(mx==t)
                               {
                                   if(val[arr[i]]<ret)
                                      ret=val[arr[i]];
                               }
                            else
                            {
                                mx=t;
                                ret=val[arr[i]];
                            }   
                       }
              }
       }
    else
       {
             for(register int i=l;i<=idx[l]*sz;i++)
                 { 
                     int t=find(l,r,arr[i]);
                    if(mx<=t)
                       {
                            if(mx==t)
                               {
                                   if(val[arr[i]]<ret)
                                      ret=val[arr[i]];
                               }
                            else
                            {
                                mx=t;
                                ret=val[arr[i]];
                            }   
                       }
                 }
             for(register int i=(idx[r]-1)*sz+1;i<=r;i++)
                 {
                       int t=find(l,r,arr[i]);
                    if(mx<=t)
                       {
                            if(mx==t)
                               {
                                   if(val[arr[i]]<ret)
                                      ret=val[arr[i]];
                               }
                            else
                            {
                                mx=t;
                                ret=val[arr[i]];
                            }   
                       }
                 }
              if(idx[l]+1<=idx[r]-1)
              {
                    int t=find(l,r,f[idx[l]+1][idx[r]-1]);
                    if(mx<=t)
                       {
                            if(mx==t)
                               {
                                   if(val[f[idx[l]+1][idx[r]-1]]<ret)
                                      ret=val[f[idx[l]+1][idx[r]-1]];
                               }
                            else
                            {
                                mx=t;
                                ret=val[f[idx[l]+1][idx[r]-1]];
                            }   
                       }
              }   
       }
    return ret;      
}
int main()
{
    n=read();m=read();
    sz=(int)sqrt(n);
    for(register int i=1;i<=n;i++)
        {
            arr[i]=read();
            tp[i]=arr[i];
            idx[i]=(i-1)/sz+1;
        }
    sort(tp+1,tp+1+n);
    for(register int i=1;i<=n;i++)
        {
            int r=disc(arr[i]);
            val[r]=arr[i];
            arr[i]=r;
            pos[r].push_back(i);
        }
    tot=(n-1)/sz+1;
    for(register int i=1;i<=tot;i++)
        pre_work(i);
    int x=0;    
    for(register int i=1;i<=m;i++)
        {
            int a,b;
            a=read();b=read();
            a=(a+x-1)%n+1;b=(b+x-1)%n+1;
            x=query((a<b?a:b),(a>b?a:b));
            printf("%d\n",x);
        }                  
    return 0;
}
View Code

猜你喜欢

转载自www.cnblogs.com/nanjolno/p/9344868.html