Problem: It is only a few days until Seollal (Korean Lunar New Year), and Jaehyun has invited his family to his garden. There are kids among the guests. To make the gathering more fun for the kids, Jaehyun is going to run a game of hide-and-seek.

The garden can be represented by a n × m grid of unit cells. Some (possibly zero) cells are blocked by rocks, and the remaining cells are free. Two cells are neighbors if they share an edge. Each cell has up to 4 neighbors: two in the horizontal direction and two in the vertical direction. 

Since the garden is represented as a grid, we can classify the cells in the garden as either "black" or "white". The top-left cell is black, and two cells which are neighbors must be different colors. Cell indices are 1-based, so the top-left corner of the garden is cell (1, 1).

Jaehyun wants to turn his garden into a maze by placing some walls between two cells. Walls can only be placed between neighboring cells. If the wall is placed between two neighboring cells a and b, then the two cells a and b are not neighboring from that point. One can walk directly between two neighboring cells if and only if there is no wall directly between them. 

A maze must have the following property. For each pair of free cells in the maze, there must be exactly one simple path between them. A simple path between cells a and b is a sequence of free cells in which the first cell is a, the last cell is b, all cells are distinct, and any two consecutive cells are neighbors which are not directly blocked by a wall.

At first, kids will gather in cell (1, 1), and start the hide-and-seek game. A kid can hide in a cell if and only if that cell is free, it is not (1, 1), and has exactly one free neighbor. Jaehyun planted roses in the black cells, so it's dangerous if the kids hide there. So Jaehyun wants to create a maze where the kids can only hide in white cells.

You are given the map of the garden as input. Your task is to help Jaehyun create a maze.

Input

Your program will be judged in multiple test cases.

The first line contains the number of test cases t. (1 ≤ t ≤ 100). Afterward, t test cases with the described format will be given.

The first line of a test contains two integers n, m (2 ≤ n, m ≤ 20), the size of the grid.

In the next n line of a test contains a string of length m, consisting of the following characters (without any whitespace): 

  * O: A free cell. 
  * X: A rock. 



It is guaranteed that the first cell (cell (1, 1)) is free, and every free cell is reachable from (1, 1). 

If t ≥ 2 is satisfied, then the size of the grid will satisfy n ≤ 10, m ≤ 10. In other words, if any grid with size n > 10 or m > 10 is given as an input, then it will be the only input on the test case (t = 1).

Output

For each test case, print the following:

If there are no possible mazes, print a single line NO.

Otherwise, print a single line YES, followed by a grid of size (2n-1) × (2m-1) denoting the found maze. The rules for displaying the maze follows. All cells are indexed in 1-base.

  * For all 1 ≤ i ≤ n, 1 ≤ j ≤ m, if the cell (i, j) is free cell, print 'O' in the cell (2i-1, 2j-1). Otherwise, print 'X' in the cell (2i-1, 2j-1). 
  * For all 1 ≤ i ≤ n, 1 ≤ j ≤ m-1, if the neighboring cell (i, j), (i, j+1) have wall blocking it, print ' ' in the cell (2i-1, 2j). Otherwise, print any printable character except spaces in the cell (2i-1, 2j). A printable character has an ASCII code in range [32, 126]: This includes spaces and alphanumeric characters. 
  * For all 1 ≤ i ≤ n-1, 1 ≤ j ≤ m, if the neighboring cell (i, j), (i+1, j) have wall blocking it, print ' ' in the cell (2i, 2j-1). Otherwise, print any printable character except spaces in the cell (2i, 2j-1) 
  * For all 1 ≤ i ≤ n-1, 1 ≤ j ≤ m-1, print any printable character in the cell (2i, 2j). 



Please, be careful about trailing newline characters or spaces. Each row of the grid should contain exactly 2m-1 characters, and rows should be separated by a newline character. Trailing spaces must not be omitted in a row.

Example

Input


4
2 2
OO
OO
3 3
OOO
XOO
OOO
4 4
OOOX
XOOX
OOXO
OOOO
5 6
OOOOOO
OOOOOO
OOOOOO
OOOOOO
OOOOOO


Output


YES
OOO
  O
OOO
NO
YES
OOOOO X
  O O  
X O O X
  O    
OOO X O
O O   O
O OOOOO
YES
OOOOOOOOOOO
  O   O   O
