#include <cstdlib>
#include <cmath>
#include <climits>
#include <cfloat>
#include <map>
#include <set>
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <sstream>
#include <complex>
#include <stack>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iterator>
#include <bitset>
#include <unordered_set>
#include <unordered_map>
#include <fstream>
#include <iomanip>
#include <cassert>
#include <utility>
#include <memory>
#include <functional>
#include <deque>
#include <cctype>
#include <ctime>
#include <numeric>
#include <list>
#include <iomanip>

#if __cplusplus >= 201103L
#include <array>
#include <tuple>
#include <initializer_list>
#include <forward_list>

#define cauto const auto&
#else

#endif

using namespace std;


typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef vector<int> vint;
typedef vector<vector<int> > vvint;
typedef vector<long long> vll;
typedef vector<vector<long long> > vvll;

#define VV(T) vector<vector< T > >

template <class T>
void initvv(vector<vector<T> > &v, int a, int b, const T &t = T()){
    v.assign(a, vector<T>(b, t));
}

template <class F, class T>
void convert(const F &f, T &t){
    stringstream ss;
    ss << f;
    ss >> t;
}

#define GET_MACRO(_1, _2, _3, NAME, ...) NAME
#define _rep(i,n) _rep2((i),0,(n))
#define _rep2(i,a,b) for(int i=(a);i<(b);++i)
#define rep(...) GET_MACRO(__VA_ARGS__, _rep2, _rep)(__VA_ARGS__)
#define ALL(v) (v).begin(),(v).end()
#define PB push_back
#define fi first
#define se second
#define mkp make_pair
#define DEBUG
#ifdef DEBUG
#define dump(x)  cout << #x << " = " << (x) << endl;
#define debug(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl;
#else
#define dump(x) 
#define debug(x) 
#endif

#define MOD 1000000007LL
#define EPS 1e-8
#define INF 0x3f3f3f3f
#define INFL 0x3f3f3f3f3f3f3f3fLL
#define maxs(x,y) x=max(x,y)
#define mins(x,y) x=min(x,y)

template<class T>
class RMQ {
public:
	vector<T> data;
	vector<vector<int> > block;
	vector<int> sblock;
	int N;
	// int lg;
	
	// if Range_MAximam_Query -> return data[b] > data[a] ? b : a
	int argMin(int a, int b) {
		return data[b] < data[a] ? b : a;
	}
	// x の下位 y bit を 0 にする
	int maskBits(int x, int y) {
		return (x >> y) << y;
	}
	RMQ(vector<T> &v) {
		data = v;
		N = data.size();

		// lg = 32 - __builtin_clz(N);

		int blockSize = ((N - 1) >> 5) + 1;
		block.assign(blockSize, vector<int>(32, 0));
		for(int i = 0; i < N; i++) {
			if(i % 32 == 0) block[i >> 5][0] = i;
			block[i >> 5][0] = argMin(block[i >> 5][0], i);
		}

		for(int j = 1; j < block[0].size(); j++) {
			for(int i = 0; i < block.size(); i++) {
				block[i][j] = argMin(block[i][j - 1], block[min(i + (1 << (j - 1)), (int)block.size() - 1)][j - 1]);
			}
		}

		sblock.assign(N, 0);
		vector<int> st(32);
		int stackSize = 0;
		for(int i = 0; i < N; i++) {
			if(i % 32 == 0) stackSize = 0;
			while(stackSize && i == argMin(i, st[stackSize - 1])) --stackSize;
			if(stackSize) {
				int g = st[stackSize - 1];
				sblock[i] = sblock[g] | (1 << (g % 32));
			}
			st[stackSize++] = i;
		}
	}

