Luogu 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[<<],*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=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
} inline bool isempty(const char& ch) {return ch<=||ch>=;}
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=,M=; 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=,ans=;
for(R i=;i<=n;++i) pos[i]=(i-)/T+;
for(R j=,L=pos[n];j<=L;++j,mx=,ans=) for(R t=j;t<=L;++t) {
memcpy(f[j][t],f[j][t-],sizeof(f[j][t-]));
for(R i=(t-)*T+,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/);
for(R i=;i<=n;++i) a[i]=g(); memcpy(b,a,sizeof(a));
sort(b+,b+n+); tot=unique(b+,b+n+)-b-;
for(R i=;i<=n;++i) a[i]=lower_bound(b+,b+tot+,a[i])-b;
memcpy(vl,b,sizeof(int)*(tot+)); PRE();
for(R i=,l,r;i<=m;++i) { R mx=,ans=;
l=(g()+lst-)%n+,r=(g()+lst-)%n+; l>r?swap(l,r):(void);
R p=pos[l]+,q=pos[r]-; 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]-)*T+;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]-)*T+;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[<<],*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=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
} inline bool isempty(const char& ch) {return ch<=||ch>=;}
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=; int n,m,T,tot,lst;
vector<int> s[N];
#define pb push_back
int f[][],cnt[N],p[N],pos[N],a[N],b[N],vl[N];
inline void PRE(int p) { R ans=,mx=; memset(cnt,,sizeof(cnt));
for(R t=p,lim=pos[n];t<=lim;++t) {
for(R i=(t-)*T+,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=,ret=;
if(pos[l]==pos[r]) { memset(cnt,,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]+][pos[r]-],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]-)*T+;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/);
for(R i=;i<=n;++i) a[i]=g();
memcpy(b,a,sizeof(a)); sort(b+,b+n+);
tot=unique(b+,b+n+)-b-; memcpy(vl,b,sizeof(int)*(tot+));
for(R i=;i<=n;++i) a[i]=lower_bound(b+,b+tot+,a[i])-b,s[a[i]].pb(i);
for(R i=;i<=n;++i) pos[i]=(i-)/T+;
for(R i=;i<=pos[n];++i) PRE(i);
for(R i=,l,r;i<=m;++i) {
l=(g()+lst-)%n+,r=(g()+lst-)%n+; l>r?swap(l,r):(void);
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[<<],*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=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
} inline bool isempty(const char& ch) {return ch<=||ch>=;}
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=,M=; 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=,mx=; memset(cnt,,sizeof(cnt));
for(R t=p,lim=pos[n];t<=lim;++t) {
for(R i=(t-)*T+,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=,ret=,p=pos[l]+,q=pos[r]-;
if(pos[l]==pos[r]) { memset(cnt,,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+;i<=r;++i) {
R tmp=P[i]+-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?:cnt[a[i]]),d[q][a[i]])-P[i]+;
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,/2.3);//好像更小一点更快（也不是越小越快）
for(R i=;i<=n;++i) a[i]=g();
memcpy(b,a,sizeof(a)); sort(b+,b+n+);
tot=unique(b+,b+n+)-b-; memcpy(vl,b,sizeof(int)*(tot+));
for(R i=;i<=n;++i) a[i]=lower_bound(b+,b+tot+,a[i])-b,s[a[i]].pb(i);
for(R i=;i<=n;++i) pos[i]=(i-)/T+; for(R i=;i<=n;++i) P[i]=++cnt[a[i]];
memset(cnt,,sizeof(cnt)); memset(h[pos[n]+],0x3f,sizeof(h[pos[n]+]));
for(R t=pos[n];t;--t) { memcpy(h[t],h[t+],sizeof(h[t+]));
for(R i=min(n,t*T);i>(t-)*T;--i) h[t][a[i]]=P[i];
} for(R t=;t<=pos[n];++t) { memcpy(d[t],d[t-],sizeof(d[t-]));
for(R i=(t-)*T+,L=t*T;i<=L;++i) d[t][a[i]]=P[i];
} for(R i=;i<=pos[n];++i) PRE(i);
for(R i=,l,r;i<=m;++i) {
l=(g()+lst-)%n+,r=(g()+lst-)%n+; l>r?swap(l,r):(void);
printf("%d\n",lst=vl[solve(l,r)]);
}
}
2019.06.28
Luogu P4168 [Violet]蒲公英 分块的更多相关文章
- luogu P4168 [Violet]蒲公英
嘟嘟嘟 分块经典题竟然是一道黑题…… 分块求区间众数的大体思想是对于询问区间[L, R],预处理出这中间的整块的众数,然后统计两边零散的数在[L, R]中出现的次数,最后取出现次数最多且最小的数. 因 ...
- 洛谷 P4168 [Violet]蒲公英 解题报告
P4168 [Violet]蒲公英 题目背景 亲爱的哥哥: 你在那个城市里面过得好吗? 我在家里面最近很开心呢.昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多 ...
- p4168 [Violet]蒲公英（分块）
区间众数的重题 和数列分块入门9双倍经验还是挺好的 然后开O2水过 好像有不带log的写法啊 之后在补就是咕咕咕 // luogu-judger-enable-o2 #include <cstd ...
- [Violet]蒲公英 分块
发现写算法专题老是写不动,,,, 所以就先把我在luogu上的题解搬过来吧! 题目大意:查询区间众数,无修改,强制在线 乍一看是一道恐怖的题,仔细一看发现并没有那么难: 大致思路是这样的,首先我们要充 ...
- P4168 [Violet]蒲公英 区间众数
$ \color{#0066ff}{ 题目描述 }$ 在乡下的小路旁种着许多蒲公英,而我们的问题正是与这些蒲公英有关. 为了简化起见,我们把所有的蒲公英看成一个长度为n的序列 \((a_1,a_2.. ...
- P4168 [Violet]蒲公英
神仙分块题?其实还是很简单的,res[i][j]表示第i块到第j块的众数,然后再用sum[i][j]表示前i块中j这个种类出现的次数,然后分块瞎搞就行了,感觉我写的十分简洁,好评( //author ...
- BZOJ2724 [Violet]蒲公英 分块
题目描述 经典区间众数题目 然而是权限题,所以题目链接放Luogu的 题解 因为太菜所以只会$O(n*\sqrt{n}+n*\sqrt{n}*log(n))$的做法 就是那种要用二分的,并不会clj那 ...
- [洛谷P4168][Violet]蒲公英
题目大意:有$n(n\leqslant4\times10^4)$个数,$m(m\leqslant5\times10^4)$个询问,每次问区间$[l,r]$内的众数,若相同输出最小的,强制在线. 题解: ...
- 洛谷 P4168 [Violet] 蒲公英
历尽千辛万苦终于AC了这道题目... 我们考虑1个区间\([l,r]\), 被其完整包含的块的区间为\([L,R]\) 那么众数的来源? 1.\([l,L)\)或\((R,r]\)中出现的数字 2.\ ...
随机推荐
- 【LOJ】#3103. 「JSOI2019」节日庆典
LOJ#3103. 「JSOI2019」节日庆典 能当最小位置的值一定是一个最小后缀,而有用的最小后缀不超过\(\log n\)个 为什么不超过\(\log n\)个,看了一下zsy的博客.. 假如\ ...
- Win7原装ISO镜像封装USB3.0&网卡驱动
Win7原装ISO镜像封装USB3.0&网卡驱动 最新购买的电脑是Windows10系统,想装回Windows7,但是装Windows7发现网络适配器没出现,如果没有USB2.0接口,US ...
- sklearn数据集
数据集划分: 机器学习一般的数据集会划分为两个部分 训练数据: 用于训练,构建模型 测试数据: 在模型检验时使用,用于评估模型是否有效 sklearn数据集划分API: 代码示例文末! scikit- ...
- python如何切割字符串
python字符串的分割方法如下 str.split():字符串分割函数 通过指定分隔符对字符串进行切片,并返回分割后的字符串列表. 语法: str.split(s, num)[n] 参数说明: s: ...
- 在win下开发的项目怎么迁移到linux下面才能正常运行?
我可以直接拷贝项目目录到linux下面直接操作吗? 答案: 可以,咋可能??? 为啥??? win开发直接拷贝过去,你不凉谁凉了,我以前也同样的单纯,如果你项目里用的绝对路径! 那恭喜你,你凉了,清楚 ...
- ingress之tls和path使用
ingress tls 上节课给大家展示了 traefik 的安装使用以及简单的 ingress 的配置方法,这节课我们来学习一下 ingress tls 以及 path 路径在 ingress 对象 ...
- NIO总结-----Buffer
Java NIO中的Buffer用于和NIO通道进行交互.如你所知,数据是从通道读入缓冲区,从缓冲区写入到通道中的. 缓冲区本质上是一块可以写入数据,然后可以从中读取数据的内存.这块内存被包装成NIO ...
- poj 2891 模数不互质的中国剩余定理
Strange Way to Express Integers Description Elina is reading a book written by Rujia Liu, which intr ...
- 4.性能下降原因和常见的Join查询
性能下降 SQL慢,执行时间长,等待时间长 1.查询语句写的烂 2.索引失效 单值索引失效 和 复合索引失效 3.关联查询太多join(设计缺陷或不得已的需求) 4.服务器调优及各个参数设置(缓冲.线 ...
- jQuery EasyUI 数据网格
jQuery EasyUI 数据网格 - 转换 HTML 表格为数据网格 本节将介绍jQuery EasyUI数据网格的运用,主要内容为如何将HTML表格转换为数据网格. 本实例演示如何转换表格(ta ...