#include <algorithm>
#include <array>
#include <bitset>
#include <cassert>
#include <chrono>
#include <cstring>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <vector>

using namespace std;

// BEGIN NO SAD
#define rep(i, a, b) for(int i = a; i < (b); ++i)
#define trav(a, x) for(auto& a : x)
#define all(x) x.begin(), x.end()
#define sz(x) (int)(x).size()
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define lb lower_bound
#define ub upper_bound
typedef vector<int> vi;
#define f first
#define s second
#define derr if(1) cerr

void __print(int x) {cerr << x;}
void __print(long x) {cerr << x;}
void __print(long long x) {cerr << x;}
void __print(unsigned x) {cerr << x;}
void __print(unsigned long x) {cerr << x;}
void __print(unsigned long long x) {cerr << x;}
void __print(float x) {cerr << x;}
void __print(double x) {cerr << x;}
void __print(long double x) {cerr << x;}
void __print(char x) {cerr << '\'' << x << '\'';}
void __print(const char *x) {cerr << '\"' << x << '\"';}
void __print(const string &x) {cerr << '\"' << x << '\"';}
void __print(bool x) {cerr << (x ? "true" : "false");}
 
template<typename T, typename V>
void __print(const pair<T, V> &x) {cerr << '{'; __print(x.first); cerr << ", "; __print(x.second); cerr << '}';}
template<typename T>
void __print(const T &x) {int f = 0; cerr << '{'; for (auto &i: x) cerr << (f++ ? ", " : ""), __print(i); cerr << "}";}
void _print() {cerr << "]\n";}
template <typename T, typename... V>
void _print(T t, V... v) {__print(t); if (sizeof...(v)) cerr << ", "; _print(v...);}
#define debug(x...) cerr << "\e[91m"<<__func__<<":"<<__LINE__<<" [" << #x << "] = ["; _print(x); cerr << "\e[39m" << flush;
// END NO SAD

template<class Fun>
class y_combinator_result {
  Fun fun_;
public:
  template<class T>
  explicit y_combinator_result(T &&fun): fun_(std::forward<T>(fun)) {}

  template<class ...Args>
  decltype(auto) operator()(Args &&...args) {
    return fun_(std::ref(*this), std::forward<Args>(args)...);
  }
};

template<class Fun>
decltype(auto) y_combinator(Fun &&fun) {
  return y_combinator_result<std::decay_t<Fun>>(std::forward<Fun>(fun));
}

template<class T>
bool updmin(T& a, T b) {
  if(b < a) {
    a = b;
    return true;
  }
  return false;
}
template<class T>
bool updmax(T& a, T b) {
  if(b > a) {
    a = b;
    return true;
  }
  return false;
}
typedef int64_t ll;

struct barrett_reduction {
    unsigned mod;
    uint64_t div;
 
    barrett_reduction(unsigned m) : mod(m), div(-1LLU / m) {}
 
    unsigned operator()(uint64_t a) const {
#ifdef __SIZEOF_INT128__
        uint64_t q = uint64_t(__uint128_t(div) * a >> 64);
        uint64_t r = a - q * mod;
        return unsigned(r < mod ? r : r - mod);
#endif
        return unsigned(a % mod);
    }
};
 
template<const int &MOD, const barrett_reduction &barrett>
struct _b_int {
    int val;
 
    _b_int(int64_t v = 0) {
        if (v < 0) v = v % MOD + MOD;
        if (v >= MOD) v %= MOD;
        val = int(v);
    }
 
    _b_int(uint64_t v) {
        if (v >= uint64_t(MOD)) v %= MOD;
        val = int(v);
    }
 
    _b_int(int v) : _b_int(int64_t(v)) {}
    _b_int(unsigned v) : _b_int(uint64_t(v)) {}
 
