// BEGIN: ../contest594/H/main.cpp #line 1 "..::contest594::H::main.cpp" // BEGIN: pch.hpp #line 3 "pch.hpp" #if defined(__GNUC__) && !defined(__clang__) #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #endif #define dump(...) #define CPP_DUMP_SET_OPTION(...) #define CPP_DUMP_DEFINE_EXPORT_OBJECT(...) #define CPP_DUMP_DEFINE_EXPORT_ENUM(...) #define CPP_DUMP_DEFINE_DANGEROUS_EXPORT_OBJECT(...) // BEGIN: template.hpp #line 3 "template.hpp" #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include // BEGIN: utilities/fast_io.hpp #line 3 "utilities::fast_io.hpp" #include #include #include #include #include #include #include #include #include #include #include #include namespace m1une { namespace utilities { namespace internal { // Detect std::begin(x), std::end(x). template struct is_range : std::false_type {}; template struct is_range())), decltype(std::end(std::declval())) >> : std::true_type {}; template inline constexpr bool is_range_v = is_range::value; template using range_reference_t = decltype(*std::begin(std::declval())); template using range_value_t = std::remove_cv_t>>; template struct range_stored_value { using type = range_value_t; }; template struct range_stored_value>::value_type>> { using type = typename std::remove_cv_t>::value_type; }; template using range_stored_value_t = typename range_stored_value::type; // Treat strings and C strings as scalar output objects, not as ranges. template struct is_char_array : std::false_type {}; template struct is_char_array : std::bool_constant, char>> {}; template struct is_string_like : std::bool_constant< std::is_same_v, std::string> || std::is_same_v, const char*> || std::is_same_v, char*> || is_char_array>::value > {}; template inline constexpr bool is_string_like_v = is_string_like::value; // ModInt-like type: x.val() is printable, and x can be assigned from long long. template struct has_val_method : std::false_type {}; template struct has_val_method().val())>> : std::true_type {}; template inline constexpr bool has_val_method_v = has_val_method::value; template struct has_static_mod_raw : std::false_type {}; template struct has_static_mod_raw< T, std::void_t()))>> : std::true_type {}; template inline constexpr bool has_static_mod_raw_v = has_static_mod_raw::value; // libstdc++ before GCC 16 does not classify __int128 as an integral type in // strict ISO modes such as -std=c++23. Keep the fast-I/O interface independent // of that implementation detail. template inline constexpr bool is_integral_v = std::is_integral_v || std::is_same_v, __int128_t> || std::is_same_v, __uint128_t>; template inline constexpr bool is_signed_v = std::is_signed_v || std::is_same_v, __int128_t>; template struct make_unsigned { using type = std::make_unsigned_t; }; template <> struct make_unsigned<__int128_t> { using type = __uint128_t; }; template <> struct make_unsigned<__uint128_t> { using type = __uint128_t; }; template using make_unsigned_t = typename make_unsigned>::type; } // namespace internal struct FastInput { static constexpr int buffer_size = 1 << 20; private: std::FILE* _stream; char _buffer[buffer_size]; int _position; int _length; bool _terminal; bool refill() { _position = 0; if (_terminal) { if (std::fgets(_buffer, buffer_size, _stream) == nullptr) { _length = 0; return false; } _length = int(std::strlen(_buffer)); } else { _length = int(std::fread(_buffer, 1, buffer_size, _stream)); } return _length != 0; } template bool read_integer_from_terminal(T& value) { if (!skip_spaces()) return false; int c = read_char_raw(); bool negative = false; if (c == '-') { negative = true; c = read_char_raw(); } if constexpr (internal::is_signed_v) { T result = 0; while ('0' <= c && c <= '9') { result = negative ? result * 10 - (c - '0') : result * 10 + (c - '0'); c = read_char_raw(); } value = result; } else { T result = 0; while ('0' <= c && c <= '9') { result = result * 10 + T(c - '0'); c = read_char_raw(); } value = negative ? T(0) - result : result; } return true; } bool prepare_number() { if (_length - _position >= 64) return true; const int remaining = _length - _position; if (remaining > 0) std::memmove(_buffer, _buffer + _position, remaining); const int added = int(std::fread(_buffer + remaining, 1, buffer_size - remaining, _stream)); _position = 0; _length = remaining + added; if (_length < buffer_size) _buffer[_length] = '\0'; return _length != 0; } public: explicit FastInput(std::FILE* stream = stdin) : _stream(stream), _position(0), _length(0), _terminal(::isatty(::fileno(stream)) != 0) {} FastInput(const FastInput&) = delete; FastInput& operator=(const FastInput&) = delete; int read_char_raw() { if (_position == _length && !refill()) return EOF; return _buffer[_position++]; } bool skip_spaces() { int c = read_char_raw(); while (c != EOF && c <= ' ') c = read_char_raw(); if (c == EOF) return false; --_position; return true; } bool read(char& value) { if (!skip_spaces()) return false; value = char(read_char_raw()); return true; } bool read(std::string& value) { if (!skip_spaces()) return false; value.clear(); int c = read_char_raw(); while (c != EOF && c > ' ') { value.push_back(char(c)); c = read_char_raw(); } return true; } bool read(bool& value) { int x; if (!read(x)) return false; value = x != 0; return true; } template std::enable_if_t< internal::is_integral_v && !std::is_same_v, bool> && !std::is_same_v, char>, bool > read(T& value) { if (_terminal) return read_integer_from_terminal(value); if (!prepare_number()) return false; int c = static_cast(_buffer[_position++]); while (c <= ' ') c = static_cast(_buffer[_position++]); bool negative = false; if (c == '-') { negative = true; c = static_cast(_buffer[_position++]); } if constexpr (internal::is_signed_v) { T result = 0; while ('0' <= c && c <= '9') { const int first = c - '0'; const int second = static_cast(_buffer[_position]) - '0'; if (0 <= second && second <= 9) { result = negative ? result * 100 - (first * 10 + second) : result * 100 + (first * 10 + second); ++_position; } else { result = negative ? result * 10 - first : result * 10 + first; } c = static_cast(_buffer[_position++]); } value = result; } else { T result = 0; while ('0' <= c && c <= '9') { const unsigned first = unsigned(c - '0'); const int second = static_cast(_buffer[_position]) - '0'; if (0 <= second && second <= 9) { result = result * 100 + T(first * 10 + unsigned(second)); ++_position; } else { result = result * 10 + T(first); } c = static_cast(_buffer[_position++]); } value = negative ? T(0) - result : result; } if (_position > _length) _position = _length; return true; } template std::enable_if_t, bool> read(T& value) { if (!skip_spaces()) return false; int c = read_char_raw(); bool negative = false; if (c == '-' || c == '+') { negative = c == '-'; c = read_char_raw(); } long double result = 0; while ('0' <= c && c <= '9') { result = result * 10 + (c - '0'); c = read_char_raw(); } if (c == '.') { long double place = 0.1L; c = read_char_raw(); while ('0' <= c && c <= '9') { result += (c - '0') * place; place *= 0.1L; c = read_char_raw(); } } if (c == 'e' || c == 'E') { c = read_char_raw(); bool exponent_negative = false; if (c == '-' || c == '+') { exponent_negative = c == '-'; c = read_char_raw(); } int exponent = 0; while ('0' <= c && c <= '9') { exponent = exponent * 10 + (c - '0'); c = read_char_raw(); } long double scale = 1; long double power = 10; while (exponent > 0) { if (exponent & 1) scale *= power; power *= power; exponent >>= 1; } result = exponent_negative ? result / scale : result * scale; } value = static_cast(negative ? -result : result); return true; } template std::enable_if_t< internal::has_val_method_v && !internal::is_integral_v && !internal::is_range_v, bool > read(T& value) { long long x; if (!read(x)) return false; if constexpr (internal::has_static_mod_raw_v) { if (x >= 0 && uint64_t(x) < uint64_t(T::mod())) { value = T::raw(uint32_t(x)); } else { value = T(x); } } else { value = T(x); } return true; } template bool read(std::pair& value) { if (!read(value.first)) return false; return read(value.second); } template std::enable_if_t< internal::is_range_v && !internal::is_string_like_v, bool > read(Range& range) { using StoredValue = internal::range_stored_value_t; constexpr bool nested = internal::is_range_v && !internal::is_string_like_v; for (auto&& value : range) { if constexpr (std::is_same_v && !nested) { bool x; if (!read(x)) return false; value = x; } else { if (!read(value)) return false; } } return true; } template bool read(First& first, Second& second, Rest&... rest) { if (!read(first)) return false; return read(second, rest...); } template FastInput& operator>>(T& value) { if (!read(value)) std::abort(); return *this; } }; struct FastOutput { static constexpr int buffer_size = 1 << 20; private: inline static const auto digit_quads = [] { std::array result{}; for (int i = 0; i < 10000; i++) { int value = i; for (int j = 3; j >= 0; j--) { result[4 * i + j] = char('0' + value % 10); value /= 10; } } return result; }(); std::FILE* _stream; char _buffer[buffer_size]; int _position; int _precision; std::chars_format _float_format; char _range_separator; public: explicit FastOutput(std::FILE* stream = stdout) : _stream(stream), _position(0), _precision(6), _float_format(std::chars_format::general), _range_separator(' ') {} FastOutput(const FastOutput&) = delete; FastOutput& operator=(const FastOutput&) = delete; ~FastOutput() { flush(); } void flush() { if (_position == 0) return; std::fwrite(_buffer, 1, _position, _stream); _position = 0; } void write_char(char c) { if (_position == buffer_size) flush(); _buffer[_position++] = c; } void write(const char* s) { while (*s != '\0') write_char(*s++); } void write(const std::string& s) { for (char c : s) write_char(c); } void write(char c) { write_char(c); } void write(bool value) { write_char(value ? '1' : '0'); } template std::enable_if_t> write(T value) { char digits[128]; auto [end, error] = std::to_chars( digits, digits + sizeof(digits), value, _float_format, _precision ); if (error != std::errc()) std::abort(); for (const char* pointer = digits; pointer != end; pointer++) { write_char(*pointer); } } template std::enable_if_t< internal::is_integral_v && !std::is_same_v, bool> && !std::is_same_v, char> > write(T value) { using Raw = std::remove_cv_t; using Unsigned = internal::make_unsigned_t; Unsigned magnitude; if constexpr (internal::is_signed_v) { if (value < 0) { write_char('-'); magnitude = Unsigned(0) - Unsigned(value); } else { magnitude = Unsigned(value); } } else { magnitude = value; } if (magnitude == 0) { write_char('0'); return; } unsigned chunks[16]; int count = 0; while (magnitude >= 10000) { const Unsigned quotient = magnitude / 10000; chunks[count++] = unsigned(magnitude - quotient * 10000); magnitude = quotient; } if (_position > buffer_size - 64) flush(); const unsigned leading = unsigned(magnitude); const char* first = digit_quads.data() + 4 * leading; int skip = leading < 10 ? 3 : leading < 100 ? 2 : leading < 1000 ? 1 : 0; for (; skip < 4; skip++) _buffer[_position++] = first[skip]; while (count--) { const char* digits = digit_quads.data() + 4 * chunks[count]; std::memcpy(_buffer + _position, digits, 4); _position += 4; } } template std::enable_if_t< internal::has_val_method_v && !internal::is_integral_v && !internal::is_range_v > write(const T& value) { write(value.val()); } template void write(const std::pair& value) { write(value.first); write_char(' '); write(value.second); } template std::enable_if_t< internal::is_range_v && !internal::is_string_like_v > write(const Range& range) { using StoredValue = internal::range_stored_value_t; constexpr bool nested = internal::is_range_v && !internal::is_string_like_v; bool first = true; for (const auto& value : range) { if (!first) write_char(nested ? '\n' : _range_separator); first = false; if constexpr (std::is_same_v && !nested) { write(static_cast(value)); } else { write(value); } } } template void print(const First& first, const Rest&... rest) { write(first); ((write_char(' '), write(rest)), ...); } void println() { write_char('\n'); } void set_precision(int precision) { _precision = precision; } void set_fixed(int precision = 6) { _float_format = std::chars_format::fixed; _precision = precision; } void set_general(int precision = 6) { _float_format = std::chars_format::general; _precision = precision; } void set_range_separator(char separator) { _range_separator = separator; } template void println(const Args&... args) { print(args...); write_char('\n'); } template FastOutput& operator<<(const T& value) { write(value); return *this; } }; } // namespace utilities } // namespace m1une // END: utilities/fast_io.hpp #line 103 "template.hpp" using namespace std; namespace m1une { namespace template_io { inline utilities::FastInput& input() { static utilities::FastInput instance; return instance; } inline utilities::FastOutput& output() { static utilities::FastOutput instance; return instance; } } // namespace template_io } // namespace m1une using ll = long long; using u32 = unsigned int; using u64 = unsigned long long; using i128 = __int128; using u128 = unsigned __int128; #ifdef __SIZEOF_FLOAT128__ using f128 = __float128; #endif template constexpr T infty = 0; template <> constexpr int infty = 1'000'000'000; template <> constexpr ll infty = ll(infty) * infty * 2; template <> constexpr u32 infty = infty; template <> constexpr u64 infty = infty; template <> constexpr i128 infty = i128(infty) * infty; template <> constexpr double infty = infty; template <> constexpr long double infty = infty; using pi = pair; using pl = pair; using vi = vector; using vl = vector; template using vc = vector; template using vvc = vector>; using vvi = vvc; using vvl = vvc; template using vvvc = vector>; template using vvvvc = vector>; template using vvvvvc = vector>; template using pqg = std::priority_queue, greater>; template using umap = unordered_map; // template // using tree = __gnu_pbds::tree, // __gnu_pbds::rb_tree_tag, // __gnu_pbds::tree_order_statistics_node_update>; #define vv(type, name, h, ...) vector> name(h, vector(__VA_ARGS__)) #define vvv(type, name, h, w, ...) \ vector>> name(h, vector>(w, vector(__VA_ARGS__))) #define vvvv(type, name, a, b, c, ...) \ vector>>> name( \ a, vector>>(b, vector>(c, vector(__VA_ARGS__)))) #define overload4(a, b, c, d, e, ...) e #define overload3(a, b, c, d, ...) d // FOR(a) := for (ll _ = 0; _ < (ll)a; ++_) // FOR(i, a) := for (ll i = 0; i < (ll)a; ++i) // FOR(i, a, b) := for (ll i = a; i < (ll)b; ++i) // FOR(i, a, b, c) := for (ll i = a; i < (ll)b; i += (c)) // FOR_R(a) := for (ll i = (a) - 1; i >= 0; --i) // FOR_R(i, a) := for (ll i = (a) - 1; i >= 0; --i) // FOR_R(i, a, b) := for (ll i = (b) - 1; i >= (ll)a; --i) #define FOR1(a) for (ll _ = 0; _ < (ll)a; ++_) #define FOR2(i, a) for (ll i = 0; i < (ll)a; ++i) #define FOR3(i, a, b) for (ll i = a; i < (ll)b; ++i) #define FOR4(i, a, b, c) for (ll i = a; i < (ll)b; i += (c)) #define FOR1_R(a) for (ll i = (a) - 1; i >= 0; --i) #define FOR2_R(i, a) for (ll i = (a) - 1; i >= 0; --i) #define FOR3_R(i, a, b) for (ll i = (b) - 1; i >= (ll)a; --i) #define FOR(...) overload4(__VA_ARGS__, FOR4, FOR3, FOR2, FOR1)(__VA_ARGS__) #define FOR_R(...) overload3(__VA_ARGS__, FOR3_R, FOR2_R, FOR1_R)(__VA_ARGS__) #define FORI1(a) for (int _ = 0; _ < (int)a; ++_) #define FORI2(i, a) for (int i = 0; i < (int)a; ++i) #define FORI3(i, a, b) for (int i = a; i < (int)b; ++i) #define FORI4(i, a, b, c) for (int i = a; i < (int)b; i += (c)) #define FORI1_R(a) for (int i = (a) - 1; i >= 0; --i) #define FORI2_R(i, a) for (int i = (a) - 1; i >= 0; --i) #define FORI3_R(i, a, b) for (int i = (b) - 1; i >= (int)a; --i) #define FORI(...) overload4(__VA_ARGS__, FORI4, FORI3, FORI2, FORI1)(__VA_ARGS__) #define FORI_R(...) overload3(__VA_ARGS__, FORI3_R, FORI2_R, FORI1_R)(__VA_ARGS__) #define FOR_subset(t, s) for (int t = (s); t >= 0; t = (t == 0 ? -1 : (t - 1) & (s))) #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() int popcnt(int x) { return __builtin_popcount(x); } int popcnt(u32 x) { return __builtin_popcount(x); } int popcnt(ll x) { return __builtin_popcountll(x); } int popcnt(u64 x) { return __builtin_popcountll(x); } int popcnt_mod_2(int x) { return __builtin_parity(x); } int popcnt_mod_2(u32 x) { return __builtin_parity(x); } int popcnt_mod_2(ll x) { return __builtin_parityll(x); } int popcnt_mod_2(u64 x) { return __builtin_parityll(x); } // (0, 1, 2, 3, 4) -> (-1, 0, 1, 1, 2) int topbit(int x) { return (x == 0 ? -1 : 31 - __builtin_clz(x)); } int topbit(u32 x) { return (x == 0 ? -1 : 31 - __builtin_clz(x)); } int topbit(ll x) { return (x == 0 ? -1 : 63 - __builtin_clzll(x)); } int topbit(u64 x) { return (x == 0 ? -1 : 63 - __builtin_clzll(x)); } // (0, 1, 2, 3, 4) -> (-1, 0, 1, 0, 2) int lowbit(int x) { return (x == 0 ? -1 : __builtin_ctz(x)); } int lowbit(u32 x) { return (x == 0 ? -1 : __builtin_ctz(x)); } int lowbit(ll x) { return (x == 0 ? -1 : __builtin_ctzll(x)); } int lowbit(u64 x) { return (x == 0 ? -1 : __builtin_ctzll(x)); } template T floor(T a, T b) { return a / b - (a % b && (a ^ b) < 0); } template T ceil(T x, T y) { return floor(x + y - 1, y); } template T bmod(T x, T y) { return x - y * floor(x, y); } template pair divmod(T x, T y) { T q = floor(x, y); return {q, x - q * y}; } template T POW(U x_, int n) { T x = x_; T ret = 1; while (n > 0) { if (n & 1) ret *= x; x *= x; n >>= 1; } return ret; } template T SUM(const vector& A) { T sm = 0; for (auto&& a : A) sm += a; return sm; } #define LB(c, x) distance((c).begin(), lower_bound(all(c), (x))) #define UB(c, x) distance((c).begin(), upper_bound(all(c), (x))) #define UNIQUE(x) sort(all(x)), x.erase(unique(all(x)), x.end()), x.shrink_to_fit() template inline bool chmax(T& a, const S& b) { return (a < b ? a = b, 1 : 0); } template inline bool chmin(T& a, const S& b) { return (a > b ? a = b, 1 : 0); } // ? は -1 vc s_to_vi(const string& S, char first_char) { vc A(S.size()); FOR(i, S.size()) { A[i] = (S[i] != '?' ? S[i] - first_char : -1); } return A; } template vector cumsum(vector& A, int off = 1) { int N = A.size(); vector B(N + 1); FOR(i, N) { B[i + 1] = B[i] + A[i]; } if (off == 0) B.erase(B.begin()); return B; } template vector argsort(const vector& A) { vector ids(A.size()); iota(all(ids), 0); sort(all(ids), [&](int i, int j) { return (A[i] == A[j] ? i < j : A[i] < A[j]); }); return ids; } // A[I[0]], A[I[1]], ... template vc rearrange(const vc& A, const vc& I) { vc B(I.size()); FOR(i, I.size()) B[i] = A[I[i]]; return B; } template constexpr auto min(T... a) { return min(initializer_list>{a...}); } template constexpr auto max(T... a) { return max(initializer_list>{a...}); } template void print(const Ts&... values) { m1une::template_io::output().println(values...); } void YESNO(bool b) { m1une::template_io::output().println(b ? "YES" : "NO"); } void YesNo(bool b) { m1une::template_io::output().println(b ? "Yes" : "No"); } void YES() { m1une::template_io::output().println("YES"); } void NO() { m1une::template_io::output().println("NO"); } void Yes() { m1une::template_io::output().println("Yes"); } void No() { m1une::template_io::output().println("No"); } // END: template.hpp #line 29 "pch.hpp" // END: pch.hpp #line 2 "..::contest594::H::main.cpp" auto& fastin = m1une::template_io::input(); auto& fastout = m1une::template_io::output(); // BEGIN: math/modint.hpp #line 3 "math::modint.hpp" #include #include #include #include #include namespace m1une { namespace math { template struct ModInt { static_assert(0 < Modulus, "Modulus must be positive"); private: uint32_t _v; public: static constexpr uint32_t mod() { return Modulus; } static constexpr ModInt raw(uint32_t v) noexcept { ModInt x; x._v = v; return x; } constexpr ModInt() noexcept : _v(0) {} template , int> = 0> constexpr ModInt(Integer v) noexcept { if constexpr (std::is_signed_v) { int64_t x = static_cast(v) % static_cast(Modulus); if (x < 0) x += Modulus; _v = static_cast(x); } else { _v = static_cast(static_cast(v) % Modulus); } } constexpr uint32_t val() const noexcept { return _v; } constexpr ModInt& operator++() noexcept { _v++; if (_v == Modulus) _v = 0; return *this; } constexpr ModInt& operator--() noexcept { if (_v == 0) _v = Modulus; _v--; return *this; } constexpr ModInt operator++(int) noexcept { ModInt res = *this; ++*this; return res; } constexpr ModInt operator--(int) noexcept { ModInt res = *this; --*this; return res; } constexpr ModInt& operator+=(const ModInt& rhs) noexcept { _v += rhs._v; if (_v >= Modulus) _v -= Modulus; return *this; } constexpr ModInt& operator-=(const ModInt& rhs) noexcept { _v -= rhs._v; if (_v >= Modulus) _v += Modulus; return *this; } constexpr ModInt& operator*=(const ModInt& rhs) noexcept { uint64_t z = _v; z *= rhs._v; _v = static_cast(z % Modulus); return *this; } constexpr ModInt& operator/=(const ModInt& rhs) noexcept { return *this *= rhs.inv(); } 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 bool operator==(const ModInt& rhs) const noexcept { return _v == rhs._v; } constexpr bool operator!=(const ModInt& rhs) const noexcept { return _v != rhs._v; } constexpr ModInt pow(long long n) const noexcept { ModInt res = raw(1), x = *this; while (n > 0) { if (n & 1) res *= x; x *= x; n >>= 1; } return res; } constexpr ModInt inv() const noexcept { int64_t a = _v, b = Modulus, u = 1, v = 0; while (b) { int64_t t = a / b; a -= t * b; std::swap(a, b); u -= t * v; std::swap(u, v); } if (u < 0) u += Modulus; return raw(static_cast(u)); } friend std::ostream& operator<<(std::ostream& os, const ModInt& rhs) { return os << rhs._v; } friend std::istream& operator>>(std::istream& is, ModInt& rhs) { long long v; is >> v; rhs = ModInt(v); return is; } }; using modint998244353 = ModInt<998244353>; using modint1000000007 = ModInt<1000000007>; template struct DynamicModInt { private: uint32_t _v; inline static uint32_t _mod = 1; public: static uint32_t mod() noexcept { return _mod; } static void set_mod(uint32_t modulus) noexcept { assert(modulus > 0); assert(modulus <= uint32_t(1) << 31); _mod = modulus; } static DynamicModInt raw(uint32_t v) noexcept { assert(v < _mod); DynamicModInt x; x._v = v; return x; } DynamicModInt() noexcept : _v(0) {} template , int> = 0> DynamicModInt(Integer v) noexcept { if constexpr (std::is_signed_v) { int64_t x = static_cast(v) % static_cast(_mod); if (x < 0) x += _mod; _v = static_cast(x); } else { _v = static_cast(static_cast(v) % _mod); } } uint32_t val() const noexcept { return _v; } DynamicModInt& operator++() noexcept { _v++; if (_v == _mod) _v = 0; return *this; } DynamicModInt& operator--() noexcept { if (_v == 0) _v = _mod; _v--; return *this; } DynamicModInt operator++(int) noexcept { DynamicModInt result = *this; ++*this; return result; } DynamicModInt operator--(int) noexcept { DynamicModInt result = *this; --*this; return result; } DynamicModInt& operator+=(const DynamicModInt& rhs) noexcept { _v += rhs._v; if (_v >= _mod) _v -= _mod; return *this; } DynamicModInt& operator-=(const DynamicModInt& rhs) noexcept { _v -= rhs._v; if (_v >= _mod) _v += _mod; return *this; } DynamicModInt& operator*=(const DynamicModInt& rhs) noexcept { _v = static_cast(uint64_t(_v) * rhs._v % _mod); return *this; } DynamicModInt& operator/=(const DynamicModInt& rhs) noexcept { return *this *= rhs.inv(); } DynamicModInt operator+(const DynamicModInt& rhs) const noexcept { return DynamicModInt(*this) += rhs; } DynamicModInt operator-(const DynamicModInt& rhs) const noexcept { return DynamicModInt(*this) -= rhs; } DynamicModInt operator*(const DynamicModInt& rhs) const noexcept { return DynamicModInt(*this) *= rhs; } DynamicModInt operator/(const DynamicModInt& rhs) const noexcept { return DynamicModInt(*this) /= rhs; } bool operator==(const DynamicModInt& rhs) const noexcept { return _v == rhs._v; } bool operator!=(const DynamicModInt& rhs) const noexcept { return _v != rhs._v; } DynamicModInt pow(long long exponent) const noexcept { assert(exponent >= 0); DynamicModInt result = raw(1 % _mod); DynamicModInt base = *this; while (exponent > 0) { if (exponent & 1) result *= base; base *= base; exponent >>= 1; } return result; } DynamicModInt inv() const noexcept { int64_t a = _v, b = _mod, u = 1, v = 0; while (b) { int64_t quotient = a / b; a -= quotient * b; std::swap(a, b); u -= quotient * v; std::swap(u, v); } assert(a == 1); u %= _mod; if (u < 0) u += _mod; return raw(static_cast(u)); } friend std::ostream& operator<<(std::ostream& os, const DynamicModInt& rhs) { return os << rhs._v; } friend std::istream& operator>>(std::istream& is, DynamicModInt& rhs) { long long value; is >> value; rhs = DynamicModInt(value); return is; } }; } // namespace math } // namespace m1une // END: math/modint.hpp #line 7 "..::contest594::H::main.cpp" using mint = m1une::math::modint998244353; void solve() { ll hc, ac, sc, hg, ag, sg; fastin.read(hc, ac, sc, hg, ag, sg); { ll t = min((hc - 1) / ag, (hg - 1) / ac); hc -= t * ag; hg -= t * ac; } mint p = 1, ans = 0; // dump(p, ans, hc, ac, sc, hg, ag, sg); if (hc > ag) { assert(hg <= ac); ll t = (hc - 1) / ag; p = mint(10).pow(t).inv(); ans = mint(1) - p; } else if (hg > ac) { assert(hc <= ag); ll t = (hg - 1) / ac; p = mint(10).pow(t).inv(); } // dump(p, ans, hc, ac, sc, hg, ag, sg); if (sc == sg) { ans += p / mint(2); } else if (sc > sg) { ans += p * mint(10) / mint(11); } else { ans += p / mint(11); } print(ans); } int main() { CPP_DUMP_SET_OPTION(max_line_width, 80); CPP_DUMP_SET_OPTION(log_label_func, cpp_dump::log_label::filename()); CPP_DUMP_SET_OPTION(enable_asterisk, true); int T = 1; cin >> T; while (T--) solve(); return 0; } // END: ../contest594/H/main.cpp