/* -*- coding: utf-8 -*- * * 2788.cc: No.2788 4-33 Hard - yukicoder */ #include #include using namespace std; /* constant */ const int SA = 4; const int SB = 33; const int MOD = 998244353; /* typedef */ using ll = long long; template struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } 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); } }; using mi = MI; /* global variables */ ll qs[SA + 1][SB + 1], xs[SA + 1][SB + 1]; mi dp[2][9][SA + 1][SB + 1]; mi invs[9], cs[9]; /* subroutines */ /* main */ int main() { for (int i = 1; i < 9; i++) invs[i] = mi(i).inv(); for (int i = 0; i <= SA; i++) for (int j = 0; j <= SB; j++) scanf("%lld", qs[i] + j); for (int i = 0; i <= SA; i++) for (int j = 0; j <= SB; j++) scanf("%lld", xs[i] + j); dp[0][0][0][0] = 1; int cur = 0, nxt = 1; for (int a = 0; a <= SA; a++) for (int b = 0; b <= SB; b++) { ll qab = qs[a][b]; cs[0] = 1; for (int i = 1; i < 9; i++) cs[i] = cs[i - 1] * (qab + 1 - i) * invs[i]; for (int j = 0; j < 9; j++) for (int ka = 0; ka <= SA; ka++) copy(dp[cur][j][ka], dp[cur][j][ka] + SB + 1, dp[nxt][j][ka]); for (int i = 1; i < 9; i++) for (int j = 0; j + i < 9; j++) for (int ka = 0; ka + a * i <= SA; ka++) for (int kb = 0; kb + b * i <= SB; kb++) dp[nxt][j + i][ka + a * i][kb + b * i] += dp[cur][j][ka][kb] * cs[i]; swap(cur, nxt); } mi sum = 0; for (int a = 0; a <= SA; a++) sum += dp[cur][8][SA - a][SB] * xs[a][0]; printf("%d\n", (int)sum); return 0; }