[violet]蒲公英题解

前几天刚学习了分块,感觉这道题用分块求解的方式挺巧妙的

既然用的是分块,那么肯定是两端暴力求解,中间要快速地处理每个块

首先我们要得到一个结论,最后求出的这一个众数必定为中间块的众数或者是两端的任意一个数,那么我们用\(nu[i][j]\)来表示第\(i\)个块到第\(j\)个块的众数,我们可以用用\(O(n\sqrt{n})\)的时间复杂度,先枚举每个块,然后枚举后面的点来求出\(nu[i][j]\),放一下这一段的代码

void pre(int x)
{
    memset(cnt,0,sizeof(cnt));
    int maxx=0,num1=0;
    for(int i=(x-1)*m+1;i<=n;i++)
    {
        cnt[a[i]]++;
        if(cnt[a[i]]>maxx||(cnt[a[i]]==maxx&&a[i]<num1))
            maxx=cnt[a[i]],num1=a[i];
        nu[x][pos[i]]=num1;
    }
}

接下来就是每一次的询问,先假定中间的块的众数为整个区间的众数,并求出在区间内的个数,再枚举两端的每一个点,求出在区间内的个数,如果比之前的多那么就更新

至于求一个数在区间内的个数,我们可以巧用\(vector\),记录每个数字出现的位置,然后用upper_bound和lower_bound求出在区间内的个数

记得离散化!

#include<bits/stdc++.h>
using namespace std;
int n,a[100003],num[100003],id,pos[100003],l,r,su,cnt[100003],nu[3003][3003],res,m,b[100003],las,o;
vector<int>l1[100003];
map<int,int>mp;
int ask1(int ll,int rr,int x)
{
    return upper_bound(l1[x].begin(),l1[x].end(),rr)-lower_bound(l1[x].begin(),l1[x].end(),ll); //求出一个数在区间内的个数
}
void pre(int x) //预处理nu[i][j]
{
    memset(cnt,0,sizeof(cnt));
    int maxx=0,num1=0;
    for(int i=(x-1)*m+1;i<=n;i++)
    {
        cnt[a[i]]++;
        if(cnt[a[i]]>maxx||(cnt[a[i]]==maxx&&a[i]<num1))
            maxx=cnt[a[i]],num1=a[i];
        nu[x][pos[i]]=num1;
    }
}
int ask(int l1,int r1)
{
    int maxx=0,num1=0;
    if(pos[l1]+1<pos[r1])
        maxx=ask1(l1,r1,nu[pos[l1]+1][pos[r1]-1]),num1=nu[pos[l1]+1][pos[r1]-1]; //假定众数为中间的块的众数
    for(int i=l1;i<=min(r1,pos[l1]*m);i++) //枚举左端的块
    {
        res=ask1(l1,r1,a[i]);
        if(res>maxx||(res==maxx&&a[i]<num1))
            maxx=res,num1=a[i];
    }
    if(pos[l1]!=pos[r1]) //枚举右端的块
        for(int i=(pos[r1]-1)*m+1;i<=r1;i++)
        {
            res=ask1(l1,r1,a[i]);
            if(res>maxx||(res==maxx&&a[i]<num1))
                maxx=res,num1=a[i];
        }
    return num[num1];
}
int main()
{
    scanf("%d%d",&n,&o);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        b[i]=a[i];
    }
    sort(b+1,b+n+1);
    for(int i=1;i<=n;i++)
        if(!mp[b[i]])
            id++,mp[b[i]]=id,num[id]=b[i]; //离散化
    for(int i=1;i<=n;i++)
        a[i]=mp[a[i]],l1[a[i]].push_back(i);
    for(int i=1;i<=id;i++)
        l1[i].push_back(n+1);
    m=pow(n,0.4);
    for(int i=1;i<=n;i++)
        pos[i]=(i-1)/m+1;
    su=pos[n];
    for(int i=1;i<=su;i++)
        pre(i);
    for(int i=1;i<=o;i++)
    {
        scanf("%d%d",&l,&r);
        l=(l+las-1)%n+1,r=(r+las-1)%n+1;
        if(l>r)
            swap(l,r);
        las=ask(l,r);
        cout<<las<<endl;
    }
    return 0;
}

原文地址:https://www.cnblogs.com/dzice/p/12195959.html

时间: 2024-07-30 03:57:30

[violet]蒲公英题解的相关文章

Luogu P4168 [Violet]蒲公英

