#include<bits/stdc++.h>
#ifndef xxx
#define dbg(...) ;
#endif
using namespace std;
template<class T> ostream &operator<<(ostream &o, const vector <T> &v) { bool f=false; for(T i:v) { f?o<<' ':o; o<<(i); f=true; } return o; }
template<class T> void sort(T &v) { std::sort(v.begin(), v.end()); }
template<class T, class C> void sort(T &v, C c) { std::sort(v.begin(), v.end(), c); }
template<class T> void reverse(T &v) { std::reverse(v.begin(), v.end()); }
template<class T> bool is_sorted(T &v) { return std::is_sorted(v.begin(), v.end()); }
template<class T> inline void _min(T &x, const T &y) { x>y?x=y:x; }
template<class T> inline void _max(T &x, const T &y) { x<y?x=y:x; }
istream &operator>>(istream &i, __int128_t &x) { x=0; short f=1; char c=0; while(!isdigit(c)) { if(c=='-')f=-1; c=i.get(); } while(isdigit(c))x=x*10+c-'0', c=i.get(); x*=f; return i; }
ostream &operator<<(ostream &o, __int128_t x) { if(x==0) { o<<0; return o; } bool f=false; string s; if(x<0)f=true, x=-x; while(x)s+=x%10+'0', x/=10; reverse(s); if(f)o<<'-'; o<<s; return o; }
mt19937 mt(time(0));
typedef double db;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
constexpr int maxn=180+5;
constexpr int mod=1e9+7;
constexpr int inf=0x3f3f3f3f;
constexpr ll INF=0x3f3f3f3f3f3f3f3fll;
constexpr double pi=acos(-1.0);
constexpr double eps=1e-9;
vector<vector<int>>e;
vector<int>composition(const vector<int> &a, const vector<int> &b)
{
	if(!a.size())return b;
	vector<int>ans;
	for(int i=1; i<=a.size(); i++)
	{
		ans.push_back(a[b[i-1]-1]);
	}
	return ans;
}
int main()//MARK: main
{
#ifndef xxx
	ios::sync_with_stdio(false), cin.tie(nullptr), cout.tie(nullptr);
#endif
	cout<<fixed<<setprecision(10);
	int n, m;
	cin>>n>>m;
	for(int i=1; i<=m; i++)
	{
		vector<int>vec(n);
		for(int &x:vec)
		{
			cin>>x;
		}
		e.push_back(vec);
	}
	function<void(int, int, vector<int> &, map<vector<int>, bool> &)>dfs=[&](int step, int lim, vector<int> &now, map<vector<int>, bool> &mp)
	{
		if(step==lim+1)
		{
			mp[now]=true;
			return;
		}
		vector<int> nex=composition(now, e[step]);
		dfs(step+1, lim, nex, mp);
		dfs(step+1, lim, now, mp);
	};
	auto getans=[&]()->int
	{
		int ans=0;
		if(n==1)
		{
			return 1;
		}
		else if(n==2)
		{
			const vector<int>v1={1, 2};
			const vector<int>v2={2, 1};
			if(find(e.begin(), e.end(), v1)!=e.end()&&find(e.begin(), e.end(), v2)!=e.end())
			{
				return 2;
			}
			else if(count(e.begin(), e.end(), v2)>=2)
			{
				return 2;
			}
			return 1;
		}
		if(e.size()<=20)
		{
			map<vector<int>, bool>mp;
			vector<int>now;
			dfs(0, e.size()-1, now, mp);
			ans=mp.size()-1;
		}
		else if(e.size()<=36)
		{
			map<vector<int>, bool>mp1, mp2, mp;
			vector<int>now1, now2;
			int mid=e.size()/2;
			dfs(0, mid, now1, mp1);
			dfs(mid+1, e.size()-1, now2, mp2);
			mp=mp1;
			for(auto &p:mp2)
				mp.insert(p);
			for(const auto &[a, f1]:mp1)
			{
				for(const auto &[b, f2]:mp2)
				{
					mp[composition(a, b)]=true;
				}
			}
			ans=mp.size()-1;
		}
		else
		{
			map<vector<int>, bool>mp1, mp2, mp3, mp4, mp;
			vector<int>now1, now2, now3, now4;
			int l1=e.size()/4;
			int l2=2*l1;
			int l3=3*l1;
			dfs(0, l1, now1, mp1);
			dfs(l1+1, l2, now2, mp2);
			dfs(l2+1, l3, now3, mp3);
			dfs(l3+1, e.size()-1, now4, mp4);
			mp=mp1;
			for(auto &p:mp2)
				mp.insert(p);
			for(auto &p:mp3)
				mp.insert(p);
			for(auto &p:mp4)
				mp.insert(p);
			for(const auto &[a, f1]:mp1)
			{
				for(const auto &[b, f2]:mp2)
					for(const auto &[c, f3]:mp3)
						for(const auto &[d, f4]:mp4)
						{
							auto comp1=composition(a, b);
							auto comp2=composition(comp1, c);
							auto comp3=composition(comp2, d);
							mp[comp3]=true;
						}
			}
			ans=mp.size()-1;
		}
		return ans;
	};
	cout<<getans()<<'\n';
	return 0;
}