Codeforces 1696D Permutation Graph

转换题意：假设p在l位置，如果[l,r]是一个最值都在端点处的区间，那么p可以移动到r，问最小移动次数是多少？

考虑一个样例：2 1 5 3 4

对于每次走，我们应该是贪心的选择最大区间来走，那对每个序列，假设pos1,pos2是两个最值的位置，并且pos1<pos2，那么[pos1,pos2]是最大区间，为什么呢？

对pos1的左边的元素，如果他是区间最值，那么pos1就不会是区间最值，而pos1位置的元素是整个序列的最值，互相矛盾，所有每次只需取[pos1,pos2]作为一次操作。此时只需要找到[l,pos1]和[pos2,r]的操作最小值即可。

对于不修改地查询序列最值，ST表效率比线段树高些，每次查询能做到的复杂度。

#include<iostream>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<ctype.h>
#include<algorithm>
#include<string>
#include<vector>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<deque>
#define ll long long
#define ull unsigned long long
#define fr first
#define se second
#define endl '\n'
using namespace std;

const int inf=0x3fffffff;
const long long INF=0x3f3f3f3f3f3f;
const double eps=1e-8;
const double pi=2*acos(0.0);

template<typename T>
T father(T x){return x&1?x-1>>1:x>>1;}
template<typename T>
T ls(T x){return x<<1;}
template<typename T>
T rs(T x){return x<<1|1;}
template<typename T>
T abs_(T a){return a<0?-a:a;}
template<typename T>
T ceil_(T a,T b){return a%b?a/b+1:a/b;}
template<typename T>
T log_(T a,T b){return log(b)/log(a);}
template<typename T>
T gcd(T a,T b){return !b?a:gcd(b,a%b);}
template<typename T>
T lcm(T a,T b){return a/gcd(a,b)*b;}
template<typename T>
void exgcd(T a,T b,T& x,T& y)
{
	if(!b){x=1;y=0;return;}
	exgcd(b,a%b,x,y);
	T temp=x;x=y;
	y=temp-a/b*y;
}
template<typename T>
T solve_function(T a,T b,T c)
{
	if(a<0){a=-a;c=-c;}
	T x,y,base=gcd(a,b);
	exgcd(a,b,x,y);
	x*=c/base;T p=b/base;
	x=(x%p+p)%p;
	if(x<0) x+=abs_(p);
	return x;
}
template<typename T>
T reverse_element(T a,T k){return solve_function(a,k,(T)1);}
template<typename T>
T slow_multi(T a,T b,T p)
{
	T ret=0;
	while(b)
	{
		if(b&1) ret=(ret+a)%p;
		a=(a<<1)%p;
		b>>=1;
	}
	return ret;
}
template<typename T>
T qpow(T a,T b)
{
	T ret=1,base=a;
	while(b)
	{
		if(b&1) ret=ret*base;
		base=base*base;
		b>>=1;
	}
	return ret;
}

int read()
{
	int ret=0,base=1;
	char ch=getchar();
	while(!isdigit(ch))
	{
		if(ch=='-') base=-1;
		ch=getchar();
	}
	while(isdigit(ch))
	{
		ret=(ret<<1)+(ret<<3)+ch-48;
		ch=getchar();
	}
	return ret*base;
}

int n,pos[500005],a[500005],st1[500005][23],st2[500005][23];

int query_min(int l,int r)
{
	int len=log_(2,r-l+1);
	return min(st1[l][len],st1[r-(1<<len)+1][len]);
}
 
int query_max(int l,int r)
{
	int len=log_(2,r-l+1);
	return max(st2[l][len],st2[r-(1<<len)+1][len]);
}

int make(int l,int r)
{
	if(l>=r) return 0;
	int pos1=pos[query_min(l,r)],pos2=pos[query_max(l,r)];
	if(pos1>pos2) swap(pos1,pos2);
	return 1+make(l,pos1)+make(pos2,r);
}

void work()
{
	n=read();
	for(int i=1;i<=n;i++) pos[st1[i][0]=st2[i][0]=a[i]=read()]=i;
	for(int j=1;j<=22;j++)
	{
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			st1[i][j]=min(st1[i][j-1],st1[i+(1<<j-1)][j-1]);
			st2[i][j]=max(st2[i][j-1],st2[i+(1<<j-1)][j-1]);
		}
	}
	printf("%d\n",make(1,n));
}

int main()
{
	int t=read();
	while(t--) work();
 	return 0;
}