Problem: You are given a set of n segments on a line [L_i; R_i]. All 2n segment endpoints are pairwise distinct integers.

The set is laminar — any two segments are either disjoint or one of them contains the other.

Choose a non-empty subsegment [l_i, r_i] with integer endpoints in each segment (L_i ≤ l_i < r_i ≤ R_i) in such a way that no two subsegments intersect (they are allowed to have common endpoints though) and the sum of their lengths (∑_{i=1}^n r_i - l_i) is maximized.

Input

The first line contains a single integer n (1 ≤ n ≤ 2 ⋅ 10^3) — the number of segments.

The i-th of the next n lines contains two integers L_i and R_i (0 ≤ L_i < R_i ≤ 10^9) — the endpoints of the i-th segment.

All the given 2n segment endpoints are distinct. The set of segments is laminar.

Output

On the first line, output the maximum possible sum of subsegment lengths.

On the i-th of the next n lines, output two integers l_i and r_i (L_i ≤ l_i < r_i ≤ R_i), denoting the chosen subsegment of the i-th segment.

Example

Input


4
1 10
2 3
5 9
6 7


Output


7
3 6
2 3
7 9
6 7

Note

The example input and the example output are illustrated below. <image> <image>

Solution (in CPP):

Solution:
//Original Code:
//#include <self/utility>
//#include <self/debug>
//using namespace std;
//int n;
//int l[2005];
//int r[2005];
//vector<int> edges[2005];
//vector<int> roots;
//int resl[2005];
//int resr[2005];
//int szlim[2005];
//
//inline void compressAddress(ll &t,int a,int b,int c,int d,int e,int ch=0)
//{
//	t=0;
//	(((((((((t|=a)<<=11)|=b)<<=11)|=c)<<=11)|=d)<<=11)|=e);
//	if(ch) t|=(1LL<<60);
//}
//
//inline int decompressAddress(ll t,int &a,int &b,int &c,int &d,int &e)
//{
//	int ch=(bool)(t&(1LL<<60));
//	int mask=(1<<11)-1;
//	e=t&mask;t>>=11;
//	d=t&mask;t>>=11;
//	c=t&mask;t>>=11;
//	b=t&mask;t>>=11;
//	a=t&mask;
//	return ch;
//}
//
//struct Value
//{
//	int a[2][2];
//	pair<ll,ll> trace[2][2];
//	int* operator [] (int x){
//		return a[x];
//	}
//	const int* operator [] (int x)const{
//		return a[x];
//	}
//}vmin;
//vector<vector<Value> > dp[2005];
//int sz[2005];
//int depth[2005];
//int fn[2005];//???
///**
// * ?????????:
// * ????????
// * ?(?????????)+(????????)+(?????????)
// */
//
//
//int getInfo(int x,int d=0)
//{
//	depth[x]=d;
//	sz[x]=1;
//	for(auto u:edges[x])
//	{
//		sz[x]+=getInfo(u,d+1);
//	}
//	fn[x]=(sz[x]==1?0:edges[x].size()+1);
//	szlim[x]=min(sz[x]-1,depth[x])+2;
//	return sz[x];
//}
//
//void merge(vector<Value> &a,vector<Value> &b,vector<Value>&c,int pm,int x,int u,int id)
//{
//	int n=a.size()-1,m=b.size()-1,t=n+m+1;
//	c.resize(min(szlim[x],t+1),vmin);
//	for(int i=n;i>=0;i--)
//	{
//		for(int j=0;j<=m;j++)
//		{
//			for(int l1=0;l1<2;l1++)
//			{
//				for(int l2=0;l2<2;l2++)
//				{
//					for(int r1=0;r1<2;r1++)
//					{
//						for(int r2=0;r2<2;r2++)
//						{
//							int forceCh=false;
//							if(l2 || r1) forceCh=true;
//							for(int ch=forceCh;ch<2;ch++)
//							{
//								if(ch) assert(pm!=0);
//								if(i+j+ch>=szlim[x]) continue;
//								// chmax(c[i+j+ch][l1][r2],a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0));
//								if(c[i+j+ch][l1][r2]<a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0))
//								{
//									// if(x==0 && id==1 && i+j+ch==1 && l1==0 && r2==0)
//									// {
//									// 	debug(ch);
//									// }
//									c[i+j+ch][l1][r2]=a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0);
//									compressAddress(c[i+j+ch].trace[l1][r2].first,x,id,i,l1,l2,ch);
//									// int ta,tb,tc,td,te;
//									// ch=decompressAddress(c[i+j+ch].trace[l1][r2].first,ta,tb,tc,td,te);
//									compressAddress(c[i+j+ch].trace[l1][r2].second,u,fn[u],j,r1,r2);
//								}
//							}
//						}
//					}
//				}
//			}
//		}
//	}
//}
//
//void dfs(int x)
//{
//	for(auto u:edges[x])
//	{
//		dfs(u);
//	}
//	int n=edges[x].size();
//	int m=1;
//	dp[x].resize(fn[x]+1);
//	dp[x][0].resize(m+1,vmin);
//	if(n==0)//??????????????
//	{
//		dp[x][0][0][0][0]=r[x]-l[x];
//		return;
//	}
//	else//???,??????0??????,???????????????
//	{
//		int u=edges[x][0],lef=l[u]-l[x];
//		m+=sz[u];
//		dp[x][1]=dp[u][fn[u]];
//		dp[x][1].resize(m+1,vmin);
//		for(int i=0;i<dp[u][fn[u]].size();i++)
//		{
//			for(int l=0;l<2;l++)
//			{
//				for(int r=0;r<2;r++)
//				{
//					compressAddress(dp[x][1][i].trace[l][r].first,u,fn[u],i,l,r);
//					dp[x][1][i].trace[l][r].second=-llinf-1;
//				}
//			}
//		}
//		//?????????x???????????
//		for(int i=m;i>=0;i--)
//		{
//			for(int r=0;r<2;r++)
//			{
//				// chmax(dp[x][1][i][1][r],dp[x][1][i][0][r]);
//				if(dp[x][1][i][1][r]<dp[x][1][i][0][r])
//				{
//					dp[x][1][i][1][r]=dp[x][1][i][0][r];
//					compressAddress(dp[x][1][i].trace[1][r].first,x,1,i,0,r);
//					dp[x][1][i].trace[1][r].second=-1-l[u];
//				}
//			}
//		}
//		//????????? 
//		for(int i=0;i<=m;i++)
//		{
//			for(int r=0;r<2;r++)
//			{
//				dp[x][1][i][1][r]+=lef;
//			}
//		}
//	}
//
//	for(int i=1;i<n;i++)
//	{
//		int u=edges[x][i],md=l[edges[x][i]]-r[edges[x][i-1]];
//		merge(dp[x][i],dp[u][fn[u]],dp[x][i+1],md,x,u,i);
//	}
//
//	//?????????,??????????????
//	vector<Value> &v=dp[x][n];int rgt=r[x]-r[edges[x][n-1]];
//	dp[x][n].resize(sz[x]+1,vmin);
//	for(int i=sz[x];i>=0;i--)
//	{
//		for(int l=0;l<2;l++)
//		{
//			// chmax(v[i][l][1],v[i][l][0]);
//			if(v[i][l][1]<v[i][l][0])
//			{
//				v[i][l][1]=v[i][l][0];
//				compressAddress(v[i].trace[l][1].first,x,n,i,l,0);
//				v[i].trace[l][1].second=-1-r[edges[x][n-1]];
//			}
//		}
//	}
//	for(int i=0;i<=sz[x];i++)
//	{
//		for(int l=0;l<2;l++)
//		{
//			v[i][l][1]+=rgt;
//		}
//	}
//	//????,??????/?????,???????????
//	//??????????????,????????-1
//	//?????????,?????????
//	dp[x][fn[x]].resize(min(szlim[x],sz[x]),vmin);
//	for(int i=0;i<=sz[x];i++)
//	{
//		for(int l=0;l<2;l++)
//		{
//			for(int r=0;r<2;r++)
//			{
//				for(int fl=l;fl<2;fl++)
//				{
//					for(int fr=r;fr<2;fr++)
//					{
//						if(i+(fl-l)+(fr-r)-1>=min(szlim[x],sz[x])) continue;
//						if(i+(fl-l)+(fr-r)-1<0) continue;
//						if(dp[x][fn[x]][i+(fl-l)+(fr-r)-1][l][r]<v[i][fl][fr])
//						{
//							dp[x][fn[x]][i+(fl-l)+(fr-r)-1][l][r]=v[i][fl][fr];
//							compressAddress(dp[x][fn[x]][i+(fl-l)+(fr-r)-1].trace[l][r].first,x,n,i,fl,fr);
//							dp[x][fn[x]][i+(fl-l)+(fr-r)-1].trace[l][r].second=-llinf;
//						}
//					}
//				}
//			}
//		}
//	}
//}
//
//vector<int> chSegs;
//
//void pushSeg(int a,int b,int c,int d,int e)
//{
//	// debug(a,b,c,d,e,dp[a][b][c][d][e]);
//	if(edges[a].size()==0)//???????????????
//	{
//		assert((b|c|d|e)==0);
//		chSegs.push_back(l[a]);
//		chSegs.push_back(r[a]);
//		return;
//	}
//	int na,nb,nc,nd,ne;
//	pair<ll,ll> &o=dp[a][b][c].trace[d][e];
//	int ch=decompressAddress(o.first,na,nb,nc,nd,ne);
//	// debug(ch);
//	// debug(a,b,c,d,e,na,nb,nc,nd,ne);
//	if(o.second>=0)
//	{
//		// if(a==2 && ch)
//		// {
//		// 	debug(b,ne,r[edges[a][b-2]]);
//		// }
//		if(ch && ne==0)
//		{
//			chSegs.push_back(r[edges[a][b-2]]);
//		}
//		pushSeg(na,nb,nc,nd,ne);
//		decompressAddress(o.second,na,nb,nc,nd,ne);
//		if(ch && nd==0)
//		{
//			chSegs.push_back(l[edges[a][b-1]]);
//		}
//		pushSeg(na,nb,nc,nd,ne);
//	}
//	else if(o.second==-llinf-1)
//	{
//		pushSeg(na,nb,nc,nd,ne);
//	}
//	else
//	{
//		// debug("FUCK");
//		// assert(na==a);
//		// assert(nb==b);
//		// assert(mabs(nc-c)<=1);
//		int pl,pr;
//		// if(o.second==-llinf)
//		// {
//		// 	debug(a,nd!=d,ne!=e);
//		// }
//		if(o.second==-llinf) pl=l[a],pr=r[a];
//		else pl=pr=-o.second-1;
//		if(nd!=d)
//		{
//			chSegs.push_back(pl);
//		}
//		if(ne!=e)
//		{
//			if(o.second==-llinf && a==3) debug(pr);
//			chSegs.push_back(pr);
//		}
//		pushSeg(na,nb,nc,nd,ne);
//	}
//}
//
//int solve(int x)
//{
//	getInfo(x);
//	dfs(x);
//	int val=dp[x][fn[x]][0][0][0];
//	pushSeg(x,fn[x],0,0,0);
//	return val;
//}
//
//int main()
//{
//	// freopen("input.txt","r",stdin);
//	for(int i=0;i<2;i++) for(int j=0;j<2;j++) vmin[i][j]=-inf;
//
//	cin>>n;
//	for(int i=0;i<n;i++)
//	{
//		cin>>l[i]>>r[i];
//	}
//	for(int i=0;i<n;i++)
//	{
//		pii fa=mp(-inf,-1);
//		for(int j=0;j<n;j++)
//		{
//			if(l[j]<l[i] && r[j]>r[i]) chmax(fa,mp(l[j],j));
//		}
//		if(fa.second!=-1) edges[fa.second].push_back(i);else roots.push_back(i);
//	}
//	// debug(roots);
//
//	for(int i=0;i<n;i++)
//	{
//		sort(all(edges[i]),[](const int &a,const int &b){return l[a]<l[b];});
//	}
//
//	vector<pii> allres;
//	int res=0;
//	for(auto x:roots)
//	{
//		res+=solve(x);
//	}
//	cout<<res<<endl;
//	sort(all(chSegs));
//	assert(chSegs.size()==n*2);
//	memset(resl,-1,sizeof(resl));
//	for(int i=0;i<n*2;i+=2)
//	{
//		int lp=chSegs[i];
//		int rp=chSegs[i+1];
//		pii o=mp(-inf,-1);
//		for(int j=0;j<n;j++)
//		{
//			if(resl[j]==-1 && l[j]<=lp && r[j]>=rp) chmax(o,mp(l[j],j));
//		}
//		assert(o.second!=-1);
//		resl[o.second]=lp;
//		resr[o.second]=rp;
//	}
//	for(int i=0;i<n;i++)
//	{
//		assert(resl[i]>=l[i] && resr[i]<=r[i]);
//		res-=resr[i]-resl[i];
//		printf("%d %d\n",resl[i],resr[i]);
//	}
//	assert(res==0);
//	// return res;
//	// debug(res);
//	// debug(chSegs);
//	// for(int i=0;i<n*2;i+=2)
//	// {
//	// 	printf("%d %d\n",chSegs[i],chSegs[i+1]);
//	// }
//	return 0;
//}

