#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pb push_back
#define ff first
#define ss second
#define _ << " " <<
#define yes cout<<"YES\n"
#define no cout<<"NO\n"
#define ull unsigned long long
#define lll __int128
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
#define BlueCrowner ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#define FOR(i, a, b) for (ll i = (a); i < (b); i++)
#define FORD(i, a, b) for (ll i = (a); i >= (b); i--)
const ll mod = 1e9 + 7;
const ll mod1 = 998244353;
const ll naim = 1e9;
const ll max_bit = 60;
const ull tom = ULLONG_MAX;
const ll MAXN = 100005;
const ll LOG = 20;
const ll NAIM = 1e18;
const ll N = 2e6 + 5;
// ---------- GCD ----------
ll gcd(ll a, ll b) {
while (b) {
a %= b;
swap(a, b);
}
return a;
}
// ---------- LCM ----------
ll lcm(ll a, ll b) {
return a / gcd(a, b) * b;
}
// ---------- Modular Exponentiation ----------
ll modpow(ll a, ll b, ll m = mod) {
ll c = 1;
a %= m;
while (b > 0) {
if (b & 1) c = c * a % m;
a = a * a % m;
b >>= 1;
}
return c;
}
// ---------- Modular Inverse (Fermat’s Little Theorem) ----------
ll modinv(ll a, ll m = mod) {
return modpow(a, m - 2, m);
}
// ---------- Factorials and Inverse Factorials ----------
ll fact[N], invfact[N];
void pre_fact(ll n = N-1, ll m = mod) {
fact[0] = 1;
for (ll i = 1; i <= n; i++) fact[i] = fact[i-1] * i % m;
invfact[n] = modinv(fact[n], m);
for (ll i = n; i > 0; i--) invfact[i-1] = invfact[i] * i % m;
}
// ---------- nCr ----------
ll nCr(ll n, ll r, ll m = mod) {
if (r < 0 || r > n) return 0;
return fact[n] * invfact[r] % m * invfact[n-r] % m;
}
// ---------- Sieve of Eratosthenes ----------
vector<ll> primes;
bool is_prime[N];
void sieve(ll n = N-1) {
fill(is_prime, is_prime + n + 1, true);
is_prime[0] = is_prime[1] = false;
for (ll i = 2; i * i <= n; i++) {
if (is_prime[i]) {
for (ll j = i * i; j <= n; j += i)
is_prime[j] = false;
}
}
for (ll i = 2; i <= n; i++)
if (is_prime[i]) primes.pb(i);
}
void solve() {
ll n; cin >> n;
vector<vector<ll>> g(n + 1);
FOR(i, 0, n - 1) {
ll u, v; cin >> u >> v;
g[u].pb(v);
g[v].pb(u);
}
vector<ll> sz(n + 1, 0);
ll ans1 = 0;
function<void(ll, ll)> dfs_size = [&](ll u, ll par) {
sz[u] = 1;
for(auto &v : g[u]) {
if(v == par) continue;
dfs_size(v, u);
ans1 += min(sz[v], n - sz[v]);
sz[u] += sz[v];
}
};
dfs_size(1, 0);
function<ll(ll, ll)> find_centroid = [&](ll u, ll par) -> ll {
for(auto &v : g[u]) {
if(v == par) continue;
if(sz[v] > n / 2) return find_centroid(v, u);
}
return u;
};
ll c = find_centroid(1, 0);
vector<ll> ord(n + 1);
ll cur = 1;
function<void(ll, ll)> dfs = [&](ll u, ll par) {
ord[cur] = u;
for(auto &v : g[u]) {
if(v == par) continue;
cur++;
dfs(v, u);
}
};
dfs(c, 0);
vector<ll> res1(n + 1, 0);
iota(all(res1), 0);
ord[0] = ord[n];
vector<ll> t(n + 1, 0);
iota(all(t), 0);
FOR(i, 1, n + 1) {
res1[ord[i]] = t[ord[(i + n / 2) % n]];
}
ans1 *= 2;
vector<ll> res2(n + 1, 0);
iota(all(res2), 0);
ll ans2 = 0;
function<void(ll, ll)> dfs2 = [&](ll u, ll par) {
for(auto &v : g[u]) {
if(v == par) continue;
dfs2(v, u);
}
if(res2[u] == u) {
swap(res2[u], res2[par]);
ans2 += 2;
}
};
dfs2(1, 0);
if(res2[1] == 0) {
res2[1] = 1;
swap(res2[res2[g[1][0]]], res2[1]);
}
cout << ans2 _ ans1 << '\n';
FOR(i, 1, n + 1) cout << res2[i] << ' '; cout << '\n';
FOR(i, 1, n + 1) cout << res1[i] << ' '; cout << '\n';
}
int main() {
BlueCrowner;
ll t = 1;
//cin >> t;
while (t--) {
solve();
}
return 0;
}
