#include <queue>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/hash_policy.hpp>
using namespace std;
int read(){
	char c=getchar();int x=0;
	while(c<48||c>57) c=getchar();
	do x=(x<<1)+(x<<3)+(c^48),c=getchar();
	while(c>=48&&c<=57);
	return x;
}
const int N=400003,M=800003,P=1000000007;
typedef long long ll;
int n,m,cnt;
int dfn[N],low[N],num;
int stk[N],tp,d[N],deg[N];
int hd[N],ver[M],nxt[M],tot;
int eu[M],ev[M],fac[M];
inline void add(int u,int v){nxt[++tot]=hd[u];hd[u]=tot;ver[tot]=v;}
inline void chmn(int &x,int v){if(x>v) x=v;}
ll res;
__gnu_pbds::gp_hash_table<int,bool> mp[N];
vector<int> e[N],vec[N];
namespace tr{
	int hd[N],ver[N<<1],nxt[N<<1],tot;
	inline void add(int u,int v){
		vec[u-n].emplace_back(v);++deg[v];
		nxt[++tot]=hd[u];hd[u]=tot;ver[tot]=v;
		nxt[++tot]=hd[v];hd[v]=tot;ver[tot]=u;
	}
	int ft[N],de[N];
	void dfs(int u,int fa){
		ft[u]=fa;de[u]=de[fa]+1;
		for(int i=hd[u];i;i=nxt[i])
			if(ver[i]!=fa) dfs(ver[i],u);
	}
	inline int finds(int x,int y){
		if(ft[x]==ft[y]) return ft[x]-n;
		if(de[x]<de[y]) return ft[y]-n;
		else return ft[x]-n;
	}
}
void tarjan(int u){
	low[u]=dfn[u]=++num;stk[++tp]=u;
	for(int i=hd[u];i;i=nxt[i]){
		int v=ver[i];
		if(dfn[v]) chmn(low[u],dfn[v]);
		else{
			tarjan(v);
			chmn(low[u],low[v]);
			if(low[v]==dfn[u]){
				int x;tr::add(++cnt,u);
				do tr::add(cnt,x=stk[tp--]);while(x!=v);
			}
		}
	}
}
bool vis[N];
bool solve(vector<int> V,vector<int> E){
	for(int x:E) ++d[eu[x]],++d[ev[x]],mp[eu[x]][ev[x]]=mp[ev[x]][eu[x]]=1;
	queue<int> que;
	for(int x:V) if(d[x]==2) que.emplace(x);
	for(int i=1;i<=int(V.size()-2);++i){
		if(que.empty()) return 1;
		int u=que.front();
		que.pop();
		int x=0,y=0;
		vis[u]=1;
		for(auto [p,c]:mp[u])
			if(c){
				if(x) y=p;
				else x=p;
			}
		mp[x][u]=mp[y][u]=0;
		if(mp[x][y]){
			if(--d[x]==2) que.emplace(x);
			if(--d[y]==2) que.emplace(y);
		}
		else mp[x][y]=mp[y][x]=1;
	}
	int a=0,b=0;
	for(int x:V) if(!vis[x]){if(a) b=x;else a=x;}
	for(int x:V) d[x]=0,mp[x].clear();
	return 0;
}
int main(){
	n=read();m=read();cnt=n;fac[0]=1;
	for(int i=1;i<=m;++i){
		fac[i]=(ll)fac[i-1]*i%P;
		eu[i]=read();
		ev[i]=read();
		add(eu[i],ev[i]);add(ev[i],eu[i]);
	}
	tarjan(1);tr::dfs(1,0);
	for(int i=1;i<=m;++i) e[tr::finds(eu[i],ev[i])].emplace_back(i);
	int res=n;
	for(int i=1;i<=n;++i) res=(ll)res*fac[deg[i]]%P;
	for(int i=1;i<=cnt-n;++i){
		if(vec[i].size()<=2lu) continue;
		if(solve(vec[i],e[i])){puts("0");return 0;}
		res<<=1;if(res>=P) res-=P;
	}
	printf("%d\n",res);
	return 0;
}