//substituted with C++ Inliner
#ifndef _SELF_UTILITY
#define _SELF_UTILITY 1
#include <numeric>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <stdio.h>
#include <stdlib.h>
#include <vector>
#include <map>
#include <queue>
#include <set>
#include <string>
#include <string.h>
#include <stack>
#include <assert.h>
#include <bitset>
#include <time.h>
#define Endl endl
#define mp make_pair
#define mt make_tuple
#define ll long long 
#define ull unsigned long long
#define pii pair<int,int>
#define over(A) {cout<<A<<endl;exit(0);}
#define all(A) A.begin(),A.end()
#define quickcin ios_base::sync_with_stdio(false);
#ifdef __TAKE_MOD
int mod;
#else
#ifdef __TAKE_CONST_MOD
const int mod=__TAKE_CONST_MOD;
#else
const int mod=1000000007;
#endif
#endif
const int gmod=3;
const int inf=1039074182;
const double eps=1e-9;
const ll llinf=2LL*inf*inf;
template <typename T1,typename T2> inline void chmin(T1 &x,T2 b) {if(b<x) x=b;}
template <typename T1,typename T2> inline void chmax(T1 &x,T2 b) {if(b>x) x=b;}
inline void chadd(int &x,int b) {x+=b-mod;x+=(x>>31 & mod);}
template <typename T1,typename T2> inline void chadd(T1 &x,T2 b) {x+=b;if(x>=mod) x-=mod;}
template <typename T1,typename T2> inline void chmul(T1 &x,T2 b) {x=1LL*x*b%mod;}
template <typename T1,typename T2> inline void chmod(T1 &x,T2 b) {x%=b,x+=b;if(x>=b) x-=b;}
template <typename T> inline T mabs(T x) {return (x<0?-x:x);}
using namespace std;
#endif
#ifndef _SELF_DEBUG
#define _SELF_DEBUG 1
#ifndef _SELF_OPERATOR
#define _SELF_OPERATOR 1
using namespace std;
template <typename T>
ostream &operator<<(ostream &cout, vector<T> vec)
{
	cout << "{";
	for (int i = 0; i < (int)vec.size(); i++)
	{
		cout << vec[i];
		if (i != (int)vec.size() - 1)
			cout << ',';
	}
	cout << "}";
	return cout;
}