    static int inv_mod(int a, int m = MOD) {
        // https://en.wikipedia.org/wiki/Extended_Euclidean_algorithm#Example
        int g = m, r = a, x = 0, y = 1;
 
        while (r != 0) {
            int q = g / r;
            g %= r; swap(g, r);
            x -= q * y; swap(x, y);
        }
 
        return x < 0 ? x + m : x;
    }
 
    explicit operator int() const { return val; }
    explicit operator unsigned() const { return val; }
    explicit operator int64_t() const { return val; }
    explicit operator uint64_t() const { return val; }
    explicit operator double() const { return val; }
    explicit operator long double() const { return val; }
 
    _b_int& operator+=(const _b_int &other) {
        val -= MOD - other.val;
        if (val < 0) val += MOD;
        return *this;
    }
 
    _b_int& operator-=(const _b_int &other) {
        val -= other.val;
        if (val < 0) val += MOD;
        return *this;
    }
 
    static unsigned fast_mod(uint64_t x) {
#if !defined(_WIN32) || defined(_WIN64)
        return barrett(x);
#endif
        // Optimized mod for Codeforces 32-bit machines.
        // x must be less than 2^32 * MOD for this to work, so that x / MOD fits in an unsigned 32-bit int.
        unsigned x_high = unsigned(x >> 32), x_low = unsigned(x);
        unsigned quot, rem;
        asm("divl %4\n"
            : "=a" (quot), "=d" (rem)
            : "d" (x_high), "a" (x_low), "r" (MOD));
        return rem;
    }
 
    _b_int& operator*=(const _b_int &other) {
        val = fast_mod(uint64_t(val) * other.val);
        return *this;
    }
 
    _b_int& operator/=(const _b_int &other) {
        return *this *= other.inv();
    }
 
    friend _b_int operator+(const _b_int &a, const _b_int &b) { return _b_int(a) += b; }
    friend _b_int operator-(const _b_int &a, const _b_int &b) { return _b_int(a) -= b; }
    friend _b_int operator*(const _b_int &a, const _b_int &b) { return _b_int(a) *= b; }
    friend _b_int operator/(const _b_int &a, const _b_int &b) { return _b_int(a) /= b; }
 
    _b_int& operator++() {
        val = val == MOD - 1 ? 0 : val + 1;
        return *this;
    }
 
    _b_int& operator--() {
        val = val == 0 ? MOD - 1 : val - 1;
        return *this;
    }
 
    _b_int operator++(int) { _b_int before = *this; ++*this; return before; }
    _b_int operator--(int) { _b_int before = *this; --*this; return before; }
 
    _b_int operator-() const {
        return val == 0 ? 0 : MOD - val;
    }
 
    friend bool operator==(const _b_int &a, const _b_int &b) { return a.val == b.val; }
    friend bool operator!=(const _b_int &a, const _b_int &b) { return a.val != b.val; }
    friend bool operator<(const _b_int &a, const _b_int &b) { return a.val < b.val; }
    friend bool operator>(const _b_int &a, const _b_int &b) { return a.val > b.val; }
    friend bool operator<=(const _b_int &a, const _b_int &b) { return a.val <= b.val; }
    friend bool operator>=(const _b_int &a, const _b_int &b) { return a.val >= b.val; }
 
    _b_int inv() const {
        return inv_mod(val);
    }
 
    _b_int pow(int64_t p) const {
        if (p < 0)
            return inv().pow(-p);
 
        _b_int a = *this, result = 1;
 
        while (p > 0) {
            if (p & 1)
                result *= a;
 
            p >>= 1;
 
            if (p > 0)
                a *= a;
        }
 
        return result;
    }
 
    friend ostream& operator<<(ostream &os, const _b_int &m) {
        return os << m.val;
    }
 
    friend istream& operator>>(istream &is, _b_int &m) {
        int64_t x;
        is >> x;
        m = x;
        return is;
    }
};
 
int MOD = 998244353;
barrett_reduction barrett(MOD);
using mnum = _b_int<MOD, barrett>;

