Problem: You are given a binary matrix a of size n × m. A binary matrix is a matrix where each element is either 0 or 1.

You may perform some (possibly zero) operations with this matrix. During each operation you can inverse the row of this matrix or a column of this matrix. Formally, inverting a row is changing all values in this row to the opposite (0 to 1, 1 to 0). Inverting a column is changing all values in this column to the opposite.

Your task is to sort the initial matrix by some sequence of such operations. The matrix is considered sorted if the array [a_{1, 1}, a_{1, 2}, ..., a_{1, m}, a_{2, 1}, a_{2, 2}, ..., a_{2, m}, ..., a_{n, m - 1}, a_{n, m}] is sorted in non-descending order.

Input

The first line of the input contains two integers n and m (1 ≤ n, m ≤ 200) — the number of rows and the number of columns in the matrix.

The next n lines contain m integers each. The j-th element in the i-th line is a_{i, j} (0 ≤ a_{i, j} ≤ 1) — the element of a at position (i, j).

Output

If it is impossible to obtain a sorted matrix, print "NO" in the first line.

Otherwise print "YES" in the first line. In the second line print a string r of length n. The i-th character r_i of this string should be '1' if the i-th row of the matrix is inverted and '0' otherwise. In the third line print a string c of length m. The j-th character c_j of this string should be '1' if the j-th column of the matrix is inverted and '0' otherwise. If there are multiple answers, you can print any.

Examples

Input


2 2
1 1
0 1


Output


YES
00
10


Input


3 4
0 0 0 1
0 0 0 0
1 1 1 1


Output


YES
010
0000


Input


3 3
0 0 0
1 0 1
1 1 0


Output


NO

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
vector<vector<int> > solve(vector<vector<int> > a) {
  int n = int(a.size());
  int m = int(a[0].size());
  vector<vector<int> > ca = a;
  for (int i = 0; i < n; ++i) {
    int row_first = ca[i][0];
    for (int j = 0; j < m; ++j) {
      ca[i][j] ^= row_first;
    }
  }
  vector<vector<int> > uca = ca;
  sort(uca.begin(), uca.end());
  uca.resize(int(unique(uca.begin(), uca.end()) - uca.begin()));
  if (uca.size() > 2) return vector<vector<int> >();
  if (uca.size() == 1) {
    vector<vector<int> > sol;
    vector<int> rf(n);
    for (int i = 0; i < n; ++i) rf[i] = a[i][0] ^ a[0][0];
    sol.push_back(rf);
    sol.push_back(a[0]);
    return sol;
  }
  for (int d = 0; d < int(ca.size()); ++d) {
    int ind = int(find(uca.begin(), uca.end(), ca[d]) - uca.begin());
    if (ind == int(uca.size())) continue;
    int cnt = 0;
    for (int i = 0; i < n && cnt < 2; ++i) {
      if (ca[i] == uca[ind]) ++cnt;
    }
    uca.erase(uca.begin() + ind);
    if (cnt > 1) continue;
    vector<int> cf = ca[(d + 1) % n];
    vector<int> r = a[d];
    int rf_d = a[d][0] ^ cf[0];
    for (int i = 0; i < m; ++i) {
      r[i] = a[d][i] ^ rf_d ^ cf[i];
    }
    vector<int> sr = r;
    sort(sr.begin(), sr.end());
    if (r != sr) continue;
    vector<int> rf(n);
    for (int i = 0; i < d; ++i) {
      rf[i] = a[i][0] ^ cf[0];
    }
    rf[d] = rf_d;
    for (int i = d + 1; i < n; ++i) {
      rf[i] = a[i][0] ^ cf[0] ^ 1;
    }
    vector<vector<int> > sol;
    sol.push_back(rf);
    sol.push_back(cf);
    return sol;
  }
  return vector<vector<int> >();
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  int n, m;
  cin >> n >> m;
  vector<vector<int> > a(n, vector<int>(m));
  for (int i = 0; i < n; ++i) {
    for (int j = 0; j < m; ++j) cin >> a[i][j];
  }
  vector<vector<int> > sol = solve(a);
  if (sol.size() == 0)
    cout << "NO" << endl;
  else {
    assert(sol.size() == 2);
    cout << "YES" << endl;
    for (int dim = 0; dim < 2; ++dim) {
      for (int i = 0; i < int(sol[dim].size()); ++i) {
        cout << sol[dim][i];
      }
      cout << endl;
    }
  }
}