template <typename T1, typename T2>
ostream &operator<<(ostream &cout, pair<T1, T2> p)
{
	cout << "(" << p.first << ',' << p.second << ")";
	return cout;
}

template <typename T, typename T2>
ostream &operator<<(ostream &cout, set<T, T2> s)
{
	vector<T> t;
	for (auto x : s)
		t.push_back(x);
	cout << t;
	return cout;
}

template <typename T, typename T2>
ostream &operator<<(ostream &cout, multiset<T, T2> s)
{
	vector<T> t;
	for (auto x : s)
		t.push_back(x);
	cout << t;
	return cout;
}

template <typename T>
ostream &operator<<(ostream &cout, queue<T> q)
{
	vector<T> t;
	while (q.size())
	{
		t.push_back(q.front());
		q.pop();
	}
	cout << t;
	return cout;
}

template <typename T1, typename T2, typename T3>
ostream &operator<<(ostream &cout, map<T1, T2, T3> m)
{
	for (auto &x : m)
	{
		cout << "Key: " << x.first << ' ' << "Value: " << x.second << endl;
	}
	return cout;
}

template <typename T1, typename T2>
void operator+=(pair<T1, T2> &x,const pair<T1, T2> y)
{
	x.first += y.first;
	x.second += y.second;
}
template <typename T1,typename T2>
pair<T1,T2> operator + (const pair<T1,T2> &x,const pair<T1,T2> &y)
{
	return make_pair(x.first+y.first,x.second+y.second);
}

