#include<bits/stdc++.h>
using namespace std;
const int mod=998244353;
int n,m,jc[1010000];
#define pb push_back
vector<int>g[1010000],G[1010000],E[1010000];
set<int>s[1010000];
int dfn[1010000],low[1010000],sta[1010000],top,ru[1010000],an,ff[1010000],N;
void dfs(int x){
	dfn[x]=low[x]=++dfn[0];
	sta[++top]=x;
	for(int v:g[x]){
		if(!dfn[v]){
			dfs(v),low[x]=min(low[x],low[v]);
			if(low[v]>=dfn[x]){
				ff[v]=++N,G[N].pb(v);
				if(sta[top]!=v)an=2*an%mod;
				while(sta[top]!=v)G[N].pb(sta[top]),ru[sta[top]]++,ff[sta[top]]=N,top--;
				G[N].pb(x);
				ff[N]=x,ru[v]++,ru[x]++,top--;
			}
		}
		low[x]=min(low[x],dfn[v]);
	}
}
int U[1010000],V[1010000],yt[1010000];
int main(){
	scanf("%d%d",&n,&m);
	N=n;
	jc[0]=1;
	for(int i=1;i<=n;i++)jc[i]=1ll*i*jc[i-1]%mod;
	for(int i=1;i<=m;i++)scanf("%d%d",&U[i],&V[i]),g[U[i]].pb(V[i]),g[V[i]].pb(U[i]);
	an=n;dfs(1);
	for(int i=1;i<=n;i++)an=1ll*an*jc[ru[i]]%mod;
	memset(ru,0,sizeof(ru));
	for(int i=1;i<=m;i++){
		int u=U[i],v=V[i];
		if(ff[u]==ff[v])E[ff[u]].pb(i);
		else if(ff[ff[u]]==v)E[ff[u]].pb(i);
		else E[ff[v]].pb(i);
	}
	for(int i=n+1;i<=N;i++)if((int)G[i].size()>2){
		for(int e:E[i])s[U[e]].insert(V[e]),s[V[e]].insert(U[e]);
		queue<int>q;for(int x:G[i]){
			if(s[x].size()==2)q.push(x);
		}
		int step=(int)G[i].size();step-=2;
		while(!q.empty()&&step){
			step--;
			int u=q.front();q.pop();
			auto i1=s[u].begin(),i2=--s[u].end();
			int l=*i1,r=*i2;
			yt[u]+=2,yt[l]++,yt[r]++;
			s[l].erase(u);
			s[r].erase(u);
			if(s[l].find(r)==s[l].end())s[l].insert(r),s[r].insert(l);
			if(s[l].size()==2)q.push(l);
			if(s[r].size()==2)q.push(r);
		}
		if(step){puts("0");return 0;}
		for(int x:G[i]){
			s[x].clear();
			if(yt[x]!=2){puts("0");return 0;}
			yt[x]=0;
		}
	}
	printf("%d",an);
	return 0;
}