#include using namespace std; using ll = long long; using P = pair; #define rep(i, a, b) for(long long i = (a); i < (b); ++i) #define rrep(i, a, b) for(long long i = (a); i >= (b); --i) constexpr long long inf = 4e18; struct SetupIO { SetupIO() { ios::sync_with_stdio(0); cin.tie(0); cout << fixed << setprecision(30); } } setup_io; template struct StaticModint { using mint = StaticModint; static constexpr uint32_t mod() { return m; } static constexpr mint raw(const uint32_t v) { mint a; a._v = v; return a; } constexpr StaticModint() : _v(0) {} template constexpr StaticModint(const T& v) { static_assert(is_integral_v); if constexpr(is_signed_v) { int64_t x = int64_t(v % int64_t(m)); if(x < 0) x += m; _v = uint32_t(x); } else _v = uint32_t(v % m); } constexpr uint32_t val() const { return _v; } constexpr mint& operator++() { return *this += 1; } constexpr mint& operator--() { return *this -= 1; } constexpr mint operator++(int) { mint res = *this; ++*this; return res; } constexpr mint operator--(int) { mint res = *this; --*this; return res; } constexpr mint& operator+=(mint rhs) { if(_v >= m - rhs._v) _v -= m; _v += rhs._v; return *this; } constexpr mint& operator-=(mint rhs) { if(_v < rhs._v) _v += m; _v -= rhs._v; return *this; } constexpr mint& operator*=(mint rhs) { return *this = *this * rhs; } constexpr mint& operator/=(mint rhs) { return *this *= rhs.inv(); } constexpr mint operator+() const { return *this; } constexpr mint operator-() const { return mint{} - *this; } constexpr mint pow(long long n) const { assert(0 <= n); if(n == 0) return 1; mint x = *this, r = 1; while(1) { if(n & 1) r *= x; n >>= 1; if(n == 0) return r; x *= x; } } constexpr mint inv() const { if constexpr(prime) { assert(_v); return pow(m - 2); } else { const auto eg = inv_gcd(_v, m); assert(eg.first == 1); return eg.second; } } friend constexpr mint operator+(mint lhs, mint rhs) { return lhs += rhs; } friend constexpr mint operator-(mint lhs, mint rhs) { return lhs -= rhs; } friend constexpr mint operator*(mint lhs, mint rhs) { return uint64_t(lhs._v) * rhs._v; } friend constexpr mint operator/(mint lhs, mint rhs) { return lhs /= rhs; } friend constexpr bool operator==(mint lhs, mint rhs) { return lhs._v == rhs._v; } friend constexpr bool operator!=(mint lhs, mint rhs) { return lhs._v != rhs._v; } friend istream& operator>>(istream& in, mint& x) { long long a; in >> a; x = a; return in; } friend ostream& operator<<(ostream& out, const mint& x) { return out << x.val(); } private: uint32_t _v = 0; static constexpr bool prime = []() -> bool { if(m == 1) return 0; if(m == 2 or m == 7 or m == 61) return 1; if(m % 2 == 0) return 0; uint32_t d = m - 1; while(d % 2 == 0) d /= 2; for(uint32_t a : {2, 7, 61}) { uint32_t t = d; mint y = mint(a).pow(t); while(t != m - 1 && y != 1 && y != m - 1) { y *= y; t <<= 1; } if(y != m - 1 && t % 2 == 0) return 0; } return 1; }(); static constexpr pair inv_gcd(const int32_t a, const int32_t b) { if(a == 0) return {b, 0}; int32_t s = b, t = a, m0 = 0, m1 = 1; while(t) { const int32_t u = s / t; s -= t * u; m0 -= m1 * u; swap(s, t); swap(m0, m1); } if(m0 < 0) m0 += b / s; return {s, m0}; } }; using modint998244353 = StaticModint<998244353>; using modint1000000007 = StaticModint<1000000007>; vector divisor(const long long n) { assert(n >= 1); vector res; for(long long i = 1; i * i <= n; ++i) { if(n % i == 0) { res.push_back(i); if(i * i != n) res.emplace_back(n / i); } } sort(res.begin(), res.end()); return res; } using mint = modint998244353; template struct Binomial { Binomial(int n) : fac(n + 1), ifac(n + 1) { fac[0] = 1; for(int i = 1; i <= n; ++i) fac[i] = fac[i - 1] * i; ifac[n] = fac[n].inv(); for(int i = n; i >= 1; --i) ifac[i - 1] = ifac[i] * i; } mint fact(int n) const { if(n < 0) return 0; return fac[n]; } mint perm(int n, int r) const { if(n < 0 or n < r or r < 0) return 0; return fac[n] * ifac[n - r]; } mint comb(int n, int r) const { if(n < 0 or n < r or r < 0) return 0; return fac[n] * ifac[n - r] * ifac[r]; } mint homo(int n, int r) const { if(n < 0 or r < 0) return 0; if(r == 0) return 1; return comb(n + r - 1, r); } mint operator()(int n, int r) const { return comb(n, r); } private: vector fac, ifac; }; Binomial binom(1000005); int main(void) { ll n, m; cin >> n >> m; auto func = [&](const ll x) -> mint { if(n % 2 != abs(x) % 2) return 0; return binom(n, (n + x) / 2) / mint(2).pow(n); }; if(m == 0) { cout << func(0) * 2 - func(0) * func(0) << '\n'; return 0; } mint ans = 0; for(const ll d : divisor(abs(m))) { ans += func(d) * func(m / d); } cout << ans * 2 << '\n'; }