template <typename T1,typename T2>
pair<T1,T2> operator - (const pair<T1,T2> &x,const pair<T1,T2> &y)
{
	return mp(x.first-y.first,x.second-y.second);
}

template <typename T1, typename T2>
pair<T1, T2> operator-(pair<T1, T2> x)
{
	return make_pair(-x.first, -x.second);
}

template <typename T>
vector<vector<T>> operator~(vector<vector<T>> vec)
{
	vector<vector<T>> v;
	int n = vec.size(), m = vec[0].size();
	v.resize(m);
	for (int i = 0; i < m; i++)
	{
		v[i].resize(n);
	}
	for (int i = 0; i < m; i++)
	{
		for (int j = 0; j < n; j++)
		{
			v[i][j] = vec[j][i];
		}
	}
	return v;
}
#endif
#include <sstream>
void print0x(int x)
{
	std::vector <int> vec;
	while(x)
	{
		vec.push_back(x&1);
		x>>=1;
	}
	std::reverse(vec.begin(),vec.end());
	for(int i=0;i<(int)vec.size();i++)
	{
		std::cout<<vec[i];
	}
	std::cout<<' ';
}

template <typename T>
void print0x(T x,int len)
{
	std::vector <int> vec;
	while(x)
	{
		vec.push_back(x&1);
		x>>=1;
	}
	reverse(vec.begin(),vec.end());
	for(int i=(int)vec.size();i<len;i++)
	{
		putchar('0');
	}
	for(size_t i=0;i<vec.size();i++)
	{
		std::cout<<vec[i];
	}
	std::cout<<' ';
}
vector<string> vec_splitter(string s) {
	s += ',';
	vector<string> res;
	while(!s.empty()) {
		res.push_back(s.substr(0, s.find(',')));
		s = s.substr(s.find(',') + 1);
	}
	return res;
}
void debug_out(
vector<string> __attribute__ ((unused)) args,
__attribute__ ((unused)) int idx, 
__attribute__ ((unused)) int LINE_NUM) { cerr << endl; } 
template <typename Head, typename... Tail>
void debug_out(vector<string> args, int idx, int LINE_NUM, Head H, Tail... T) {
	if(idx > 0) cerr << ", "; else cerr << "Line(" << LINE_NUM << ") ";
	stringstream ss; ss << H;
	cerr << args[idx] << " = " << ss.str();
	debug_out(args, idx + 1, LINE_NUM, T...);
}
#define debug(...) debug_out(vec_splitter(#__VA_ARGS__), 0, __LINE__, __VA_ARGS__)
#endif
using namespace std;
int n;
int l[2005];
int r[2005];
vector<int> edges[2005];
vector<int> roots;
int resl[2005];
int resr[2005];
int szlim[2005];

