//Author: Kevin
#include<bits/stdc++.h>
//#pragma GCC optimize("O2")
using namespace std;
#define ll long long
#define ull unsigned ll
#define pb emplace_back
#define mp make_pair
#define ALL(x) (x).begin(),(x).end()
#define rALL(x) (x).rbegin(),(x).rend()
#define srt(x) sort(ALL(x))
#define rev(x) reverse(ALL(x))
#define rsrt(x) sort(rALL(x))
#define sz(x) (int)(x.size())
#define inf 0x3f3f3f3f
#define pii pair<int,int>
#define lb(v,x) (int)(lower_bound(ALL(v),x)-v.begin())
#define ub(v,x) (int)(upper_bound(ALL(v),x)-v.begin())
#define uni(v) v.resize(unique(ALL(v))-v.begin())
#define longer __int128_t
void die(string S){puts(S.c_str());exit(0);}
const ll mod=998244353;
mt19937_64 rng(20210448);
using perm=vector<int>;
int n,m;
perm I;
perm inv(const perm &p)
{
	vector<int> vec(n);
	for(int i=0;i<n;i++)
		vec[p[i]]=i;
	return vec;
}
perm comp(const perm &p,const perm &q)
{
	vector<int> vec(n);
	for(int i=0;i<n;i++)
		vec[i]=p[q[i]];
	return vec;
}
perm P[33];
vector<perm> vect[33];
vector<perm> E;
vector<pii> G[33];
int flag[50500];
int vis[33];
perm ind[33];
int u[50500],v[50500];
void dfs(int u)
{
	vis[u]=1;
	for(auto pr:G[u])
		if(!vis[pr.first])
		{
			flag[pr.second]=flag[pr.second^1]=1;
			ind[pr.first]=comp(E[pr.second],ind[u]);
			dfs(pr.first);
		}
}
ll f[33][33];
int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m;
	for(int i=0;i<n;i++)
		I.pb(i);
	for(int i=1;i<=m;i++)
	{
		for(int j=0;j<n;j++)
		{
			int x;
			cin>>x;
			x--;
			P[i].pb(x);
		}
	}
	vector<perm> pool;
	for(int i=1;i<=m;i++)
		pool.pb(P[i]);
	for(int i=0;i<n;i++)
	{
		for(int j=i;j<n;j++)
			G[j].clear();
		int tot=0;
		E.clear();
		for(auto p:pool)
			for(int j=i;j<n;j++)
			{
				G[p[j]].pb(j,tot);
				u[tot]=p[j];
				v[tot]=j;
				E.pb(inv(p));
				tot++;
				G[j].pb(p[j],tot);
				E.pb(p);
				tot++;
			}
		memset(vis,0,sizeof(vis));
		memset(flag,0,sizeof(flag));
		ind[i]=I;
		dfs(i);
		for(int j=i;j<n;j++)
			if(vis[j])
				vect[i].pb(ind[j]);
		vector<perm> npool;
		for(int j=0;j<tot;j+=2)
			if(!flag[j]&&vis[u[j]])
				npool.pb(comp(comp(inv(ind[v[j]]),E[j]),ind[u[j]]));
		pool.clear();
		for(int j=0;j<n+n;j++)
		{
			perm p=I;
			for(int k=0;k<=sz(npool)*2;k++)
				p=comp(p,npool[rng()%sz(npool)]);
			pool.pb(p);
		}
	}
	for(int i=0;i<n;i++)
		for(int j=i+1;j<n;j++)
			f[i][j]=1;
	for(int i=n-1;i>=0;i--)
	{
		ll g[33][33];
		memset(g,0,sizeof(g));
		for(int a=0;a<n;a++)
			for(int b=0;b<n;b++)
				for(auto p:vect[i])
					g[p[a]][p[b]]=(g[p[a]][p[b]]+f[a][b])%mod;
		memcpy(f,g,sizeof(f));
	}
	ll ans=0;
	for(int i=0;i<n;i++)
		for(int j=0;j<i;j++)
			ans=(ans+f[i][j])%mod;
	cout<<ans<<endl;
	return 0;
}