#include <bits/stdc++.h>

using namespace std;


#define ll long long
#define flt double
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define rep(i, n) for(int i = 0; i < n; i += 1)


const ll inf = 1e18;
mt19937 rnd(234);
const ll mod = 998244353;

const int maxn = 1000;
int n, m;
char g[maxn][maxn];


void swap_vertices(int u, int v) {
    if (u == v) {
        return;
    }
    rep(i, n) {
        swap(g[i][u], g[i][v]);
    }
    rep(i, n) {
        swap(g[u][i], g[v][i]);
    }
}


void make_flex(int u, int v) {
    assert(u != v);
    rep(i, n) g[v][i] ^= g[u][i];
    rep(i, n) g[i][v] ^= g[i][u];
}


int32_t main() {
    if (1) {
        ios::sync_with_stdio(0);
        cin.tie(0);
        cout.tie(0);
    }
    cin >> n >> m;
    rep(i, n) rep(j, n) g[i][j] = 0;
    for (int i = 0; i < m; i += 1) {
        int u, v;
        cin >> u >> v;
        --u;
        --v;
        g[u][v] = 1;
    }
    for (int i = 0; i < n - 1; i += 1) {
        int to = -1;
        for (int j = i + 1; j < n; j += 1) {
            if (g[i][j] ^ g[j][i]) {
                to = j;
                break;
            }
        }
        if (to == -1) {
            continue;
        }
        for (int j = to + 1; j < n; j += 1) {
            if (g[i][j] ^ g[j][i]) {
                make_flex(to, j);
                assert(g[i][j] ^ g[j][i] == 0);
            }
        }
        swap_vertices(to, i + 1);
    }
    array<array<ll, 2>, 2> dp;
    rep(i, 2) rep(j, 2) dp[i][j] = 0;
    dp[0][0] = 1;
    for (int i = 0; i < n; i += 1) {
        int d = 0;
        if (i > 0) {
            d = g[i][i - 1] ^ g[i - 1][i];
        }
        array<array<ll, 2>, 2> ndp;
        rep(x, 2) rep(y, 2) ndp[x][y] = 0;
        rep(prev, 2) rep(prev_sum, 2) rep(cur, 2) {
            auto& new_val = ndp[(prev_sum + cur * g[i][i] + d * cur * prev) % 2][cur];
            new_val += dp[prev_sum][prev];
            if (new_val >= mod) {
                new_val -= mod;
            }
        }
        dp.swap(ndp);
    }
    ll result = (dp[0][0] + dp[0][1]) % mod;
    cout << result << "\n";
    return 0;
}

/*
5 5
1 2
2 3
3 4
4 5
1 5

*/