OOO OOO OOO
O   O   O  
OOO OOO OOO
  O   O   O
OOO OOO OOO
O   O   O  
OOO OOO OOO

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
struct MI {
 private:
  char bb[1 << 14];
  FILE* f;
  char *bs, *be;
  char e;
  bool o, l;

 public:
  MI() : f(stdin), bs(0), be(0) {}
  MI(FILE* f) : f(f), bs(0), be(0) {}
  inline operator bool() { return !l; }
  inline char get() {
    if (o) {
      o = 0;
      return e;
    }
    if (bs == be) be = (bs = bb) + fread(bb, 1, sizeof(bb), f);
    if (bs == be) {
      l = 1;
      return -1;
    };
    return *bs++;
  }
  inline void unget(char c) {
    o = 1;
    e = c;
  }
  template <class T>
  inline T read() {
    T r;
    *this > r;
    return r;
  }
  template <class T>
  inline MI& operator>(T&);
};
template <class T>
struct Q {
  const static bool U = T(-1) >= T(0);
  inline void operator()(MI& t, T& r) const {
    r = 0;
    char c;
    bool y = 0;
    if (U)
      for (;;) {
        c = t.get();
        if (c == -1) goto E;
        if (isdigit(c)) break;
      }
    else
      for (;;) {
        c = t.get();
        if (c == -1) goto E;
        if (c == '-') {
          c = t.get();
          if (isdigit(c)) {
            y = 1;
            break;
          };
        } else if (isdigit(c))
          break;
        ;
      };
    for (;;) {
      if (c == -1) goto E;
      if (isdigit(c))
        r = r * 10 + (c ^ 48);
      else
        break;
      c = t.get();
    }
    t.unget(c);
  E:;
    if (y) r = -r;
  }
};
template <>
struct Q<char> {
  inline void operator()(MI& t, char& r) {
    int c;
    for (;;) {
      c = t.get();
      if (c == -1) {
        r = -1;
        return;
      };
      if (!isspace(c)) {
        r = c;
        return;
      };
    }
  }
};
template <class T>
inline MI& MI::operator>(T& t) {
  Q<T>()(*this, t);
  return *this;
}
template <class T>
std::ostream& operator<(std::ostream& out, const T& t) {
  return out << t;
}
using std::cout;
MI cin;
const int $n = 25;
const int $p = $n * $n;
const int $e = 770;
int tp;
bool has[$e];
void extend_init();
bool can_add_A(int x);
bool can_add_B(int x);
void pre_A(int old);
bool ok_for_A(int now);
void pre_B(int now);
bool ok_for_B(int old);
void on_extend();
inline void matroid_intersection() {
  static int q[$e], fr[$e];
  static bool fin[$e];
  int hd, tl;
  auto push = [&](int x, int y) -> bool {
    fr[y] = x;
    q[++tl] = y;
    if (fin[y]) {
      for (int c = y; ~c; c = fr[c]) has[c] ^= 1;
      return 1;
    }
    return 0;
  };
  auto extend = [&]() -> bool {
    extend_init();
    hd = 1;
    tl = 0;
    for (int __r = (tp), i = (1); i <= __r; ++i) {
      fr[i] = fin[i] = 0;
      if (has[i]) continue;
      const auto A = can_add_A(i), B = can_add_B(i);
      if (A && B) return has[i] = 1;
      if (A)
        fr[q[++tl] = i] = -1;
      else
        fin[i] = B;
    }
    while (hd <= tl) {
      const auto x = q[hd++];
      if (has[x]) {
        pre_A(x);
        for (int __r = (tp), y = (1); y <= __r; ++y)
          if (!has[y] && !fr[y] && ok_for_A(y) && push(x, y)) return 1;
      } else {
        pre_B(x);
        for (int __r = (tp), y = (1); y <= __r; ++y)
          if (has[y] && !fr[y] && ok_for_B(y) && push(x, y)) return 1;
      }
    }
    return 0;
  };
  while (extend()) on_extend();
}
struct point {
  int x, y;
} pos[$e];
int n, m, dn, dm, es[$e][2], black[$e], id[$n][$n], deg[$p], f[$p], ps, add;
bool bl[$n][$n];
char out[$n * 2][$n * 2];
int getf(int x) { return f[x] == x ? x : (f[x] = getf(f[x])); }
inline void build() {
  for (int __r = (ps), i = (1); i <= __r; ++i) f[i] = i;
  for (int __r = (tp), i = (1); i <= __r; ++i)
    if (has[i]) f[getf(es[i][0])] = getf(es[i][1]);
}
void extend_init() {
  build();
  for (int __r = (ps), i = (1); i <= __r; ++i) deg[i] = 0;
  for (int __r = (tp), i = (1); i <= __r; ++i)
    if (has[i]) ++deg[es[i][0]], ++deg[es[i][1]];
}
bool can_add_A(int x) { return getf(es[x][0]) != getf(es[x][1]); }
bool can_add_B(int x) { return black[x] == 1 || deg[black[x]] != 2; }
void pre_A(int old) {
  has[old] = 0;
  build();
  has[old] = 1;
}
bool ok_for_A(int now) { return getf(es[now][0]) != getf(es[now][1]); }
void pre_B(int now) { add = black[now]; }
bool ok_for_B(int old) {
  return black[old] == add || add == 1 || deg[add] != 2;
}
void on_extend() {}
inline void work() {
  cin > n > m;
  dn = n * 2 - 1;
  dm = m * 2 - 1;
  for (int __r = (dn), i = (1); i <= __r; ++i)
    for (int __r = (dm), j = (1); j <= __r; ++j) out[i][j] = ' ';
  ps = 0;
  for (int __r = (n), i = (1); i <= __r; ++i) {
    cin.unget(cin.read<char>());
    for (int __r = (m), j = (1); j <= __r; ++j) {
      bl[i][j] = (out[i * 2 - 1][j * 2 - 1] = cin.get()) == 'X';
      if (!bl[i][j])
        id[i][j] = ++ps;
      else
        id[i][j] = 0;
    }
  }
  tp = 0;
  for (int __r = (n), i = (1); i <= __r; ++i)
    for (int __r = ((m)-1), j = ((i == 1) + 1); j <= __r; ++j)
      if (id[i][j] && id[i][j + 1])
        pos[++tp] = {i * 2 - 1, j * 2}, es[tp][0] = id[i][j],
        es[tp][1] = id[i][j + 1], black[tp] = es[tp][(i ^ j) & 1], has[tp] = 0;
  for (int __r = (m), j = (1); j <= __r; ++j)
    for (int __r = ((n)-1), i = ((j == 1) + 1); i <= __r; ++i)
      if (id[i][j] && id[i + 1][j])
        pos[++tp] = {i * 2, j * 2 - 1}, es[tp][0] = id[i][j],
        es[tp][1] = id[i + 1][j], black[tp] = es[tp][(i ^ j) & 1], has[tp] = 0;
  matroid_intersection();
  for (int __r = (ps), i = (1); i <= __r; ++i) deg[i] = 0;
  for (int __r = (tp), i = (1); i <= __r; ++i)
    if (has[i]) ++deg[black[i]];
  for (int __r = (n), i = (1); i <= __r; ++i)
    for (int __r = (m), j = (1); j <= __r; ++j)
      if (id[i][j] && (i != 1 || j != 1) && !((i ^ j) & 1) &&
          deg[id[i][j]] != 2)
        return (void)(cout < "NO\n");
  if (id[1][2])
    pos[++tp] = {1, 2}, es[tp][0] = 1, es[tp][1] = id[1][2], has[tp] = 0;
  if (id[2][1])
    pos[++tp] = {2, 1}, es[tp][0] = 1, es[tp][1] = id[2][1], has[tp] = 0;
  build();
  for (int __r = (tp), i = (1); i <= __r; ++i)
    if (getf(es[i][0]) != getf(es[i][1]))
      has[i] = 1, f[getf(es[i][0])] = f[getf(es[i][1])];
  for (int __r = (tp), i = (1); i <= __r; ++i)
    if (has[i]) out[pos[i].x][pos[i].y] = 'O';
  cout < "YES\n";
  for (int __r = (dn), i = (1); i <= __r; ++i) {
    for (int __r = (dm), j = (1); j <= __r; ++j) cout < out[i][j];
    cout < ('\n');
  }
}
int main() {
  for (int __r = ((cin.read<int>())), _ = (1); _ <= __r; ++_) work();
}