/* -*- coding: utf-8 -*- * * 2428.cc: No.2428 Returning Shuffle - yukicoder */ #include #include #include #include using namespace std; /* constant */ const int MAX_N = 1000000; const int MAX_P = 1000000; const int MOD = 998244353; /* typedef */ typedef vector vi; typedef pair pii; typedef vector vpii; template struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) {} MI(long long _v): v(_v % MOD) {} MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; typedef MI mi; /* global variables */ bool primes[MAX_P + 1]; int ps[MAX_N], ss[MAX_N], qs[MAX_N]; bool used[MAX_N]; /* subroutines */ int gen_primes(int maxp, vi &pnums) { fill(primes, primes + maxp + 1, true); primes[0] = primes[1] = false; int p; for (p = 2; p * p <= maxp; p++) if (primes[p]) { pnums.push_back(p); for (int q = p * p; q <= maxp; q += p) primes[q] = false; } for (; p <= maxp; p++) if (primes[p]) pnums.push_back(p); return (int)pnums.size(); } bool prime_decomp(int n, vi &pnums, vpii& pds) { pds.clear(); int pn = pnums.size(); for (int i = 0; i < pn; i++) { int pi = pnums[i]; if (pi * pi > n) { if (n > 1) pds.push_back(pii(n, 1)); return true; } if (n % pi == 0) { int fi = 0; while (n % pi == 0) n /= pi, fi++; pds.push_back(pii(pi, fi)); } } return false; } /* main */ int main() { vi pnums; gen_primes(MAX_P, pnums); int n, m; scanf("%d%d", &n, &m); for (int i = 0; i < n; i++) ps[i] = i; for (int i = 0; i < m; i++) { int ti; scanf("%d", &ti); for (int j = 0; j < ti; j++) scanf("%d", ss + j), ss[j]--; for (int j = 0; j < ti; j++) qs[j] = ps[ss[j]]; for (int j = 0; j < ti; j++) ps[ss[(j + 1) % ti]] = qs[j]; } //for (int i = 0; i < n; i++) printf("%d ", ps[i]); putchar('\n'); vpii lpds; for (int st = 0; st < n; st++) if (! used[st]) { int l = 0; for (int u = st; ! used[u]; u = ps[u]) used[u] = true, l++; if (l > 1) { vpii pds, lpds1; prime_decomp(l, pnums, pds); int i = 0, j = 0, n0 = lpds.size(), n1 = pds.size(); while (i < n0 && j < n1) { if (lpds[i].first < pds[j].first) lpds1.push_back(lpds[i++]); else if (lpds[i].first > pds[j].first) lpds1.push_back(pds[j++]); else { int e = max(lpds[i].second, pds[j].second); lpds1.push_back(pii(lpds[i].first, e)); i++, j++; } } while (i < n0) lpds1.push_back(lpds[i++]); while (j < n1) lpds1.push_back(pds[j++]); swap(lpds, lpds1); } } mi p = 1; for (auto &pd: lpds) p *= mi(pd.first).pow(pd.second); printf("%d\n", p.v); return 0; }