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;
}
