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[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 JACKfreopen("NOIPAK++.in","r",stdin);OUT;
#endifn=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 JACKfreopen("NOIPAK++.in","r",stdin);OUT;
#endifn=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 JACKfreopen("NOIPAK++.in","r",stdin);OUT;
#endifn=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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/483032.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

一文掌握明年物联网传感器市场!2022中国AIoT产业全景图谱报告新鲜出炉

来源&#xff1a;传感器专家网物联网产业是传感器应用最广泛的领域之一&#xff0c;研发新型传感器&#xff0c;做传感器市场&#xff0c;都不能不考虑物联网产业的需求。2021年随着COVID-19 病毒的存在已常态化&#xff0c;防疫抗疫进入了拉锯阶段&#xff1b;波及全球的“芯片…

spring学习笔记06-spring整合junit(出现的问题,解决的思路)

spring学习笔记06-spring整合junit&#xff08;出现的问题&#xff0c;解决的思路&#xff09; 文章目录spring学习笔记06-spring整合junit&#xff08;出现的问题&#xff0c;解决的思路&#xff09;3.1测试类中的问题和解决思路3.1.1 问题3.2.1 第一步&#xff1a;pom.xml文件…

phpstorm 不能自动打开上次的历史文件

问题产生的原因&#xff1a;可能是电脑非正常关机&#xff0c;导致phpstrom 无法正常关闭。 一开始我以为是配置上那里出现了问题&#xff0c;所以直接就把配置删除了&#xff0c;然后生成了默认配置。发现还是无法解决问题。然后一个个配置查阅过去。因为我的ide是英文&#x…

EUV光刻机内部揭秘!

转载自: ittbank来源&#xff1a;半导体行业观察PatrickWhelan正在透过他的洁净室服面板凝视着事情的进展。在他面前是一块闪闪发光的玻璃&#xff0c;大约有一个烤箱那么大&#xff0c;上面刻有许多挖出的部分以减轻重量&#xff0c;看起来像一个外星图腾。Whelan 的团队正在将…

数据结构与算法——搜索

文章目录1.内容概述2.岛屿数量2.1 题目描述2.2 DFS深度搜索算法思路2.3 BFS宽度搜索算法思路2.4 C代码实现3.单词接龙3.1 题目描述3.2 算法思路3.3 C代码实现4.单词接龙 II4.1 题目描述4.2 算法思路5.火柴拼正方形5.1 题目描述5.2 算法思路5.3 代码实现5.4 算法思路25.5 代码实…

mybatis学习笔记-01什么是mybatis

该视频为狂神说java视频配套笔记&#xff08;博主自己手打223&#xff0c;日后做个参考223&#xff09;&#xff0c;b站连接&#xff1a;Mybatis最新完整教程IDEA版【通俗易懂】-01-什么是mybatis 文章目录1.简介1.1是什么1.2持久层1.3持久层1.4为什么需要mybatis&#xff1f;1…

小程序基础能力~网络

&#xff08;1&#xff09;网络-使用说明 网络 在小程序/小游戏中使用网络相关的 API 时&#xff0c;需要注意下列问题&#xff0c;请开发者提前了解。 1. 服务器域名配置 每个微信小程序需要事先设置通讯域名&#xff0c;小程序只可以跟指定的域名与进行网络通信。包括普通 HT…

mybatis学习笔记-02-第一个mybatis程序

该视频为狂神说java视频配套笔记&#xff08;博主自己手打223&#xff0c;日后做个参考223&#xff09;&#xff0c;b站连接&#xff1a;Mybatis最新完整教程IDEA版【通俗易懂】-02-第一个mybatis程序) 文章目录2.第一个mybatis程序2.1搭建环境2.2 创建一个模块2.3编写代码2.4测…

徐波所长专访 | 人工智能:从“作坊式”走向“工业化”新时代

来源&#xff1a;人民网人工智能创新不断 “一体两翼”快速发展人民网&#xff1a;当前&#xff0c;人工智能技术创新不断&#xff0c;应用层出不穷。它究竟走到了哪一步&#xff1f;能否谈谈您是如何看待我国人工智能技术发展现状的&#xff1f;徐波&#xff1a;人工智能是一个…

linux——线程(1)

文章目录1.线程概念1.1 什么是线程1.2 线程和进程区别1.3 线程实现原理1.4 三级映射1.5 线程共享资源1.6 线程非共享资源1.7 线程优、缺点2.线程控制原语2.1 pthread_self 函数2.2 pthread_create 函数3.线程与共享3.1 线程共享全局变量4.线程退出4.1 pthread_exit 函数1.线程概…

mybatis学习笔记-03-CRUD操作

该视频为狂神说java视频配套笔记&#xff08;博主自己手打223&#xff0c;日后做个参考223&#xff09;&#xff0c;b站连接&#xff1a;Mybatis最新完整教程IDEA版【通俗易懂】-03-CRUD操作) 3.CRUD 1.namespace ​ namespace中的包名要和 Dao/mapper 接口的包名一致 2.se…

Redis压缩列表原理与应用分析

摘要 Redis是一款著名的key-value内存数据库软件&#xff0c;同时也是一款卓越的数据结构服务软件。它支持字符串、列表、哈希表、集合、有序集合五种数据结构类型&#xff0c;同时每种数据结构类型针对不同的应用场景又支持不同的编码方式。这篇文章主要介绍压缩列表编码&…

Science公布2021年度十大科学突破,AI这项前所未有的突破上榜

来源&#xff1a; 学术头条撰文&#xff1a;朱哼哼编审&#xff1a;王哈哈排版&#xff1a;李雪薇2021 年 11 月 17 日&#xff0c;Science 杂志公布了 2021 年的年度科学突破榜单&#xff0c;AlphaFold 和 RoseTTA-fold 两种基于人工智能预测蛋白质结构的技术位列榜首。除此之…

linux——线程(2)

文章目录1.pthread_join 函数2.pthread_detach 函数3.pthread_cancel 函数4.终止线程方式5.控制原语对比6.线程属性7.线程属性初始化8.线程的分离状态9.线程使用注意事项1.pthread_join 函数 阻塞等待线程退出&#xff0c;获取线程退出状态 其作用&#xff0c;对应进程中 waitp…

mybatis学习笔记-04-常见错误排查

4.错误排查 1.标签问题 mapper.xml中的namespace&#xff0c;中间是用小圆点隔开的 resource绑定 mapper&#xff0c;需要使用路径!&#xff08;使用/隔开&#xff09; 3.程序配置文件必须符合要求 4.NullPointerException &#xff1a;没有注册到资源 5.输出的xml文件中…

linux——线程通信(1)

文章目录1.同步概念2.线程同步3.数据混乱原因&#xff1a;4.互斥量 mutex5.mutex 主要应用函数&#xff1a;5.1 pthread_mutex_init 函数5.2 pthread_mutex_destroy 函数5.3 pthread_mutex_lock 函数5.4 pthread_mutex_unlock 函数5.5 pthread_mutex_trylock 函数5.6 加锁与解锁…

OpenAI教GPT-3学会上网,「全知全能」的AI模型上线了

来源&#xff1a;机器学习研究组订阅它被命名为 WebGPT&#xff0c;OpenAI 认为浏览网页的方式提高了 AI 解答问题的准确性。如果 AI 学会上网&#xff0c;那么它就拥有了无限获取知识的方式&#xff0c;之后会发生什么就不太好预测了。于是著名 AI 研究机构 OpenAI 教那个开启…

201671030111 李蓉 实验十四 团队项目评审课程学习总结

项目内容这个作业属于哪个课程软件工程这个作业的要求在哪里实验十四 团队项目评审&课程学习总结作业学习目标掌握软件项目评审会流程&#xff0c;反思总结课程学习内容。任务一&#xff1a;结合本学期课程学习内容&#xff0c;对比《实验一 软件工程准备》的任务5你所提的…

linux——线程通信(2)

文章目录1.条件变量1.1 条件变量函数&#xff1a;1.2 pthread_cond_init 函数1.3 pthread_cond_destroy 函数1.4 pthread_cond_wait 函数1.5 pthread_cond_timedwait 函数1.6 pthread_cond_signal 函数1.7 pthread_cond_broadcast 函数2.生产者消费者模型3.条件变量的优点&…

Hadoop(2):常见的MapReduce[在Ubuntu中运行!]

1 以词频统计为例子介绍 mapreduce怎么写出来的 弄清楚MapReduce的各个过程&#xff1a; 将文件输入后&#xff0c;返回的<k1,v1>代表的含义是&#xff1a;k1表示偏移量&#xff0c;即v1的第一个字母在文件中的索引&#xff08;从0开始数的&#xff09;&#xff1b;v1表…