#include #include #include #include #include #include #include #include #include using mint = atcoder::static_modint<943718401>; mint S(const int n, const mint& x, std::vector a, std::vector b, std::vector c) { assert(n == std::ssize(a) - 1 && n == std::ssize(b) - 1 && n == std::ssize(c) - 1); // ijk \neq 0 のとき A_i B_j C_k x^{ijk} = 0 が成り立つ。 // ijk = 0 \iff j = 0 または k = 0 に注意する。 if (x == 0) { mint s = 0; // j = 0 のとき s += std::inner_product(a.begin(), a.end(), b.begin(), mint::raw(0)) * c[0]; // j > 0 かつ k = 0 のとき s += std::inner_product(std::next(a.begin()), a.end(), std::next(c.begin()), mint::raw(0)) * b[0]; return s; } const auto T = [](const int i) -> std::int64_t { return std::int64_t{i - 1} * i * (i + 1) / 3; }; // A_i x^{t_i} std::ranges::for_each( a, [&x, T, i = 0](mint& e) mutable -> void { e *= x.pow(T(i++)); }); // B_j / x^{t_j}, C_k / x^{t_k} const mint inv = x.inv(); std::ranges::for_each( b, [&inv, T, i = 0](mint& e) mutable -> void { e *= inv.pow(T(i++)); }); std::ranges::for_each( c, [&inv, T, i = 0](mint& e) mutable -> void { e *= inv.pow(T(i++)); }); const std::vector s = atcoder::convolution(b, c); return std::inner_product(a.begin(), a.end(), s.begin(), mint::raw(0)); } int main() { constexpr int kMaxN = 200000; const auto GetMint = []() -> mint { int x; std::cin >> x; assert(0 <= x && x < mint::mod()); return mint::raw(x); }; int n; std::cin >> n; assert(0 <= n && n <= kMaxN); const mint x = GetMint(); std::vector a(n + 1), b(n + 1), c(n + 1); std::ranges::for_each(a, [GetMint](mint& x) -> void { x = GetMint(); }); std::ranges::for_each(b, [GetMint](mint& x) -> void { x = GetMint(); }); std::ranges::for_each(c, [GetMint](mint& x) -> void { x = GetMint(); }); std::cout << S(n, x, a, b, c).val() << '\n'; return 0; }