P4168 [Violet]蒲公英 题意 题目背景 亲爱的哥哥: 你在那个城市里面过得好吗? 我在家里面最近很开心呢.昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多小朋友也被它杀掉了.我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢.不过奶奶说他是很难受的时候才做出这样的事的-- 最近村子里长出了一大片一大片的蒲公英.一刮风,这些蒲公英就能飘到好远的地方了呢.我觉得要是它们能飘到那个城市里面,让哥哥看看就好了呢! 哥哥你要快点回来哦! 爱你的妹妹 \(V

[Luogu P4168] [Violet]蒲公英 (分块)

题面 洛咕 Solution 题目要求求出区间众数,强制在线. 区间众数是一个比较尴尬的问题,我们无法用区间数据结构来处理这个问题,因为我们没法很好的合并区间众数的答案. 既然区间数据结构解决不了这个问题,我们可以考虑一下使用基于分块的算法,例如莫队. 这题用莫队非常好处理,不幸的是,这题要求强制在线. 因此我们考虑使用分块算法. 分块算法的核心在于把一整个块的信息压缩起来以便快速处理. 我们要查询一段区间的众数,我们可以考虑这样搞:对于这个区间内连续的块,我们先快速地查询这个连续的块中的众数,

题解【luogu4168 [Violet]蒲公英】

Description 给出一个长度为 \(n\) 序列 \(a\) ,\(m\) 次询问,每次询问区间 \([l,r]\) 里的众数(出现次数最多的数).若有多个,输出最小的. \(a_i \leq 10^9, n \leq 40000, m \leq 50000\),强制在线. Solution \(a_i \leq 10^9\) ,先离散化.然后 算法一:暴力 \(n ^ 2\) ,预计得分 20 : 实际得分 20 (开了 O2 直接变成 85 啥操作) 算法二:\(n \leq 400

[Violet]蒲公英

题意: 给出一个长度为 \(n\) 序列\(a\) ,\(m\)次询问,每次询问区间 \(l,r\) 里的众数(出现次数最多的数).若有多个,输出最小的. \(a_i \leq 10^9, n \leq 40000, m \leq 50000\),强制在线. 题解: 看了题解才懂的.根据https://www.cnblogs.com/acfunction/p/10051345.html hzwer给出了更巧妙的方法http://hzwer.com/3582.html \(a_i \le 10^9

[Violet]蒲公英 分块

Code: #include<cstdio> #include<string> #include<iostream> #include<algorithm> #include<vector> #include<cstring> #include<cmath> using namespace std; void setIO(string a){ freopen((a+".in").c_str(),&quo

p4168 [Violet]蒲公英(分块)

区间众数的重题 和数列分块入门9双倍经验还是挺好的 然后开O2水过 好像有不带log的写法啊 之后在补就是咕咕咕 // luogu-judger-enable-o2 #include <cstdio> #include <algorithm> #include <cstring> #include <vector> #include <map> #include <cmath> using namespace std; int m,b

P4168 [Violet]蒲公英 分块

这道题算是好好写了.写了三种方法. 有一个好像是$qwq$$N\sqrt(N)$的方法,,但是恳请大佬们帮我看看为什么这么慢$qwq$(后面的第三种) 注:$pos[i]$表示$i$属于第$pos[i]$块. 第一种是统计所有可能的块组成的区间中(第i块到第j块),每个数出现的次数,记做$f[i][j][k]$,和所有可能的块组成的区间的答案,记做$h[i][j]$. 然后每次先把整块的答案作为初始答案,然后对于散块中的每个值$vl$,暴力修改对应的$f[i][j][vl]$,更新答案. 当块长

分块简单入门

分块简单入门 按 树状数组虽超级快,但是很僵硬,不灵活:线段树虽快,但是却不直观,代码量大:所以,速度较慢但直观形象.代码量小的分块大法不实为线段树的替代品. 网络上关于分块的教程不知道为什么很少,虽然早有hzwer大神的分块九讲,但是还是少了入门级详解教程.此篇将分为三个阶段,保证初学者在有意识地阅读后基本掌握分块. 1.简单的入门题 问题引入 给定长度为N(N<=1e5)的正数数列A,然后输入Q(Q<=1e5)行操作命令,指令形如Q l r,表示统计数列中第l~r个数中的最大值 思路 首先

BZOJ2724: [Violet 6]蒲公英

2724: [Violet 6]蒲公英 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 795  Solved: 248[Submit][Status] Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1 Output Sample Input Sample Output HINT 修正下: n <= 40000, m <= 50000 S