inline void compressAddress(ll &t,int a,int b,int c,int d,int e,int ch=0)
{
	t=0;
	(((((((((t|=a)<<=11)|=b)<<=11)|=c)<<=11)|=d)<<=11)|=e);
	if(ch) t|=(1LL<<60);
}

inline int decompressAddress(ll t,int &a,int &b,int &c,int &d,int &e)
{
	int ch=(bool)(t&(1LL<<60));
	int mask=(1<<11)-1;
	e=t&mask;t>>=11;
	d=t&mask;t>>=11;
	c=t&mask;t>>=11;
	b=t&mask;t>>=11;
	a=t&mask;
	return ch;
}

struct Value
{
	int a[2][2];
	pair<ll,ll> trace[2][2];
	int* operator [] (int x){
		return a[x];
	}
	const int* operator [] (int x)const{
		return a[x];
	}
}vmin;
vector<vector<Value> > dp[2005];
int sz[2005];
int depth[2005];
int fn[2005];//???
/**
 * ?????????:
 * ????????
 * ?(?????????)+(????????)+(?????????)
 */


int getInfo(int x,int d=0)
{
	depth[x]=d;
	sz[x]=1;
	for(auto u:edges[x])
	{
		sz[x]+=getInfo(u,d+1);
	}
	fn[x]=(sz[x]==1?0:edges[x].size()+1);
	szlim[x]=min(sz[x]-1,depth[x])+2;
	return sz[x];
}

