#include <bits/stdc++.h>
using namespace std;

const int maxn = 4010;
int n, m, tot, s[maxn], t[maxn], cnt[maxn][maxn];
int ans, fa[maxn], dep[maxn], deg[maxn], dist[maxn][maxn];
bool rev[maxn], chk[maxn], chke[maxn];
vector<int> G[maxn], path[maxn];
vector<pair<int, int>> cov[maxn];

int main() {
    scanf("%d %d", &n, &m);
    for (int i = 1, u, v; i < n; i++) {
        scanf("%d %d", &u, &v), deg[u]++, deg[v]++;
        G[u].push_back(v), G[v].push_back(u);
    }
    function<void(int)> dfs = [&](int v) {
        for (int u : G[v]) if (u ^ fa[v]) {
            fa[u] = v, dep[u] = dep[v] + 1, dfs(u);
        }
    };
    dfs(1);
    for (int i = 1; i <= n; i++) {
        function<void(int, int)> dfs = [&](int v, int fa) {
            for (int u : G[v]) if (u ^ fa) {
                dist[i][u] = dist[i][v] + 1, dfs(u, v);
            }
        };
        dfs(i, 0);
    }
    auto add = [&](int id, int u, int v) {
        s[id] = u, t[id] = v;
        if (u == v) return;
        vector<int> L, R;
        while (u ^ v) {
            if (dep[u] > dep[v]) L.push_back(u), u = fa[u];
            else R.push_back(v), v = fa[v];
        }
        for (int v : L) path[id].push_back(v);
        path[id].push_back(u), reverse(R.begin(), R.end());
        for (int v : R) path[id].push_back(v);
        for (int i = 0; i < path[id].size(); i++) {
            cov[path[id][i]].push_back({id, i});
            if (i == path[id].size() - 1) continue;
            int u = path[id][i], v = path[id][i + 1];
            if (u > v) swap(u, v);
            if (cnt[u][v]++ < 2) ans++;
        }
    };
    for (int i = 1, u, v; i <= m; i++) {
        scanf("%d %d", &u, &v), add(++tot, u, v);
    }
    function<void()> solve = [&]() {
        int u = 0, v = 0;
        for (int i = 1; i <= n; i++) {
            if (!chk[i] && deg[i] == 1) { u = i; break; }
        }
        if (!u) return;
        for (int i : G[u]) {
            if (!chk[i]) { v = i; break; }
        }
        chk[u] = 1, deg[v]--;
        auto query = [&](int v, int pos) {
            return (~pos && pos < path[v].size()) ? path[v][pos] : 0;
        };
        auto get = [&](int v, int pos, int t) {
            return query(v, pos - 1) == t || query(v, pos + 1) == t;
        };
        pair<int, int> pu = {0, 0}, pv = {0, 0};
        for (auto p : cov[u]) {
            if (!chke[p.first] && get(p.first, p.second, v)) {
                if (!pu.first) pu = p;
                else if (!pv.first) pv = p;
            }
        }
        if (!pv.first) { solve(); return; }
        auto reverse = [&](int v, int x, int y) {
            if (rev[v]) swap(x, y);
            if ((dist[s[v]][x] + dist[y][t[v]] + 1) == dist[s[v]][t[v]]) rev[v] ^= 1;
        };
        reverse(pu.first, u, v), reverse(pv.first, v, u);
        chke[pu.first] = chke[pv.first] = 1;
        auto find = [&](int v, int pos, int x) {
            return query(v, pos + 1) ^ x ? path[v][0] : path[v].back();
        };
        add(++tot, find(pu.first, pu.second, v), find(pv.first, pv.second, v));
        int tmp = tot; solve();
        rev[pu.first] ^= rev[tmp], rev[pv.first] ^= rev[tmp];
    };
    printf("%d\n", ans), solve();
    for (int i = 1; i <= m; i++) {
        if (rev[i]) printf("%d %d\n", t[i], s[i]);
        else printf("%d %d\n", s[i], t[i]);
    }
    return 0;
}

Measured memory_bytes: 95616