//#define NDEBUG #include #include #include #include #include #include namespace n91 { using i8 = std::int_fast8_t; using i32 = std::int_fast32_t; using i64 = std::int_fast64_t; using u8 = std::uint_fast8_t; using u32 = std::uint_fast32_t; using u64 = std::uint_fast64_t; using isize = std::ptrdiff_t; using usize = std::size_t; struct rep { struct itr { usize i; constexpr itr(const usize i) noexcept : i(i) {} void operator++() noexcept { ++i; } constexpr usize operator*() const noexcept { return i; } constexpr bool operator!=(const itr x) const noexcept { return i != x.i; } }; const itr f, l; constexpr rep(const usize f, const usize l) noexcept : f(std::min(f, l)), l(l) {} constexpr auto begin() const noexcept { return f; } constexpr auto end() const noexcept { return l; } }; struct revrep { struct itr { usize i; constexpr itr(const usize i) noexcept : i(i) {} void operator++() noexcept { --i; } constexpr usize operator*() const noexcept { return i; } constexpr bool operator!=(const itr x) const noexcept { return i != x.i; } }; const itr f, l; constexpr revrep(const usize f, const usize l) noexcept : f(l - 1), l(std::min(f, l) - 1) {} constexpr auto begin() const noexcept { return f; } constexpr auto end() const noexcept { return l; } }; template auto md_vec(const usize n, const T &value) { return std::vector(n, value); } template auto md_vec(const usize n, Args... args) { return std::vector(n, md_vec(args...)); } template constexpr T difference(const T &a, const T &b) noexcept { return a < b ? b - a : a - b; } template void chmin(T &a, const T &b) noexcept { if (b < a) a = b; } template void chmax(T &a, const T &b) noexcept { if (a < b) a = b; } template class fix_point : private F { public: explicit constexpr fix_point(F &&f) : F(std::forward(f)) {} template constexpr decltype(auto) operator()(Args &&... args) const { return F::operator()(*this, std::forward(args)...); } }; template constexpr decltype(auto) make_fix(F &&f) { return fix_point(std::forward(f)); } template T scan() { T ret; std::cin >> ret; return ret; } } // namespace n91 #include namespace n91 { template class modint { using u64 = std::uint_fast64_t; public: using value_type = u64; static constexpr u64 mod = Modulus; private: static_assert(mod < static_cast(1) << 32, "Modulus must be less than 2**32"); u64 v; constexpr modint &negate() noexcept { if (v != 0) v = mod - v; return *this; } public: constexpr modint(const u64 x = 0) noexcept : v(x % mod) {} constexpr u64 &value() noexcept { return v; } constexpr const u64 &value() const noexcept { return v; } constexpr modint operator+() const noexcept { return modint(*this); } constexpr modint operator-() const noexcept { return modint(*this).negate(); } constexpr modint operator+(const modint rhs) const noexcept { return modint(*this) += rhs; } constexpr modint operator-(const modint rhs) const noexcept { return modint(*this) -= rhs; } constexpr modint operator*(const modint rhs) const noexcept { return modint(*this) *= rhs; } constexpr modint operator/(const modint rhs) const noexcept { return modint(*this) /= rhs; } constexpr modint &operator+=(const modint rhs) noexcept { v += rhs.v; if (v >= mod) v -= mod; return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (v < rhs.v) v += mod; v -= rhs.v; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { v = v * rhs.v % mod; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = mod - 2; while (exp) { if (exp % 2 != 0) *this *= rhs; rhs *= rhs; exp /= 2; } return *this; } constexpr bool operator==(const modint rhs) const noexcept { return v == rhs.v; } constexpr bool operator!=(const modint rhs) const noexcept { return v != rhs.v; } }; template constexpr typename modint::u64 modint::mod; } // namespace n91 #include #include namespace n91 { template > constexpr T power(T base, U exp, const Operate &oper = Operate(), T iden = 1) { while (exp != 0) { if (exp % 2 != 0) { iden = oper(iden, base); } exp /= 2; base = oper(base, base); } return iden; } } // namespace n91 using namespace std; template void UFWT(vector &a, T rev2) { int n = a.size(); for (int d = 1; d < n; d <<= 1) { for (int m = d << 1, i = 0; i < n; i += m) { for (int j = 0; j < d; j++) { T x = a[i + j], y = a[i + j + d]; a[i + j] = (x + y) * rev2; a[i + j + d] = (x - y) * rev2; } } } } namespace n91 { void main_() { /* std::ios::sync_with_stdio(false); std::cin.tie(nullptr); //*/ constexpr i64 MOD = 998244353; const auto mod = [&](const i64 x) { return (x % MOD + MOD) % MOD; }; using mint = modint<998244353>; const usize n = scan(); std::vector c(1 << n); for (auto &e : c) std::cin >> e; UFWT(c, 1); mint ans = 0; for (const usize i : rep(0, 1 << n)) { const i64 l = scan(); const i64 h = scan(); if (c[i] < 0) { ans += mod(c[i]) * mod(l); } else { ans += mod(c[i]) * mod(h); } } for (const usize i : rep(0, n)) { ans /= 2; } std::cout << ans.value() << std::endl; } } // namespace n91 int main() { n91::main_(); return 0; }