#include <bits/stdc++.h>
using namespace std;
const int mm = 100005;
const int mod = 1000000007;
using namespace std;
vector<int> G[mm];
vector<int> c[20][mm];
int par[20][mm];
int d[mm];
int n, m, q;
int lg[mm];
void dfs(int rt) {
  for (int i = 0; 1 << i < d[rt]; i++) {
    int &fa = par[i][rt];
    if (fa) par[i + 1][rt] = par[i][fa];
    vector<int> &vec = c[i + 1][rt];
    vector<int> &vec1 = c[i][rt];
    vector<int> &vec2 = c[i][fa];
    vec.resize(vec1.size() + vec2.size());
    merge(vec1.begin(), vec1.end(), vec2.begin(), vec2.end(), vec.begin());
    if (vec.size() > 10) vec.resize(10);
  }
  for (int i = 0; i < G[rt].size(); i++)
    if (G[rt][i] != par[0][rt]) {
      int &u = G[rt][i];
      par[0][u] = rt;
      d[u] = d[rt] + 1;
      dfs(u);
    }
}
int LCA(int u, int v) {
  while (d[u] != d[v]) {
    if (d[u] > d[v]) swap(u, v);
    int i = lg[d[v] - d[u]];
    v = par[i][v];
  }
  if (u == v) return v;
  for (int t = lg[d[u]]; t >= 0; t) {
    int i = 0;
    while (par[i][u] != par[i][v]) i++;
    if (i == 0) return par[0][u];
    u = par[i - 1][u];
    v = par[i - 1][v];
  }
  return v;
}
vector<int> get(int u, int up, int a, vector<int> &v) {
  vector<int> res[2];
  vector<int> *pre = &res[0], *cur = &res[1];
  while (up > 0) {
    int i = lg[up];
    vector<int> &vec = c[i][u];
    pre->resize(cur->size() + vec.size());
    merge(cur->begin(), cur->end(), vec.begin(), vec.end(), pre->begin());
    if (pre->size() > a) pre->resize(a);
    swap(cur, pre);
    u = par[i][u];
    up -= 1 << i;
  }
  v = *cur;
  return *cur;
}
void solve(int u, int v, int a) {
  int anc = LCA(u, v);
  vector<int> v1;
  get(u, d[u] - d[anc] + 1, a, v1);
  vector<int> v2;
  get(v, d[v] - d[anc], a, v2);
  vector<int> v3(v1.size() + v2.size());
  merge(v1.begin(), v1.end(), v2.begin(), v2.end(), v3.begin());
  a = min(a, (int)v3.size());
  printf("%d", a);
  for (int i = 0; i < a; i++) {
    printf(" %d", v3[i]);
  }
  cout << endl;
}
vector<int> A;
int main() {
  cin >> n >> m >> q;
  for (int i = 1; i < n; i++) {
    int x, y;
    scanf("%d%d", &x, &y);
    G[x].push_back(y);
    G[y].push_back(x);
  }
  for (int i = 0; i < m; i++) {
    int x;
    scanf("%d", &x);
    c[0][x].push_back(i + 1);
  }
  for (int i = 1; i <= n; i++) {
    sort(c[0][i].begin(), c[0][i].end());
    if (c[0][i].size() > 10) c[0][i].resize(10);
  }
  lg[0] = -1;
  for (int i = 1; i <= n; i++)
    lg[i] = (i & (i - 1)) ? lg[i - 1] : lg[i - 1] + 1;
  int root = 1;
  d[root] = 1;
  dfs(root);
  for (int i = 0; i < q; i++) {
    int u, v, a;
    scanf("%d%d%d", &u, &v, &a);
    solve(u, v, a);
  }
  return 0;
}
