Problem: The Committee for Research on Binary Viruses discovered a method of replication for a large family of viruses whose genetic codes are sequences of zeros and ones. Each virus originates from a single gene; for simplicity genes are denoted by integers from 0 to G - 1. At each moment in time a virus is a sequence of genes. When mutation occurs, one of the genes from the sequence is replaced by a certain sequence of genes, according to the mutation table. The virus stops mutating when it consists only of genes 0 and 1.

For instance, for the following mutation table: $$$ 2 → ⟨ 0\ 1 ⟩ \\\ 3 → ⟨ 2\ 0\ 0⟩\\\ 3 → ⟨ 1\ 3⟩\\\ 4 → ⟨ 0\ 3\ 1\ 2⟩\\\ 5 → ⟨ 2\ 1⟩\\\ 5 → ⟨ 5⟩  a virus that initially consisted of a single gene 4, could have mutated as follows:  ⟨ 4 ⟩ → ⟨ \underline{0\ 3\ 1\ 2} ⟩ → ⟨ 0\ \underline{2\ 0\ 0}\ 1\ 2 ⟩ → ⟨ 0\ \underline{0\ 1}\ 0\ 0\ 1\ 2 ⟩ → ⟨ 0\ 0\ 1\ 0\ 0\ 1\ \underline{0\ 1} ⟩  or in another way:  ⟨ 4 ⟩ → ⟨ \underline{0\ 3\ 1\ 2} ⟩ → ⟨ 0\ \underline{1\ 3}\ 1\ 2 ⟩ → ⟨ 0\ 1\ 3\ 1\ \underline{0\ 1} ⟩ → ⟨ 0\ 1\ \underline{2\ 0\ 0}\ 1\ 0\ 1 ⟩ → ⟨ 0\ 1\ \underline{0\ 1}\ 0\ 0\ 1\ 0\ 1 ⟩ $$$

Viruses are detected by antibodies that identify the presence of specific continuous fragments of zeros and ones in the viruses' codes. For example, an antibody reacting to a fragment ⟨ 0\ 0\ 1\ 0\ 0 ⟩ will detect a virus ⟨ 0\ 0\ 1\ 0\ 0\ 1\ 0\ 1 ⟩, but it will not detect a virus ⟨ 0\ 1\ 0\ 1\ 0\ 0\ 1\ 0\ 1 ⟩.

For each gene from 2 to G-1, the scientists are wondering whether a given set of antibodies is enough to detect all viruses that can emerge through mutations from this gene. If not, they want to know the length of the shortest virus that cannot be detected.

It may happen that sometimes scientists don't have any antibodies. Then of course no virus can be detected, so the scientists are only interested in the length of the shortest possible virus that can emerge from the gene mutations.

Input

The first line of the input will contain three integers G, N and M (G > 2, N ≥ G - 2, M ≥ 0) specifying the number of genes, the number of rows in the mutation table, and the number of antibodies.

The following N lines contain descriptions of rows of the mutation table; each line begins with two integers a and k (2 ≤ a < G, k ≥ 1), followed by a sequence of k integers b_1, b_2, …, b_k (0 ≤ b_i < G), that encode the row $$$ a → ⟨ b_1\ b_2\ …\ b_k ⟩ $$$

The sum of all values k does not exceed 100. Every integer from 2 to G - 1 appears in the table as a at least once.

The next M lines contain descriptions of the antibodies; each such line begins with an integer ℓ (ℓ ≥ 1), followed by a sequence of ℓ integers c_1, c_2, …, c_ℓ (0 ≤ c_i ≤ 1), describing the antibody. The sum of all values ℓ does not exceed 50.

Output

Your program needs to output exactly G - 2 lines, containing the answers for the subsequent genes from 2 to G - 1.

If all viruses that can mutate from this single gene can be detected by the given set of antibodies, you need to print the word "YES". You also need to print this if there are no viruses that could originate from this gene (it can happen when the sequences never stop mutating).

Otherwise you need to print the word "NO", followed by an integer denoting the minimal length of the undetectable virus. You can assume that for all the prepared input data this value will be smaller than 2^{63}.

Scoring

Subtasks: 

  1. (11 points) No antibodies (M = 0) 
  2. (14 points) N = G - 2 
  3. (25 points) One antibody (M = 1) 
  4. (32 points) The sum of all values ℓ does not exceed 10 
  5. (18 points) No further constraints 

Example

Input


6 6 2
2 2 0 1
3 3 2 0 0
3 2 1 3
4 4 0 3 1 2
5 2 2 1
5 1 5
2 1 1
5 0 0 1 0 0


Output


