#include <bits/stdc++.h>
using namespace std;
int n, k;
const int mx = 2 * 100007;
vector<int> adj[mx];
bool T[mx];
int S[mx];
int sons(int v, int p) {
  S[v] = T[v];
  for (auto x : adj[v])
    if (x != p) S[v] += sons(x, v);
  return S[v];
}
int cen = 0;
void dfs(int v, int p) {
  bool R = 1;
  for (auto x : adj[v]) {
    if (x == p) continue;
    R &= (S[x] <= k);
    dfs(x, v);
  }
  if (p != 0) R &= ((2 * k - S[v]) <= k);
  if (R) cen = v;
}
vector<int> SS[mx];
int y = 0;
vector<pair<int, int> > res;
priority_queue<pair<int, int> > F;
void dfs1(int v, int p) {
  if (T[v]) SS[y].push_back(v);
  if (v == cen && T[v]) {
    F.push({SS[y].size(), y});
    y++;
  }
  for (auto x : adj[v]) {
    if (x == p) continue;
    dfs1(x, v);
    if (v == cen && !SS[y].empty()) {
      F.push({SS[y].size(), y});
      y++;
    }
  }
}
int main() {
  scanf("%d %d", &n, &k);
  for (int i = 0; i < n - 1; i++) {
    int u, v;
    scanf("%d %d", &u, &v);
    adj[u].push_back(v);
    adj[v].push_back(u);
  }
  fill(T, T + n + 1, false);
  for (int i = 0; i < 2 * k; i++) {
    int u;
    scanf("%d", &u);
    T[u] = 1;
  }
  sons(1, 0);
  dfs(1, 0);
  printf("1\n%d\n", cen);
  y = 0;
  dfs1(cen, 0);
  while (!F.empty()) {
    int x1 = F.top().second;
    int y1 = F.top().first;
    F.pop();
    int x2 = F.top().second;
    int y2 = F.top().first;
    F.pop();
    res.push_back({SS[x1].back(), SS[x2].back()});
    SS[x1].pop_back();
    SS[x2].pop_back();
    if (!SS[x1].empty()) F.push({SS[x1].size(), x1});
    if (!SS[x2].empty()) F.push({SS[x2].size(), x2});
  }
  for (int i = 0; i < res.size(); i++)
    printf("%d %d %d\n", res[i].first, res[i].second, cen);
  return 0;
}
