#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 2;
int N, sz[MAXN];
vector<int> adj[MAXN], ord, neword;
long long sum;
int dfs(int u, int pa) {
  sz[u] = 1;
  for (int v : adj[u]) {
    if (v != pa) {
      sz[u] += dfs(v, u);
    }
  }
  return sz[u];
}
int getCentroid(int u, int p, int tot) {
  for (int v : adj[u]) {
    if (v != p)
      if (sz[v] * 2 > tot) {
        return getCentroid(v, u, tot);
      }
  }
  return u;
}
void dfs2(int u, int pa) {
  ord.push_back(u);
  for (int v : adj[u]) {
    if (v != pa) {
      dfs2(v, u);
    }
  }
}
int main() {
  cin.sync_with_stdio(0);
  cin.tie(0);
  cin >> N;
  neword.resize(N, 0);
  for (int i = 1, u, v; i < N; i++) {
    cin >> u >> v;
    adj[u].push_back(v);
    adj[v].push_back(u);
  }
  int C = getCentroid(1, 0, dfs(1, 0));
  dfs2(C, 0);
  for (int i = 1; i <= N; i++) sum += 2 * min(sz[i], N - sz[i]);
  cout << sum << '\n';
  for (int i = 0; i < N; i++) {
    neword[ord[i]] = ord[(i + N / 2) % N];
  }
  for (int i = 1; i <= N; i++) {
    cout << neword[i] << ' ';
  }
  return 0;
}