NO 2
NO 4
NO 9
YES

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
bool cur1;
int rd() {
  int res = 0, c, f = 0;
  while (!isdigit(c = getchar())) f = c == '-';
  do res = (res << 1) + (res << 3) + (c ^ 48);
  while (isdigit(c = getchar()));
  return f ? -res : res;
}
int G, n, m;
bool tomin(long long &a, long long b) {
  return (a == -1 || a > b) && ~b && (a = b, 1);
}
namespace p11 {
vector<int> v[105], g[105];
queue<int> q;
int inq[105];
long long dis[105];
void solve() {
  memset(dis, -1, sizeof dis);
  for (int i = 1, k; i <= n; i++) {
    v[i].push_back(rd()), k = rd();
    while (k--) v[i].push_back(rd());
    for (int j = 1; j < (int)v[i].size(); j++) g[v[i][j]].push_back(i);
  }
  memset(dis, -1, sizeof dis);
  dis[0] = dis[1] = inq[0] = inq[1] = 1, q.push(0), q.push(1);
  while (!q.empty()) {
    int x = q.front();
    inq[x] = 0, q.pop();
    for (int j = 0; j < (int)g[x].size(); j++) {
      int t = g[x][j];
      long long cur = 0;
      for (int k = 1; k < (int)v[t].size(); k++)
        if (~dis[v[t][k]] && ~cur)
          cur += dis[v[t][k]];
        else
          cur = -1;
      if (~cur && tomin(dis[v[t][0]], cur) && !inq[v[t][0]])
        q.push(v[t][0]), inq[v[t][0]] = 1;
    }
  }
  for (int i = 2; i < G; i++)
    if (~dis[i])
      printf("NO %lld\n", dis[i]);
    else
      puts("YES");
}
}  // namespace p11
namespace p25 {
struct Mat {
  int n, m;
  long long a[55][55];
  void init(int x, int y) { n = x, m = y, memset(a, -1, sizeof a); }
  Mat operator*(const Mat &b) {
    Mat c;
    c.init(n, b.m);
    for (int k = 0; k <= m; k++)
      for (int i = 0; i <= n; i++)
        if (~a[i][k])
          for (int j = 0; j <= b.m; j++)
            if (~b.a[k][j]) tomin(c.a[i][j], a[i][k] + b.a[k][j]);
    return c;
  }
  Mat operator+(const Mat &b) {
    Mat c;
    c.init(n, m);
    for (int i = 0; i <= n; i++)
      for (int j = 0; j <= m; j++)
        if (~a[i][j] && ~b.a[i][j]) c.a[i][j] = a[i][j] + b.a[i][j];
    return c;
  }
  void Out() {
    for (int i = 0; i <= n; i++, puts(""))
      for (int j = 0; j <= m; j++) printf("%lld ", a[i][j]);
    puts("");
  }
} t[105];
int ind[105];
vector<int> g[105], v[105];
queue<int> q;
int s[105], son[105][2], tot, flg[105], fal[105];
void Insert(int n) {
  int cur = 0;
  for (int i = 1; i <= n; i++) {
    int &nxt = son[cur][s[i]];
    if (!nxt) nxt = ++tot;
    cur = nxt;
  }
  flg[cur] = 1;
}
void init() {
  queue<int> q;
  for (int i = 0; i < 2; i++)
    if (son[0][i]) q.push(son[0][i]), fal[son[0][i]] = 0;
  while (!q.empty()) {
    int x = q.front();
    q.pop();
    for (int i = 0; i < 2; i++)
      if (son[x][i])
        q.push(son[x][i]), fal[son[x][i]] = son[fal[x]][i];
      else
        son[x][i] = son[fal[x]][i];
  }
  for (int i = 0; i < G; i++)
    t[i].init(tot + 1, tot + 1), t[i].a[0][0] = i < 2 ? -1 : 0;
  for (int i = 0; i <= tot; i++) {
    if (flg[son[i][0]]) son[i][0] = tot + 1;
    if (flg[son[i][1]]) son[i][1] = tot + 1;
    t[0].a[i][son[i][0]] = 1;
    t[1].a[i][son[i][1]] = 1;
  }
  t[0].a[tot + 1][tot + 1] = 1;
  t[1].a[tot + 1][tot + 1] = 1;
}
void solve() {
  for (int i = 1; i <= n; i++) {
    int a = rd(), k = rd();
    while (k--) {
      int x = rd();
      v[a].push_back(x);
      g[x].push_back(a);
    }
    ind[a] += v[a].size();
  }
  for (int i = 1; i <= m; i++) {
    int l = rd();
    for (int j = 1; j <= l; j++) s[j] = rd();
    Insert(l);
  }
  init();
  q.push(0), q.push(1);
  while (!q.empty()) {
    int x = q.front();
    q.pop();
    if (x >= 2)
      for (int i = 0; i < (int)v[x].size(); i++) t[x] = t[x] * t[v[x][i]];
    for (int i = 0; i < (int)g[x].size(); i++)
      if (--ind[g[x][i]] == 0) q.push(g[x][i]);
  }
  for (int i = 2; i < G; i++) {
    long long res = -1;
    for (int j = 0; j <= tot; j++) tomin(res, t[i].a[0][j]);
    if (~res)
      printf("NO %lld\n", res);
    else
      puts("YES");
  }
}
}  // namespace p25
namespace p100 {
struct Mat {
  int n, m;
  long long a[55][55];
  void init(int x, int y) { n = x, m = y, memset(a, -1, sizeof a); }
  Mat operator*(const Mat &b) {
    Mat c;
    c.init(n, b.m);
    for (int k = 0; k <= m; k++)
      for (int i = 0; i <= n; i++)
        if (~a[i][k])
          for (int j = 0; j <= b.m; j++)
            if (~b.a[k][j]) tomin(c.a[i][j], a[i][k] + b.a[k][j]);
    return c;
  }
  Mat operator+(const Mat &b) {
    Mat c;
    c.init(n, m);
    for (int i = 0; i <= n; i++)
      for (int j = 0; j <= m; j++) {
        c.a[i][j] = a[i][j];
        tomin(c.a[i][j], b.a[i][j]);
      }
    return c;
  }
  bool operator<(const Mat &b) {
    for (int i = 0; i <= n; i++)
      for (int j = 0; j <= m; j++)
        if (~a[i][j] && (b.a[i][j] == -1 || a[i][j] < b.a[i][j])) return 1;
    return 0;
  }
  void Out() {
    for (int i = 0; i <= n; i++, puts(""))
      for (int j = 0; j <= m; j++) printf("%lld ", a[i][j]);
    puts("");
  }
} t[105];
vector<int> g[105], v[105];
queue<int> q;
int inq[105], vis[105];
int s[105], son[105][2], tot, flg[105], fal[105];
void Insert(int n) {
  int cur = 0;
  for (int i = 1; i <= n; i++) {
    int &nxt = son[cur][s[i]];
    if (!nxt) nxt = ++tot;
    cur = nxt;
  }
  flg[cur] = 1;
}
void init() {
  queue<int> q;
  for (int i = 0; i < 2; i++)
    if (son[0][i]) q.push(son[0][i]), fal[son[0][i]] = 0;
  while (!q.empty()) {
    int x = q.front();
    q.pop();
    for (int i = 0; i < 2; i++)
      if (son[x][i]) {
        flg[son[x][i]] |= flg[x];
        fal[son[x][i]] = son[fal[x]][i];
        q.push(son[x][i]);
      } else
        son[x][i] = son[fal[x]][i];
  }
  for (int i = 1; i <= tot; i++) {
    int x = i;
    while (x) flg[i] |= flg[x], x = fal[x];
  }
  for (int i = 0; i < G; i++) t[i].init(tot, tot);
  for (int i = 0; i <= tot; i++) {
    if (!flg[son[i][0]]) t[0].a[i][son[i][0]] = 1;
    if (!flg[son[i][1]]) t[1].a[i][son[i][1]] = 1;
  }
}
void solve() {
  for (int i = 1, k; i <= n; i++) {
    v[i].push_back(rd()), k = rd();
    while (k--) v[i].push_back(rd());
    for (int j = 0; j < (int)v[i].size(); j++) g[v[i][j]].push_back(i);
  }
  for (int i = 1; i <= m; i++) {
    int l = rd();
    for (int j = 1; j <= l; j++) s[j] = rd();
    Insert(l);
  }
  init();
  q.push(0), q.push(1), inq[0] = inq[1] = vis[0] = vis[1] = 1;
  while (!q.empty()) {
    int x = q.front();
    inq[x] = 0, q.pop();
    for (int i = 0; i < (int)g[x].size(); i++) {
      int p = g[x][i], vv = 1;
      Mat cur;
      cur.init(tot, tot);
      for (int j = 0; j <= tot; j++) cur.a[j][j] = 0;
      for (int k = 1; k < (int)v[p].size(); k++)
        cur = cur * t[v[p][k]], vv &= vis[v[p][k]];
      if (vv && cur < t[v[p][0]]) {
        t[v[p][0]] = t[v[p][0]] + cur, vis[v[p][0]] = 1;
        if (!inq[v[p][0]]) q.push(v[p][0]), inq[v[p][0]] = 1;
      }
    }
  }
  for (int i = 2; i < G; i++) {
    long long res = -1;
    for (int j = 0; j <= tot; j++) tomin(res, t[i].a[0][j]);
    if (~res)
      printf("NO %lld\n", res);
    else
      puts("YES");
  }
}
}  // namespace p100
void my_main() {
  G = rd(), n = rd(), m = rd();
  p100::solve();
}
bool cur2;
int main() {
  my_main();
  return 0;
}