void merge(vector<Value> &a,vector<Value> &b,vector<Value>&c,int pm,int x,int u,int id)
{
	int n=a.size()-1,m=b.size()-1,t=n+m+1;
	c.resize(min(szlim[x],t+1),vmin);
	for(int i=n;i>=0;i--)
	{
		for(int j=0;j<=m;j++)
		{
			for(int l1=0;l1<2;l1++)
			{
				for(int l2=0;l2<2;l2++)
				{
					for(int r1=0;r1<2;r1++)
					{
						for(int r2=0;r2<2;r2++)
						{
							int forceCh=false;
							if(l2 || r1) forceCh=true;
							for(int ch=forceCh;ch<2;ch++)
							{
								if(ch) assert(pm!=0);
								if(i+j+ch>=szlim[x]) continue;
								// chmax(c[i+j+ch][l1][r2],a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0));
								if(c[i+j+ch][l1][r2]<a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0))
								{
									// if(x==0 && id==1 && i+j+ch==1 && l1==0 && r2==0)
									// {
									// 	debug(ch);
									// }
									c[i+j+ch][l1][r2]=a[i][l1][l2]+b[j][r1][r2]+(ch?pm:0);
									compressAddress(c[i+j+ch].trace[l1][r2].first,x,id,i,l1,l2,ch);
									// int ta,tb,tc,td,te;
									// ch=decompressAddress(c[i+j+ch].trace[l1][r2].first,ta,tb,tc,td,te);
									compressAddress(c[i+j+ch].trace[l1][r2].second,u,fn[u],j,r1,r2);
								}
							}
						}
					}
				}
			}
		}
	}
}

