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]$,更新答案。

当块长取$N^\frac{2}{3}$,时间复杂度$O(N^\frac{5}{3})$级。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ull unsigned long long
#define ll long long
#define R register int
#define pause (for(R i=1;i<=10000000000;++i))
#define OUT freopen("out.out","w",stdout);
using namespace std;
namespace Fread {
    static char B[1<<15],*S=B,*D=B;
    #define getchar() (S==D&&(D=(S=B)+fread(B,1,1<<15,stdin),S==D)?EOF:*S++)
    inline int g() {
        R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch==‘-‘?-1:fix;
        do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
    } inline bool isempty(const char& ch) {return ch<=36||ch>=127;}
    inline void gs(char* s) {register char ch; while(isempty(ch=getchar())); do *s++=ch; while(!isempty(ch=getchar()));}
}using Fread::g; using Fread::gs;
const int N=40010,M=37; int n,m,tot,T,lst;
int f[M][M][N],h[M][M],vl[N],a[N],b[N],pos[N];
inline void PRE() { R mx=0,ans=0;
    for(R i=1;i<=n;++i) pos[i]=(i-1)/T+1;
    for(R j=1,L=pos[n];j<=L;++j,mx=0,ans=0) for(R t=j;t<=L;++t) {
        memcpy(f[j][t],f[j][t-1],sizeof(f[j][t-1]));
        for(R i=(t-1)*T+1,LL=min(t*T,n);i<=LL;++i) ++f[j][t][a[i]];
        for(R i=tot;i;--i) if(f[j][t][i]>=mx) mx=f[j][t][i],ans=i;
        h[j][t]=ans;
    }
}
signed main() {
#ifdef JACK
    freopen("NOIPAK++.in","r",stdin);
    OUT;
#endif
    n=g(),m=g(),T=n/pow(n,1.0/3);
    for(R i=1;i<=n;++i) a[i]=g(); memcpy(b,a,sizeof(a));
    sort(b+1,b+n+1); tot=unique(b+1,b+n+1)-b-1;
    for(R i=1;i<=n;++i) a[i]=lower_bound(b+1,b+tot+1,a[i])-b;
    memcpy(vl,b,sizeof(int)*(tot+1)); PRE();
    for(R i=1,l,r;i<=m;++i) { R mx=0,ans=0;
        l=(g()+lst-1)%n+1,r=(g()+lst-1)%n+1; l>r?swap(l,r):(void)0;
        R p=pos[l]+1,q=pos[r]-1; ans=h[p][q],mx=f[p][q][ans];
        if(pos[l]==pos[r]) {
            for(R i=l;i<=r;++i) { ++f[p][q][a[i]];
                if(f[p][q][a[i]]>mx||(f[p][q][a[i]]==mx&&a[i]<ans))    mx=f[p][q][a[i]],ans=a[i];
            } for(R i=l;i<=r;++i) --f[p][q][a[i]];
        } else { ans=h[p][q],mx=f[p][q][ans];
            for(R i=l,L=pos[l]*T;i<=L;++i) { ++f[p][q][a[i]];
                if(f[p][q][a[i]]>mx||(f[p][q][a[i]]==mx&&a[i]<ans))    mx=f[p][q][a[i]],ans=a[i];
            } for(R i=(pos[r]-1)*T+1;i<=r;++i) { ++f[p][q][a[i]];
                if(f[p][q][a[i]]>mx||(f[p][q][a[i]]==mx&&a[i]<ans))    mx=f[p][q][a[i]],ans=a[i];
            } for(R i=l,L=pos[l]*T;i<=L;++i) --f[p][q][a[i]];
            for(R i=(pos[r]-1)*T+1;i<=r;++i) --f[p][q][a[i]];
        } printf("%d\n",lst=vl[ans]);
    }
}

第二种是预处理出所有可能的块组成的区间中（第$i$块到第$j$块）的答案$f[i][j]$,并且拿一个$vector$存每个数$vl$出现的位置$s[vl][1...n]$。

答案初始化为整块的答案，然后对于散块中的每个数$vl$,在$s[vl]$中二分出$[l,r]$的最小和最大的位置的下标，相减就是$[l,r]$有多少个$vl$,然后更新答案。

