//#pragma GCC optimize("Ofast,unroll-loops")
#include <bits/stdc++.h>
using namespace std;
typedef long long llo;
#define mp make_pair
#define pb push_back
#define a first
#define b second
#define endl '\n'
#include "books.h"
llo ans=0;
llo n;
vector<llo> adj[1000001];
llo vis[1000001];
pair<llo,llo> cot[1000001];
llo par[1000001];
llo find(llo no){
if(par[no]==no){
return no;
}
par[no]=find(par[no]);
return par[no];
}
vector<llo> ss;
void dfs(int no){
vis[no]=1;
ss.pb(no);
for(auto j:adj[no]){
if(vis[j]==0){
dfs(j);
}
}
}
vector<llo> adj2[1000001];
llo ind[1000001];
pair<llo,llo> pp;
vector<pair<llo,llo>> tt;
llo dp[1001][1001];
pair<llo,llo> rr[1000001];
void dfs2(int no,int par=-1){
pp.a=min(pp.a,tt[no].a);
pp.b=max(pp.b,tt[no].b);
vis[no]=1;
for(auto j:adj2[no]){
if(vis[j]==0){
dfs2(j);
}
}
}
long long minimum_walk(std::vector<int> p, int s) {
while(p.size()){
if(p.back()+1==p.size()){
p.pop_back();
}
else{
break;
}
}
n=p.size();
//cout<<p.size()<<":"<<endl;
if(n==0){
return 0;
}
for(int i=0;i<n;i++){
adj[i].pb(p[i]);
}
//llo co=0;
llo ans=0;
for(int i=0;i<n;i++){
if(vis[i]==0){
ss.clear();
dfs(i);
/*if(ss.size()==0){
continue;
}*/
llo mi=ss[0];
llo ma=ss[0];
for(auto j:ss){
mi=min(mi,j);
ma=max(ma,j);
}
for(auto j:ss){
cot[j]={mi,ma};
}
for(auto j:ss){
ind[j]=tt.size();
}
tt.pb({mi,ma});
//co++;
for(int j=0;j+1<ss.size();j++){
ans+=abs(ss[j]-ss[j+1]);
}
/* for(auto j:ss){
cout<<j<<",,";
}
cout<<endl;*/
ans+=(abs(ss[0]-ss.back()));
}
}
/* if(s==0){
sort(tt.begin(),tt.end());
vector<pair<llo,llo>> tt2;
for(auto j:tt){
if(tt2.size()==0){
tt2.pb(j);
}
else{
if(tt2.back().b<j.a){
tt2.pb(j);
}
else{
tt2[tt2.size()-1].b=max(tt2[tt2.size()-1].b,j.b);
}
}
}
llo xx=tt2.size();
ans+=2*(xx-1);
//ans+=2*(co-1);
return ans;
}*/
for(int i=0;i<tt.size();i++){
for(int j=0;j<n;j++){
if(ind[j]!=i){
if(j>=tt[i].a and j<=tt[i].b){
adj2[i].pb(ind[j]);
//cout<<i<<",,"<<ind[j]<<endl;
}
}
}
}
for(int j=0;j<tt.size();j++){
for(int i=0;i<tt.size();i++){
vis[i]=0;
}
pp=tt[j];
dfs2(j);
for(int i=0;i<n;i++){
if(ind[i]==j){
rr[i]=pp;
}
}
//cout<<pp.a<<","<<pp.b<<endl;
}
pair<llo,llo> cur3={rr[s].a,rr[s].b};
while(cur3.a>0 or cur3.b<n-1){
if(cur3.b==n-1){
ans+=2;
cur3={rr[cur3.a-1].a,max(cur3.b,rr[cur3.a-1].b)};
continue;
}
else if(cur3.a==0){
ans+=2;
cur3={min(rr[cur3.b+1].a,cur3.a),rr[cur3.b+1].b};
continue;
}
else{
ans+=2;
if(rr[cur3.a-1].b>cur3.b){
cur3={rr[cur3.a-1].a,max(cur3.b,rr[cur3.a-1].b)};
continue;
}
else if(rr[cur3.b+1].a<cur3.a){
cur3={min(rr[cur3.b+1].a,cur3.a),rr[cur3.b+1].b};
continue;
}
cur3={rr[cur3.a-1].a,max(cur3.b,rr[cur3.a-1].b)};
}
}
return ans;
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
dp[i][j]=-1;
}
}
dp[rr[s].a][rr[s].b]=0;
for(llo ii=0;ii<n;ii++){
for(llo i=0;i+ii<n;i++){
llo j=i+ii;
if(dp[i][j]==-1){
continue;
}
if(i>0){
pair<llo,llo> ac={rr[i-1].a,max(rr[i-1].b,j)};
if(dp[ac.a][ac.b]==-1){
dp[ac.a][ac.b]=dp[i][j]+2;
}
dp[ac.a][ac.b]=min(dp[ac.a][ac.b],dp[i][j]+2);
}
if(j<n-1){
pair<llo,llo> ac={min(rr[j+1].a,i),rr[j+1].b};
if(dp[ac.a][ac.b]==-1){
dp[ac.a][ac.b]=dp[i][j]+2;
}
dp[ac.a][ac.b]=min(dp[ac.a][ac.b],dp[i][j]+2);
}
}
}
for(int i=0;i<n;i++){
for(int j=i+1;j<n;j++){
if(rr[i].b<rr[j].a){
continue;
}
if(rr[j].b<rr[i].a){
continue;
}
if(rr[j].a<=rr[i].a and rr[i].b<=rr[j].b){
continue;
}
if(rr[i].a<=rr[j].a and rr[j].b<=rr[i].b){
continue;
}
while(true){
continue;
}
}
}
/*for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
cout<<dp[i][j]<<",";
}
cout<<endl;
}
*/
//ans+=2*(co-1);
//cout<<ans<<":"<<endl;
return ans+dp[0][n-1];
}