void identityPerm(vector<int>& v, int n) {
  v.resize(n);
  iota(all(v), 0);
}
void compose(vector<int>& ret, const vector<int>& lhs, const vector<int>& rhs) {
  ret.resize(sz(lhs));
  for(int i = 0; i < sz(ret); i++) ret[i] = rhs[lhs[i]];
}
void invertPerm(vector<int>& ret, const vector<int>& v) {
  ret.resize(sz(v));
  for(int i = 0; i < sz(ret); i++) ret[v[i]] = i;
}

void genPerm(int n, vector<vector<vector<int>>>& ret, vector<vector<int>>& src) {
  ret.resize(n);
  vector<vector<int>> currperms = src;
  if(sz(currperms) == 0) {
    currperms.resize(1);
    identityPerm(currperms[0], n);
  }
  for(int i = n-1; i >= 0; i--) {
    vector<vector<int>> curr;
    vector<vector<int>> dp(n);
    vector<int> st;
    for(const auto& cand: currperms) {
      if(sz(dp[cand[i]]) == 0) {
        dp[cand[i]] = cand;
        st.pb(cand[i]);
      }
      while(sz(st)) {
        int ci = st.back(); st.pop_back();
        for(const auto& ncand: currperms) {
          int ni = ncand[ci];
          if(sz(dp[ni]) == 0) {
            compose(dp[ni], dp[ci], ncand);
            st.pb(ni);
          }
        }
      }
    }
    vector<vector<int>> perms;
    for(int j = 0; j < n; j++) {
      if(sz(dp[j]) > 0) perms.pb(dp[j]);
    }
    vector<vector<int>> nextperms;
    vector<vector<vector<int>>> cache(n);
    for(auto& x: cache) {
      x.resize(n);
    }
    for(int j = 0; j < n; j++) {
      if(sz(dp[j]) == 0) continue;
      for(const auto& p: currperms) {
        vector<int> cand;
        vector<int> rhs;
        invertPerm(rhs, dp[p[j]]);
        compose(cand, dp[j], p);
        vector<int> ncand;
        compose(ncand, cand, rhs);
        for(int k = i-1; k >= 0; k--) {
          if(sz(cache[k][ncand[k]]) > 0) {
            vector<int> ncandinv;
            invertPerm(ncandinv, ncand);
            compose(ncand, cache[k][ncand[k]], ncandinv);
          }
          else {
            cache[k][ncand[k]] = ncand;
            nextperms.pb(ncand);
            break;
          }
        }
      }
    }
    ret[i] = perms;
    swap(currperms, nextperms);
  }
}

void solve() {
  int c, r;
  cin >> c >> r;
  vector<vector<int>> g(r);
  for(auto& x: g) {
    x.resize(c);
    for(auto& y: x) {
      cin >> y; y--;
    }
  }
  vector<vector<vector<int>>> allv;
  vector<vector<int>> perms;
  for(const auto& p: g) {
    vector<int> cand(c);
    for(int i = 0; i < c; i++) cand[i] = p[i];
    perms.pb(cand);
  }
  genPerm(c, allv, perms);
  vector<vector<mnum>> dp(c);
  for(int i = 0; i < c; i++) {
    dp[i].resize(c);
    for(int j = i+1; j < c; j++) dp[i][j] = 1;
  }
  for(const auto& curr: allv) {
    vector<vector<mnum>> ndp(c);
    for(auto& x: ndp) x.resize(c);
    for(const auto& p: curr) {
      for(int a = 0; a < c; a++) for(int b = 0; b < c; b++) ndp[p[a]][p[b]] += dp[a][b];
    }
    swap(dp, ndp);
  }
  mnum ret = 0;
  for(int i = 0; i < c; i++) for(int j = 0; j < i; j++) ret += dp[i][j];
  cout << ret << "\n";
}

// what would chika do
// are there edge cases (N=1?)
// are array sizes proper (scaled by proper constant, for example 2* for koosaga tree)
// integer overflow?
// DS reset properly between test cases
// are you doing geometry in floating points
// are you not using modint when you should

int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  solve();
}