/*

                                _/                                      _/
       _/_/        _/_/_/      _/_/_/        _/_/_/      _/  _/_/               _/_/_/
    _/    _/    _/            _/    _/    _/    _/      _/_/          _/       _/    _/
   _/    _/    _/            _/    _/    _/    _/      _/            _/       _/    _/
    _/_/        _/_/_/      _/    _/      _/_/_/      _/            _/       _/    _/

*/
#include<bits/stdc++.h>
#define int long long

using namespace std;
const int mod=998244353;

void solve(){
    int n,m;
    cin>>n>>m;
    vector<int>p(n+1),pp(n+1);
    auto cal=[&](int x)->int{
        return p[x%n]*pp[x/n]%mod;
    };
    p[0]=pp[0]=1;
    for(int i=1;i<=n;++i) p[i]=p[i-1]*2%mod;
    for(int i=1;i<=n;++i) pp[i]=pp[i]*p[n]%mod;
    vector<int>head(n,-1),to(2*m),nxt(2*m,-1);
    int cnt=-1;
    auto add=[&](int u,int v)->void{
        to[++cnt]=v;
        nxt[cnt]=head[u];
        head[u]=cnt;
    };
    for(int i=0;i<m;++i){
        int u,v;
        cin>>u>>v;
        --u,--v;
        add(u,v);
        add(v,u);
    }
    vector<int>dfn(n),low(n);
    int idx=0;
    stack<int>st;
    vector<int>id(n);
    auto tarjan=[&](auto tarjan,int u,int in)->void{
        dfn[u]=low[u]=++cnt;
        st.push(u);
        for(int i=head[u];~i;i=nxt[i]){
            int v=to[i];
            if(!dfn[v]){
                tarjan(tarjan,v,i^1);
                low[u]=min(low[u],low[v]);
            }
            else if(i!=in) low[u]=min(low[u],dfn[v]);
        }
        if(dfn[u]==low[u]){
            int x=st.top();st.pop();
            id[x]=idx;
            while(x!=u){
                x=st.top();st.pop();
                id[x]=idx;
            }
            idx++;
        }
    };
    tarjan(tarjan,0,-1);
    vector<int>sz(idx),c(idx);
    vector<vector<int>>e(idx);
    for(int i=0;i<n;++i){
        int u=id[i];
        sz[u]++;
        for(int j=head[i];~j;j=nxt[j]){
            int v=id[to[j]];
            if(u!=v) e[u].push_back(v);
            else c[u]++;
        }
    }
    vector f(idx,vector(n+1,0ll));
    auto dfs=[&](auto dfs,int u,int fa)->void{
        f[u][sz[u]]=cal((sz[u]-1)*sz[u]/2-c[u]/2);
        for(auto v:e[u]){
            if(v==fa) continue;
            dfs(dfs,v,u);
            vector<int>g(sz[u]+sz[v]+1,0ll);
            for(int i=0;i<=sz[u];++i){
                for(int j=0;j<=sz[v];++j){
                    g[i+j]+=f[u][i]*f[v][j]%mod*cal(i*j-1)%mod;
                    if(g[i+j]>=mod) g[i+j]-=mod;
                    g[i]+=mod-f[u][i]*f[v][j]%mod;
                    if(g[i]>=mod) g[i]-=mod;
                }
            }
            swap(f[u],g);
            sz[u]+=sz[v];
        }
    };
    dfs(dfs,0,-1);
    int res=0;
    for(int i=0;i<=n;++i){
        res+=f[0][i];
        // cerr<<f[0][i]<<endl;
        if(res>=mod) res-=mod;
    }
    cout<<res<<"\n";
}

signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    solve();
}