void dfs(int x)
{
	for(auto u:edges[x])
	{
		dfs(u);
	}
	int n=edges[x].size();
	int m=1;
	dp[x].resize(fn[x]+1);
	dp[x][0].resize(m+1,vmin);
	if(n==0)//??????????????
	{
		dp[x][0][0][0][0]=r[x]-l[x];
		return;
	}
	else//???,??????0??????,???????????????
	{
		int u=edges[x][0],lef=l[u]-l[x];
		m+=sz[u];
		dp[x][1]=dp[u][fn[u]];
		dp[x][1].resize(m+1,vmin);
		for(int i=0;i<dp[u][fn[u]].size();i++)
		{
			for(int l=0;l<2;l++)
			{
				for(int r=0;r<2;r++)
				{
					compressAddress(dp[x][1][i].trace[l][r].first,u,fn[u],i,l,r);
					dp[x][1][i].trace[l][r].second=-llinf-1;
				}
			}
		}
		//?????????x???????????
		for(int i=m;i>=0;i--)
		{
			for(int r=0;r<2;r++)
			{
				// chmax(dp[x][1][i][1][r],dp[x][1][i][0][r]);
				if(dp[x][1][i][1][r]<dp[x][1][i][0][r])
				{
					dp[x][1][i][1][r]=dp[x][1][i][0][r];
					compressAddress(dp[x][1][i].trace[1][r].first,x,1,i,0,r);
					dp[x][1][i].trace[1][r].second=-1-l[u];
				}
			}
		}
		//????????? 
		for(int i=0;i<=m;i++)
		{
			for(int r=0;r<2;r++)
			{
				dp[x][1][i][1][r]+=lef;
			}
		}
	}

	for(int i=1;i<n;i++)
	{
		int u=edges[x][i],md=l[edges[x][i]]-r[edges[x][i-1]];
		merge(dp[x][i],dp[u][fn[u]],dp[x][i+1],md,x,u,i);
	}

	//?????????,??????????????
	vector<Value> &v=dp[x][n];int rgt=r[x]-r[edges[x][n-1]];
	dp[x][n].resize(sz[x]+1,vmin);
	for(int i=sz[x];i>=0;i--)
	{
		for(int l=0;l<2;l++)
		{
			// chmax(v[i][l][1],v[i][l][0]);
			if(v[i][l][1]<v[i][l][0])
			{
				v[i][l][1]=v[i][l][0];
				compressAddress(v[i].trace[l][1].first,x,n,i,l,0);
				v[i].trace[l][1].second=-1-r[edges[x][n-1]];
			}
		}
	}
	for(int i=0;i<=sz[x];i++)
	{
		for(int l=0;l<2;l++)
		{
			v[i][l][1]+=rgt;
		}
	}
	//????,??????/?????,???????????
	//??????????????,????????-1
	//?????????,?????????
	dp[x][fn[x]].resize(min(szlim[x],sz[x]),vmin);
	for(int i=0;i<=sz[x];i++)
	{
		for(int l=0;l<2;l++)
		{
			for(int r=0;r<2;r++)
			{
				for(int fl=l;fl<2;fl++)
				{
					for(int fr=r;fr<2;fr++)
					{
						if(i+(fl-l)+(fr-r)-1>=min(szlim[x],sz[x])) continue;
						if(i+(fl-l)+(fr-r)-1<0) continue;
						if(dp[x][fn[x]][i+(fl-l)+(fr-r)-1][l][r]<v[i][fl][fr])
						{
							dp[x][fn[x]][i+(fl-l)+(fr-r)-1][l][r]=v[i][fl][fr];
							compressAddress(dp[x][fn[x]][i+(fl-l)+(fr-r)-1].trace[l][r].first,x,n,i,fl,fr);
							dp[x][fn[x]][i+(fl-l)+(fr-r)-1].trace[l][r].second=-llinf;
						}
					}
				}
			}
		}
	}
}

vector<int> chSegs;