当块长取$\sqrt{\frac{N}{logN}}$,时间复杂度$O(N\sqrt{NlogN})$。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ull unsigned long long
#define ll long long
#define R register int
#define pause (for(R i=1;i<=10000000000;++i))
#define OUT freopen("out.out","w",stdout);
using namespace std;
namespace Fread {
    static char B[1<<15],*S=B,*D=B;
    #define getchar() (S==D&&(D=(S=B)+fread(B,1,1<<15,stdin),S==D)?EOF:*S++)
    inline int g() {
        R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch==‘-‘?-1:fix;
        do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
    } inline bool isempty(const char& ch) {return ch<=36||ch>=127;}
    inline void gs(char* s) {register char ch; while(isempty(ch=getchar())); do *s++=ch; while(!isempty(ch=getchar()));}
}using Fread::g; using Fread::gs;
map<int,int> mp;
const int N=40010; int n,m,T,tot,lst;
vector<int> s[N];
#define pb push_back
int f[10010][10010],cnt[N],p[N],pos[N],a[N],b[N],vl[N];
inline void PRE(int p) { R ans=0,mx=0; memset(cnt,0,sizeof(cnt));
    for(R t=p,lim=pos[n];t<=lim;++t) {
        for(R i=(t-1)*T+1,lim=min(n,t*T);i<=lim;++i) {
            if(++cnt[a[i]]>mx||(cnt[a[i]]==mx&&a[i]<ans)) mx=cnt[a[i]],ans=a[i];
        } f[p][t]=ans;
    }
}
inline int calc(int l,int r,int x) {return upper_bound(s[x].begin(),s[x].end(),r)-lower_bound(s[x].begin(),s[x].end(),l);}
inline int solve(int l,int r) { R mx=0,ret=0;
    if(pos[l]==pos[r]) { memset(cnt,0,sizeof(cnt));
        for(R i=l;i<=r;++i) if(++cnt[a[i]]>mx||(cnt[a[i]]==mx&&a[i]<ret)) ret=a[i],mx=cnt[a[i]];
    } else { ret=f[pos[l]+1][pos[r]-1],mx=calc(l,r,ret);
        for(R i=l,lim=pos[l]*T;i<=lim;++i) { R t=calc(l,r,a[i]);
            if(t>mx||(t==mx&&a[i]<ret)) ret=a[i],mx=t;
        } for(R i=(pos[r]-1)*T+1;i<=r;++i) { R t=calc(l,r,a[i]);
            if(t>mx||(t==mx&&a[i]<ret)) ret=a[i],mx=t;
        }
    }
    return ret;
}
signed main() {
#ifdef JACK
    freopen("NOIPAK++.in","r",stdin);
    OUT;
#endif
    n=g(),m=g();//,T=n/sqrt(n*log2(n));
    T=qpow(n,1.0/4);
    for(R i=1;i<=n;++i) a[i]=g();
    memcpy(b,a,sizeof(a)); sort(b+1,b+n+1);
    tot=unique(b+1,b+n+1)-b-1; memcpy(vl,b,sizeof(int)*(tot+1));
    for(R i=1;i<=n;++i) a[i]=lower_bound(b+1,b+tot+1,a[i])-b,s[a[i]].pb(i);
    for(R i=1;i<=n;++i) pos[i]=(i-1)/T+1;
    for(R i=1;i<=pos[n];++i) PRE(i);
    for(R i=1,l,r;i<=m;++i) {
        l=(g()+lst-1)%n+1,r=(g()+lst-1)%n+1; l>r?swap(l,r):(void)0;
        printf("%d\n",lst=vl[solve(l,r)]);
    }
}

第三种按理说是复杂度最优秀的，但是跑的不是很快$qwq$。

同第二种，预处理出所有可能的块组成的区间中（第i块到第j块）的答案$f[i][j]$,并且拿一个$vector$存每个数$vl$出现的位置$s[vl][1...n]$。

然后预处理$a[i]$是整个数列中的第几个$a[i]$,出第$1$到第$i$块中最靠后的$vl$是第几个$vl$，记为$d[i][vl]$，预处理出第$n$块到第$i$块中最靠前的$vl$是第几个$vl$,记为$h[i][vl]$。

对于左边散块中的$vl$，查一下$d[pos[r]-1][vl]$,和右边散块中是否有更靠后的$vl$(可以事先用一个数组存起来)，然后可以求出这个区间有多少个$vl$(每一个$vl$是第几个$vl$已经知道了),更新答案。

