#include #include #include #include int main() { using namespace std; using namespace atcoder; int n; cin >> n; vector has(n, vector(10, 0)); vector cmp(n, vector()); for (int i = 0; i < n; i++) { for (int j = 0; j < 6; j++) { int a; cin >> a; has[i][a] = 1; } for (int j = 0; j < 10; j++) { if (has[i][j]) { cmp[i].push_back(j); } } } using mint = modint998244353; // mf_graph g(n + 11); // for (int i = 1; i < 10; i++) { // g.add_edge(i, 0, 0); // } // for (int i = 1; i < n; i++) { // g.add_edge(10, 11 + i, 1); // for (int j = 1; j < 10; j++) { // if (has[i][j]) { // g.add_edge(11 + i, j, 1); // } // } // } vector app(10); for (int i = 0; i < n; ++i) { for (int j = 0; j < 10; j++) { app[j] += has[i][j]; } } int nn = 100; vector fact(nn, 1), invfact(nn, 1); for (int i = 0; i < nn - 1; ++i) { fact[i + 1] = fact[i] * (i + 1); } invfact[nn - 1] = fact[nn - 1].inv(); for (int i = nn - 2; i >= 0; --i) { invfact[i] = invfact[i + 1] * (i + 1); } mint ans = 0; vector count(10); // auto print_vec = [](vector v) { // int l = v.size(); // for (int i = 0; i < l; i++) { // std::cout << v[i] << (i < l - 1 ? ' ' : '\n'); // } // }; auto dfs = [&](auto f, int d, int rem, mint invs) -> void { if (d == 10) { if (rem) { return; } mf_graph g(n + 11); for (int i = 1; i < 10; i++) { // g.add_edge(i, 0, 0); if (count[i]) { g.add_edge(i, 0, count[i]); } } for (int i = 0; i < n; i++) { g.add_edge(10, 11 + i, 1); for (int j: cmp[i]) { g.add_edge(11 + i, j, 1); } } int flow = g.flow(10, 0); // cout << "\t\t" << flow << endl; if (flow < n) { return; } ans += invs; return; } for (int u = 0; u <= rem && u <= app[d]; ++u) { // g.change_edge(d - 1, u, 0); count[d] = u; f(f, d + 1, rem - u, invs * invfact[u]); } }; dfs(dfs, 1, n, fact[n]); std::cout << ans.val() << endl; }