void pushSeg(int a,int b,int c,int d,int e)
{
	// debug(a,b,c,d,e,dp[a][b][c][d][e]);
	if(edges[a].size()==0)//???????????????
	{
		assert((b|c|d|e)==0);
		chSegs.push_back(l[a]);
		chSegs.push_back(r[a]);
		return;
	}
	int na,nb,nc,nd,ne;
	pair<ll,ll> &o=dp[a][b][c].trace[d][e];
	int ch=decompressAddress(o.first,na,nb,nc,nd,ne);
	// debug(ch);
	// debug(a,b,c,d,e,na,nb,nc,nd,ne);
	if(o.second>=0)
	{
		// if(a==2 && ch)
		// {
		// 	debug(b,ne,r[edges[a][b-2]]);
		// }
		if(ch && ne==0)
		{
			chSegs.push_back(r[edges[a][b-2]]);
		}
		pushSeg(na,nb,nc,nd,ne);
		decompressAddress(o.second,na,nb,nc,nd,ne);
		if(ch && nd==0)
		{
			chSegs.push_back(l[edges[a][b-1]]);
		}
		pushSeg(na,nb,nc,nd,ne);
	}
	else if(o.second==-llinf-1)
	{
		pushSeg(na,nb,nc,nd,ne);
	}
	else
	{
		// debug("FUCK");
		// assert(na==a);
		// assert(nb==b);
		// assert(mabs(nc-c)<=1);
		int pl,pr;
		// if(o.second==-llinf)
		// {
		// 	debug(a,nd!=d,ne!=e);
		// }
		if(o.second==-llinf) pl=l[a],pr=r[a];
		else pl=pr=-o.second-1;
		if(nd!=d)
		{
			chSegs.push_back(pl);
		}
		if(ne!=e)
		{
			if(o.second==-llinf && a==3) debug(pr);
			chSegs.push_back(pr);
		}
		pushSeg(na,nb,nc,nd,ne);
	}
}

int solve(int x)
{
	getInfo(x);
	dfs(x);
	int val=dp[x][fn[x]][0][0][0];
	pushSeg(x,fn[x],0,0,0);
	return val;
}

int main()
{
//	// freopen("input.txt","r",stdin);
	for(int i=0;i<2;i++) for(int j=0;j<2;j++) vmin[i][j]=-inf;

	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>l[i]>>r[i];
	}
	for(int i=0;i<n;i++)
	{
		pii fa=mp(-inf,-1);
		for(int j=0;j<n;j++)
		{
			if(l[j]<l[i] && r[j]>r[i]) chmax(fa,mp(l[j],j));
		}
		if(fa.second!=-1) edges[fa.second].push_back(i);else roots.push_back(i);
	}
	// debug(roots);

	for(int i=0;i<n;i++)
	{
		sort(all(edges[i]),[](const int &a,const int &b){return l[a]<l[b];});
	}

	vector<pii> allres;
	int res=0;
	for(auto x:roots)
	{
		res+=solve(x);
	}
	cout<<res<<endl;
	sort(all(chSegs));
	assert(chSegs.size()==n*2);
	memset(resl,-1,sizeof(resl));
	for(int i=0;i<n*2;i+=2)
	{
		int lp=chSegs[i];
		int rp=chSegs[i+1];
		pii o=mp(-inf,-1);
		for(int j=0;j<n;j++)
		{
			if(resl[j]==-1 && l[j]<=lp && r[j]>=rp) chmax(o,mp(l[j],j));
		}
		assert(o.second!=-1);
		resl[o.second]=lp;
		resr[o.second]=rp;
	}
	for(int i=0;i<n;i++)
	{
		assert(resl[i]>=l[i] && resr[i]<=r[i]);
		res-=resr[i]-resl[i];
		printf("%d %d\n",resl[i],resr[i]);
	}
	assert(res==0);
	// return res;
	// debug(res);
	// debug(chSegs);
	// for(int i=0;i<n*2;i+=2)
	// {
	// 	printf("%d %d\n",chSegs[i],chSegs[i+1]);
	// }
	return 0;
}