NOI模拟(3.3)螺旋序列(出题人一定是月厨)
Description
S也想寻求真正的智慧,然而由于“抑制力”的存在,她必须先解决一系列询
问。
有一个长度为n的序列a,一个长度为m序列b被称为螺旋序列当且仅当
b1=bm且对于1<=i<=m有bi<=b1。
S需要回答q个询问,每个询问用l,r两个参数描述,表示询问区间[l,r]的最长
连续子螺旋序列的长度。
Input
第一行两个整数n,q,表示序列长度和询问数。
第二行n个整数ai表示序列a。
以下q行,每行两个整数l,r表示一次询问。
Output
对每次询问输出一行一个整数表示最大连续螺旋序列的长度。
Sample Input
5 3
3 2 1 2 3
1 4
2 5
1 5
Sample Output
3
3
5
Data Constraint
本题采用捆绑测试,只有通过一个子任务的全部数据才能得到该子任务分
数,否则不得分。
子任务1(20分):n,m<=10
子任务2(20分):n,m<=1000
子任务3(20分):n,m<=2*10^5,1<=ai<=10
子任务4(40分):n,m<=5*10^5;|ai|<=10^9
Solution
30分的裸暴力
离散化相等的数字,做一个链表,st表处理区间最大值,对每个询问都暴力去跳
#include <map>
#include <vector>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rg register template<class T> inline T dma(rg T x,rg T y) {return x>y?x:y;}
template<class T> inline void read(rg T &x)
{
rg int c=getchar();rg bool b=0;
for(;c<48||c>57;c=getchar())
if(c==45)b=1;
for(x=0;c>47&&c<58;c=getchar())
x=(x<<1)+(x<<3)+c-48;
if(b)x=-x;
} const int N=500010;
std::vector<int>vc;
std::map<int,int>hc;
int n,q,seq[N],fir[N],nex[N],st[20][N],bin[20],lgg[20],ans; void sequence_table()
{
lgg[0]=-1;for(rg int i=1;i<=n;i++)lgg[i]=lgg[i>>1]+1;
for(rg int i=0;i<20;i++)bin[i]=1<<i;
for(rg int i=1;i<=n;i++)st[0][i]=seq[i];
for(rg int t=1;t<20;t++)
for(rg int i=1;i<=n;i++)
if(i+bin[t]-1<=n)st[t][i]=dma(st[t-1][i],st[t-1][i+bin[t-1]]);
} int sequence_query(rg int x,rg int y)
{
rg int t=lgg[y-x+1];
return dma(st[t][x],st[t][y-bin[t]+1]);
} void jump(rg int x,rg int lim)
{
for(rg int rig=x,lef=nex[x];lef>=lim;lef=nex[lef])
{
while(rig>lef&&sequence_query(lef,rig)>seq[lef])rig=nex[rig];
ans=dma(ans,rig-lef+1);
}
} int main()
{
freopen("sequence.in","r",stdin);
freopen("sequence.out","w",stdout);
read(n),read(q);
for(rg int i=1;i<=n;i++)
{
read(seq[i]);
vc.push_back(seq[i]);
}
std::sort(vc.begin(),vc.end());
vc.erase(std::unique(vc.begin(),vc.end()),vc.end());
for(rg int i=0;i<vc.size();i++)
hc[vc[i]]=i+1;
for(rg int i=1;i<=n;i++)seq[i]=hc[seq[i]];
memset(fir,-1,sizeof fir);
for(rg int i=1;i<=n;i++)
{
nex[i]=fir[seq[i]];
fir[seq[i]]=i;
}
sequence_table();
for(rg int x,y;q;q--)
{
ans=1;
read(x),read(y);
for(rg int i=y;i>=x;i--)
jump(i,x);
printf("%d\n",ans);
}
fclose(stdin);
fclose(stdout);
return 0;
}
出题人一定是月厨,有爱~
考虑将这些进行jump的操作序列转化成链,离线后,用链来更新询问
对于所有链长度小于s的,用树状数组区间最值+扫描线做
对于所有链长度大于s的,这样的链只有n/s条,对每条链的下标预处理每个下标i所在链上的位置,对每个询问算出最优值并更新,复杂度为n*n/s
取s=sqrt(n/logn),可以得到O(n*sqrt(n*logn))的复杂度
#include <map>
#include <vector>
#include <math.h>
#include <stdio.h>
#include <string.h>
#include <algorithm> template<class T> T mex(T x,T y) {return x>y?x:y;}
template<class T> void read(T &x)
{
int c=getchar();bool b=0;
for(;c<48||c>57;c=getchar())if(c==45)b=1;
for(x=0;c>47&&c<58;c=getchar())x=(x<<1)+(x<<3)+c-48;
if(b)x=-x;
} const int maxn=500010;
bool vis[maxn];
int n,m,lim,A[maxn],st[20][maxn],lgg[maxn],L[maxn],R[maxn],ans[maxn],nex[maxn],big[maxn],cur[maxn],tr[maxn]; std::map<int,int> pool;
std::vector<int> seq;
std::vector<std::pair<int,int> > Q[maxn]; inline void Add(int x,int c)
{
for(;x<=n;x+=x&-x)tr[x]=mex(tr[x],c);
} inline int Get(int x)
{
int ret=0;
for(;x;x-=x&-x)ret=mex(ret,tr[x]);
return ret;
} inline int maxi(int x,int y)
{
int t=lgg[++y-x];
return mex(st[t][x],st[t][y-(1<<t)]);
} int main()
{
freopen("sequence.in","r",stdin);
freopen("sequence.out","w",stdout);
read(n),read(m);
lgg[0]=-1;for(int i=1;i<=n;i++)
{
read(A[i]);
st[0][i]=A[i];
lgg[i]=lgg[i>>1]+1;
}
for(int i=1;i<=m;i++)
{
read(L[i]),read(R[i]);
Q[L[i]].push_back(std::make_pair(R[i],i));
ans[i]=1;
}
for(int t=1;t<20;t++)
for(int i=1;i+(1<<t)-1<=n;i++)
st[t][i]=mex(st[t-1][i],st[t-1][i+(1<<t-1)]);
for(int i=n,j;i;i--)
{
if(pool.count(A[i]))
{
j=pool[A[i]];
if(maxi(i,j)<=A[i])nex[i]=j;
}
pool[A[i]]=i;
}
lim=static_cast<int>(sqrt(n/log(n)/log(2)));
for(int i=1;i<=n;i++)
if(!vis[i])
{
seq.clear();
seq.push_back(0);
for(int j=i;j;j=nex[j])
{
seq.push_back(j);
vis[j]=true;
}
if(seq.size()>=lim)
{
for(int j=1,k=0;j<=n;j++)
{
for(;k<seq.size()&&j>=seq[k];k++)big[seq[k]]=1;
cur[j]=k-1;
}
for(int j=1,x,y;j<=m;j++)
if(cur[L[j]-1]+1<seq.size())
{
x=seq[cur[L[j]-1]+1];
y=seq[cur[R[j]]];
ans[j]=mex(ans[j],y-x+1);
}
}
}
for(int i=n;i;i--)
{
if(!big[i])
for(int j=i;j;j=nex[j])
Add(j,j-i+1);
for(int j=0;j<Q[i].size();j++)
ans[Q[i][j].second]=mex(ans[Q[i][j].second],Get(Q[i][j].first));
}
for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
fclose(stdin);
fclose(stdout);
return 0;
}
NOI模拟(3.3)螺旋序列(出题人一定是月厨)的更多相关文章
- Java实现 LeetCode 521 最长特殊序列 Ⅰ(出题人:“就是喜欢看你们不敢相信那么简单,又不敢提交的样子。”)
521. 最长特殊序列 Ⅰ 给定两个字符串,你需要从这两个字符串中找出最长的特殊序列.最长特殊序列定义如下:该序列为某字符串独有的最长子序列(即不能是其他字符串的子序列). 子序列可以通过删去字符串中 ...
- NOI 模拟赛 #2
得分非常惨惨,半个小时写的纯暴力 70 分竟然拿了 rank 1... 如果 OYJason 和 wxjor 在可能会被爆踩吧 嘤 T1 欧拉子图 给一个无向图,如果一个边集的导出子图是一个欧拉回路, ...
- D.出题人的手环
链接:https://ac.nowcoder.com/acm/contest/358/D 题意: 出题人的妹子送了出题人一个手环,这个手环上有 n 个珠子,每个珠子上有一个数. 有一天,出题人和妹子分 ...
- ACM_出题人这样不好吧
出题人这样不好吧 Time Limit: 2000/1000ms (Java/Others) Problem Description: 作为编协的第一次月赛,肯定是要有防AK(ALL KILL)的题目 ...
- 牛客 545C 出题人的数组 (贪心)
出题人有两个数组A,B,请你把两个数组归并起来使得$cost=\sum i c_i$最小. 归并要求原数组的数的顺序在新数组中不改变. 贪心水题 对于一段序列$A_i,A_{i+1},...,A_r$ ...
- 6.6 NOI 模拟
\(T1\)括号序列 --那是,朝思夜想也未尝得到的自由 一个比较常见的转化,考虑如何判断前一段和后一段能够拼成一个合法的括号序列 充要条件: 前半部分,'('看为\(1\), ')'看为\(-1\) ...
- 5.6 NOI模拟
\(5.6\ NOI\)模拟 明天就母亲节了,给家里打了个电话(\(lj\ hsez\)断我电话的电,在宿舍打不了,只能用教练手机打了) 其实我不是很能看到自己的\(future,\)甚至看不到高三的 ...
- 5.4 NOI模拟
\(5.4\ NOI\)模拟 \(T1\) 想到分讨,但是暴力输出一下方案之后有很多特别的情况要讨论,就弃了... 假设\(a\)是原序列,\(b\)是我们得到的序列 设\(i\)是最长公共前缀,\( ...
- 牛客练习赛38 D 出题人的手环
链接 [https://ac.nowcoder.com/acm/contest/358/D] 题意 链接:https://ac.nowcoder.com/acm/contest/358/D 来源:牛客 ...
随机推荐
- Oracle - RMAN备份 之 incarnation的实验和小结
对于RMAN恢复我一直都不是很明白,因为,老是搞不清楚不完全恢复该怎么进行,今天,通过这个实验即是对不完全恢复的实践,也是希望搞清楚incarnation到底是怎么工作的.很可惜,本人对Oracle的 ...
- 《windows核心编程系列》十五谈谈windows线程栈
谈谈windows线程栈. 当系统创建线程时会为线程预订一块地址空间区域,注意仅仅是预订.默认情况下预定的这块区域的大小是1MB,虽然预订这么多,但是系统并不会给全部区域调拨物理存储器.默认情况下,仅 ...
- [SDOI2016]墙上的句子
题目描述 考古学家发现了一堵写有未知语言的白色墙壁,上面有一个n行m列的格子,其中有些格子内被填入了某个A至Z的大写字母,还有些格子是空白的. 一直横着或竖着的连续若干个字母会形成一个单词,且每一行的 ...
- 数论 HDOJ 5407 CRB and Candies
题目传送门 题意:求LCM (C(N,0),C(N,1),...,C(N,N)),LCM是最小公倍数的意思,C函数是组合数. 分析:先上出题人的解题报告 好吧,数论一点都不懂,只明白f (n + 1) ...
- CoreText的使用方法
- (void)draw { CGContextRef context = UIGraphicsGetCurrentContext(); NSMutableAttributedString *attr ...
- E. Comments dfs模拟
http://codeforces.com/contest/747/problem/E 首先,把字符串变成这个样子. hello,2,ok,0,bye,0,test,0,one,1,two,2,a,0 ...
- rhel7安装oracle 11gR2,所需的依赖包
binutils-2.23.52.0.1-30.el7.x86_64 compat-libstdc++-33-3.2.3-61.x86_64compat-libstdc++-33-3.2.3-61.i ...
- C++中的数学函数汇总
math.h 数学函数库,一些数学计算的公式的具体实现是放在math.h里,具体有: 1 三角函数 double sin (double); double cos (double); double t ...
- Oracle中的表空间
表空间是什么? Oracle数据库包含逻辑结构和物理结构. 数据库的物理结构是指构成数据库的一组操作系统文件. 数据库的逻辑结构是指描述数据组织方式的一组逻辑概念及它们之间的关系. 表空间是数据库数据 ...
- hihocoder offer收割编程练习赛12 B 一面砖墙
思路: 就是求哪个长度出现的次数最多. 实现: #include <iostream> #include <cstdio> #include <algorithm> ...