标签：路径 维护 位置 有用 排列 oid 证明 快速 移动
给出个排列\(a\)，并且给出\(m\)个区间\([L_i,R_i]\)。
定义一个排列\(b\)对于某区间\([L,R]\)合法，当且仅当对于所有\(i,j\in [L,R],a_i<a_j\)，有\(b_i<b_j\)。对于某个区间集合合法，当且仅当对于每个区间都合法。
前\(k\)个区间的答案为：构造一个有向图\(G\)，使得如果排列\(b\)对于这\(k\)个区间形成的集合合法，当且仅当对于所有\((u,v)\in E_G\)，有\(b_u<b_v\)。问最小的\(|E_G|\)。
对\(k\in [1,m]\)分别求答案。
\(n\le 2.5*10^4,q\le 10^5\)
普通的想法：先对每个区间中，如果有\(a_u<a_v\)就连边\((u,v)\)。构出图之后，删掉一些边，使得仍然满足大小关系。如果\((u,v)\)保留，当且仅当不存在\(x\neq u,v\)，存在路径\(u\to x\to v\)。
结合区间的性质思考问题：假设\(u<v\)。如果\((u,v)\)存在，那么存在同时包含\(u,v\)的区间。在这个条件下，如果\((u,v)\)被删去：如果\([u,v]\)内存在点\(t\)满足\(a_u<a_t<a_v\)，显然要删去；否则，存在一条路径\(u\to p_1\to p_2\to \dots p_k\to v\)，必定有\(p_k<u\)，所以\(p_k,u,v\)在一个区间内，于是就只需要考虑\(k=1\)的情况。也就是说，找到区间\(i\)满足\([u,v]\subseteq [L_i,R_i]\)，且\(L_i\)最小，在\([L_i,v]\)中找到\(t\)满足\(a_t>a_u\)，最小化\(a_t\)，然后判断是否有\(a_u<a_t<a_v\)。
进一步地可以推断出：从某个点\(u\)发出的有向边至多两条，一条向左一条向右。以向右的为例，找到包含\(u\)的区间最大右端点\(maxr(u)\)，在\([u,maxr_u]\)中找到\(t\)满足\(a_t>a_u\)，\(a_t\)最小（记作\(suc_r(u)\)）。向左类似地定义\(suc_l(u)\)。\((u,suc_l(u)),(u,suc_r(u))\)都是边的候选。
一种快速的做法：类似地，换成\(a_t<a_u\)且\(a_t\)最大定义\(pre_r(u),pre_l(u)\)。有结论：假设\(u<v\)，\((u,v)\)被保留当且仅当\(suc_r(u)=v \and pre_l(v)=u\)。证明：考虑上面描述的\((u,v)\)被删去的情况，可以发现在那情况中，\(t\)可以替换\(pre_l(v)\)，所以\(pre_l(v)\neq u\)。于是就可以搞出\(O(n^2)\)的简单做法。由于\(n\)小并且开7s，这种方法可以随便过。
题解做法：
以右边为例，对于每个\(u\)，找到所有可能成为\(suc_r(u)\)的点，即每个\(v\)是区间\([u,v]\)满足\(a_v>a_u\)的\(a_v\)最小的\(v\)。显然\(v_i<v_{i+1},a_v>a_{i+1}\)。
显然，随着时间增加，\(suc_r(u)\)右移。
考虑某个\((u,v_i)\)的存在时间，由以下三个参数决定：
其存在时间为\([t_1,\min(t_2,t_3))\)。
\(u,v_i\)第一次被同时包含的时间可以搞个二维偏序求出。至于\(t_2\)，分别求出左右的\(v_i\)序列，用类似归并的方式搞。对于右边的某个\(v_i\)来说，可能的\(t\)是左边序列（位置递减）中一段后缀，只需要取后缀的第一个作为\(t\)。于是\(v_i\)对应的\(t\)是一定的，接下来要找到最小的\(j\)满足\(L_j\le t,v_i\le R_j\)，也可以二维偏序解决之。
有用的\((u,v)\)只有\(O(n\sqrt q)\)个。证明考虑对于某个区间\([l,r]\)，如果\(u,v\)有连边，则\(a_u,a_v\)在区间内的权值离散化后连续。把它变成\([l\pm 1,r \pm 1]\)，发现总量的变化量是\(O(1)\)的。假设跑个莫队，那么总数不会超过\(O(n\sqrt q)\)。
如何找到有用的\((u,v)\)：直接跑莫队，用set维护相对顺序，可以做到\(O(n\sqrt q\lg n)\)；也可以回滚莫队，把当前块以及之后的点排序，用链表维护，然后倒着删当前块后面的点（得知每个点加入时的前驱后继）。一开始左端点固定在块头（平常是块尾），右端点右移，在链表中加点，对于每个询问移动左端点过来（每次删掉左端点），询问后再把左端点移动回去。那么这样做除了排序之外是\(O(n\sqrt q)\)。
于是总时间\(O(n\sqrt q\lg n)\)。
using namespace std;
#include <bits/stdc++.h>
#define N 25005
#define M 100005
int n,m;
int a[N];
int lmin[N],rmax[N];
int suc[N][2],pre[N][2];
int ans;
void work(int l,int r){
for (int i=l;i<=r;++i){
for (int &j=rmax[i];j<=r;++j)
if (a[j]>a[i]){
int &t=suc[i][0];
if (a[j]<a[t])
ans+=(pre[j][1]==i)-(pre[t][1]==i),t=j;
}
else{
int &t=pre[i][0];
if (a[j]>a[t])
ans+=(suc[j][1]==i)-(suc[t][1]==i),t=j;
}
for (int &j=lmin[i];j>=l;--j)
if (a[j]>a[i]){
int &t=suc[i][1];
if (a[j]<a[t])
ans+=(pre[j][0]==i)-(pre[t][0]==i),t=j;
}
else{
int &t=pre[i][1];
if (a[j]>a[t])
ans+=(suc[j][0]==i)-(suc[t][0]==i),t=j;
}
}
}
int main(){
scanf("%d%d",&n,&m);
for (int i=1;i<=n;++i)
scanf("%d",&a[i]);
a[0]=n+1,a[n+1]=0;
for (int i=1;i<=n;++i){
lmin[i]=i-1,rmax[i]=i+1;
suc[i][0]=suc[i][0]=0;
pre[i][0]=pre[i][1]=n+1;
}
for (int i=1;i<=m;++i){
int l,r;
scanf("%d%d",&l,&r);
work(l,r);
printf("%d\n",ans);
}
return 0;
}
标签：路径 维护 位置 有用 排列 oid 证明 快速 移动
原文地址：https://www.cnblogs.com/jz-597/p/14695626.html