当块长取$\sqrt{n}$时，时间复杂度为O(N\sqrt{N})级(不知道推没推错)。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ull unsigned long long
#define ll long long
#define R register int
#define pause (for(R i=1;i<=10000000000;++i))
#define OUT freopen("out.out","w",stdout);
using namespace std;
namespace Fread {
    static char B[1<<15],*S=B,*D=B;
    #define getchar() (S==D&&(D=(S=B)+fread(B,1,1<<15,stdin),S==D)?EOF:*S++)
    inline int g() {
        R ret=0,fix=1; register char ch; while(!isdigit(ch=getchar())) fix=ch==‘-‘?-1:fix;
        do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret*fix;
    } inline bool isempty(const char& ch) {return ch<=36||ch>=127;}
    inline void gs(char* s) {register char ch; while(isempty(ch=getchar())); do *s++=ch; while(!isempty(ch=getchar()));}
}using Fread::g; using Fread::gs;
const int N=40010,M=1000; int n,m,T,tot,lst;
vector<int> s[N];
#define pb push_back
int f[M][M],cnt[N],P[N],pos[N],a[N],b[N],vl[N],h[M][N],d[M][N],c[M][M];
inline void PRE(int p) { R ans=0,mx=0; memset(cnt,0,sizeof(cnt));
    for(R t=p,lim=pos[n];t<=lim;++t) {
        for(R i=(t-1)*T+1,lim=min(n,t*T);i<=lim;++i) {
            if(++cnt[a[i]]>mx||(cnt[a[i]]==mx&&a[i]<ans)) mx=cnt[a[i]],ans=a[i];
        } f[p][t]=ans,c[p][t]=mx;
    }
}
inline int solve(int l,int r) { R mx=0,ret=0,p=pos[l]+1,q=pos[r]-1;
    if(pos[l]==pos[r]) { memset(cnt,0,sizeof(cnt));
        for(R i=l;i<=r;++i) if(++cnt[a[i]]>mx||(cnt[a[i]]==mx&&a[i]<ret)) ret=a[i],mx=cnt[a[i]];
    } else { ret=f[p][q],mx=c[p][q]; memset(cnt,0x3f,sizeof(cnt));
        for(R i=l,L=pos[l]*T;i<=L;++i) if(cnt[a[i]]==0x3f3f3f3f) cnt[a[i]]=P[i];
        for(R i=q*T+1;i<=r;++i) {
            R tmp=P[i]+1-min(cnt[a[i]],h[p][a[i]]);
            if(tmp>mx||(tmp==mx&&a[i]<ret)) ret=a[i],mx=tmp;
            cnt[a[i]]=P[i];
        }
        for(R i=l,L=pos[l]*T;i<=L;++i) {
            R tmp=max((cnt[a[i]]==0x3f3f3f3f?0:cnt[a[i]]),d[q][a[i]])-P[i]+1;
            if(tmp>mx||(tmp==mx&&a[i]<ret)) ret=a[i],mx=tmp;
        }
    } return ret;
}
signed main() {
#ifdef JACK
    freopen("NOIPAK++.in","r",stdin);
    OUT;
#endif
    n=g(),m=g(); T=pow(n,1/2.3);//好像更小一点更快（也不是越小越快）
    for(R i=1;i<=n;++i) a[i]=g();
    memcpy(b,a,sizeof(a)); sort(b+1,b+n+1);
    tot=unique(b+1,b+n+1)-b-1; memcpy(vl,b,sizeof(int)*(tot+1));
    for(R i=1;i<=n;++i) a[i]=lower_bound(b+1,b+tot+1,a[i])-b,s[a[i]].pb(i);
    for(R i=1;i<=n;++i) pos[i]=(i-1)/T+1; for(R i=1;i<=n;++i) P[i]=++cnt[a[i]];
    memset(cnt,0,sizeof(cnt)); memset(h[pos[n]+1],0x3f,sizeof(h[pos[n]+1]));
    for(R t=pos[n];t;--t) { memcpy(h[t],h[t+1],sizeof(h[t+1]));
        for(R i=min(n,t*T);i>(t-1)*T;--i) h[t][a[i]]=P[i];
    } for(R t=1;t<=pos[n];++t) { memcpy(d[t],d[t-1],sizeof(d[t-1]));
        for(R i=(t-1)*T+1,L=t*T;i<=L;++i) d[t][a[i]]=P[i];
    } for(R i=1;i<=pos[n];++i) PRE(i);
    for(R i=1,l,r;i<=m;++i) {
        l=(g()+lst-1)%n+1,r=(g()+lst-1)%n+1; l>r?swap(l,r):(void)0;
        printf("%d\n",lst=vl[solve(l,r)]);
    }
}

2019.06.28

原文地址：https://www.cnblogs.com/Jackpei/p/11105253.html

时间： 2024-11-06 03:34:37

P4168 [Violet]蒲公英 分块的相关文章

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

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

Luogu P4168 [Violet]蒲公英

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

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

[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

【BZOJ2724】[Violet 6]蒲公英 分块+二分

[BZOJ2724][Violet 6]蒲公英 Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1 Output Sample Input 6 3 1 2 3 2 1 2 1 5 3 6 1 5 Sample Output 1 2 1 HINT 修正下: n <= 40000, m <= 50000 题解:分块还是练脑子啊~ 结论:一个区间的众数要么是区间中一个块的众数,要么是块外的任意

【bzoj2724】[Violet 6]蒲公英 分块+STL-vector

题目描述 输入 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1 输出 样例输入 6 3 1 2 3 2 1 2 1 5 3 6 1 5 样例输出 1 2 1 题解 分块+STL-vector 一个显而易见的结论:区间众数一定是一段连续的块的众数或块外的数,证明略(逃 先把数据离散化,然后分块预处理出f[i][j],表示从块i到块j的众数位置.具体实现的话直接开个桶存一下就好了. 然后考虑询问,整块的直接拿出来求一下出现次

[日常摸鱼]bzoj2724蒲公英-分块

区间众数经典题~ http://begin.lydsy.com/JudgeOnline/problem.php?id=4839这里可以提交~ 题意大概就是没有修改的询问区间众数,如果有一样的输出最小的,强制在线,$n \leq 4*10^4,a_i \leq 10^9$. log数据结构脑补一遍好像没什么可以做的,数据范围我们可以分块! 不过分块之前肯定要离散化一下,而且还要保存离散化前的数据(因为要回答的是出现最多的数),离散化的方法在上一篇博客里面~ 假设分成$L$块,每块大小$s=\lfl

题解【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