	// min{data[i] : i in [l .. r]}
	T query(int l, int r) {
		if(l == r) return data[l];
		int l1 = (l >> 5) + 1;
		int r1 = (r >> 5) - 1;
		int ret = l;

		if(l1 <= r1) { // calc sparse table
			int d = r1 - l1 + 1;
			if(d <= 2) {
				ret = argMin(ret, argMin(block[l1][0], block[r1][0]));
			}
			else {
				int lg2 = 32 - __builtin_clz(d) - 1;
				ret = argMin(ret, argMin(block[l1][lg2], block[r1 - (1 << lg2) + 1][lg2]));
			}
		}

		if(l1 - r1 == 2) { // same block
			int t = maskBits(sblock[r], l % 32);
			if(t == 0) {
				ret = argMin(ret, r);
			}
			else {
				ret = argMin(ret, __builtin_ctz(t) + ((l1 - 1) << 5));
			}
		}
		else { // other block
			int t1 = maskBits(sblock[(l1 << 5) - 1], l % 32);
			if(t1 == 0) {
				ret = argMin(ret, (l1 << 5) - 1);
			}
			else {
				ret = argMin(ret, __builtin_ctz(t1) + ((l1 - 1) << 5));
			}

			int t2 = sblock[r];
			if(t2 == 0) {
				ret = argMin(ret, r);
			}
			else {
				ret = argMin(ret, __builtin_ctz(t2) + ((r1 + 1) << 5));
			}
		}

		return data[ret];
	}
};

void mainmain(){
	int n;
	cin>>n;
	vint v(n);
	rep(i,n) cin>>v[i], v[i]--;
	vint rank(n);
	rep(i,n) rank[v[i]] = i;
	vint odd(n, INF);
	vint even(n, INF);
	int mn = INF;
	rep(i,n){
		if(i&1) odd[i] = v[i];
		else even[i] = v[i],mins(mn, v[i]);
	}
	RMQ<int> oRMQ(odd);
	RMQ<int> eRMQ(even);
	using P = pair<int,pii>;
	priority_queue<P> pq;
	pq.push(mkp(-mn, pii(0,n-1)));
	vint ans;
	while(pq.size()){
		int val1;
		pii seg;
		tie(val1, seg) = pq.top();
		pq.pop();
		val1 *= -1;
		// cout<<val1<<" "<<seg.fi<<" "<<seg.se<<endl;
		int pos1 = rank[val1];
		int val2;
		if(seg.fi&1){
			val2 = eRMQ.query(pos1, seg.se);
		}
		else{
			val2 = oRMQ.query(pos1, seg.se);
		}
		int pos2 = rank[val2];
		// cout<<"val "<<val1<<" "<<val2<<endl;
		// cout<<"pos "<<pos1<<" "<<pos2<<endl;
		ans.PB(val1);
		ans.PB(val2);
		if(seg.fi != pos1){
			P tmp;
			tmp.se = pii(seg.fi, pos1-1);
			tmp.fi = -(tmp.se.fi%2==0?eRMQ.query(tmp.se.fi, tmp.se.se): oRMQ.query(tmp.se.fi, tmp.se.se));
			pq.push(tmp);
		}
		if(pos1+1 != pos2){
			P tmp;
			tmp.se = pii(pos1+1, pos2-1);
			tmp.fi = -(tmp.se.fi%2==0?eRMQ.query(tmp.se.fi, tmp.se.se): oRMQ.query(tmp.se.fi, tmp.se.se));
			pq.push(tmp);	
		}
		if(pos2 != seg.se){
			P tmp;
			tmp.se = pii(pos2+1, seg.se);
			tmp.fi = -(tmp.se.fi%2==0?eRMQ.query(tmp.se.fi, tmp.se.se): oRMQ.query(tmp.se.fi, tmp.se.se));
			pq.push(tmp);
		}
	}
	rep(i,n){
		if(i) cout<<" ";
		cout<<ans[i]+1;
	}
	cout<<endl;
}


signed main() {
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout<<fixed<<setprecision(20);
    mainmain();
}

Measured memory_bytes: 11296