#include #include #include using namespace std; typedef long long ll; const int MOD = 998244353; struct mint { int val; mint() : val(0) {} mint(long long v) { if (abs(v) >= MOD) { v %= MOD; } if (v < 0) { v += MOD; } val = v; } mint &operator++() { val++; if (val == MOD) { val = 0; } return *this; } mint &operator--() { if (val == 0) { val = MOD; } val--; return *this; } mint &operator+=(const mint &x) { val += x.val; if (val >= MOD) { val -= MOD; } return *this; } mint &operator-=(const mint &x) { val -= x.val; if (val < 0) { val += MOD; } return *this; } mint &operator*=(const mint &x) { val = (int)((long long)val * x.val % MOD); return *this; } mint &operator/=(const mint &x) { *this *= x.inv(); return *this; } mint operator-() { return mint() - *this; } mint pow(long long n) const { mint x = 1, r = *this; while (n) { if (n & 1) { x *= r; } r *= r; n >>= 1; } return x; } mint inv() const { return pow(MOD - 2); } friend mint operator+(const mint &x, const mint &y) { return mint(x) += y; } friend mint operator-(const mint &x, const mint &y) { return mint(x) -= y; } friend mint operator*(const mint &x, const mint &y) { return mint(x) *= y; } friend mint operator/(const mint &x, const mint &y) { return mint(x) /= y; } friend bool operator==(const mint &x, const mint &y) { return x.val == y.val; } friend bool operator!=(const mint &x, const mint &y) { return x.val != y.val; } friend std::ostream &operator<<(std::ostream &os, const mint &x) { return os << x.val; } friend std::istream &operator>>(std::istream &is, mint &x) { int v; is >> v; x = mint(v); return is; } }; mint modpow(mint x, long long n) { mint res = 1, r = x; while (n) { if (n & 1) { res *= r; } r *= r; n >>= 1; } return res; } struct combination_big { long long maxn; int maxk; std::vector nP; std::vector inP; std::vector fact; std::vector ifact; combination_big(long long maxn, int maxk) : maxn(maxn), maxk(maxk), nP(maxk + 1), inP(maxk + 1), fact(maxk + 1), ifact(maxk + 1) { nP[0] = 1; inP[0] = 1; fact[0] = 1; ifact[0] = 1; for (int i = 1; i <= maxk; i++) { nP[i] = nP[i - 1] * (maxn + 1 - i); inP[i] = inP[i - 1] * modpow(maxn + 1 - i, MOD - 2); fact[i] = fact[i - 1] * i; ifact[i] = ifact[i - 1] * modpow(i, MOD - 2); } } mint operator()(long long n, int k) { assert(n <= maxn); assert(maxn - n + k <= maxk); if (k > n || k < 0) { return 0; } return nP[maxn - n + k] * inP[maxn - n] * ifact[k]; } }; struct combination_fix_n { long long n; int maxk; std::vector nP; std::vector fact; std::vector ifact; combination_fix_n(long long n, int maxk) : n(n), maxk(maxk) { nP.resize(maxk + 1); fact.resize(maxk + 1); ifact.resize(maxk + 1); nP[0] = 1; fact[0] = 1; for (int i = 1; i <= maxk; i++) { nP[i] = nP[i - 1] * (n + 1 - i); fact[i] = fact[i - 1] * i; } ifact[maxk] = modpow(fact[maxk], MOD - 2); for (int i = maxk - 1; i >= 0; i--) { ifact[i] = ifact[i + 1] * (i + 1); } } mint operator()(int k) { assert(k <= maxk); if (k > n || k < 0) { return 0; } return nP[k] * ifact[k]; } }; int main() { ll n, m; cin >> n >> m; combination_big com(n + 120, m + 300); mint ans = modpow(2, n); for (int i = 0; i < m; i++) { ans -= com(n, i); } cout << ans << endl; }