#pragma region Macros #include using namespace std; using namespace __gnu_pbds; // using namespace __gnu_cxx; #define TO_STRING(var) # var #define pb emplace_back #define int ll #define endl '\n' using ll = long long; using ld = long double; const ld PI = acos(-1); const ld EPS = 1e-10; const int INF = 1 << 30; const ll INFL = 1LL << 62; const int MOD = 998244353; // const int MOD = 1000000007; int ceil(int x, int y) {return (x > 0 ? (x + y - 1) / y : x / y);} int POW(int x, int y) { if (y != (int)(y) or y < 0 or x != 0 && x != 1 && y > 64) {cout << "Error" << endl;return 0;} if (y == 0) return 1; if (y % 2 == 0) return POW(x * x, y / 2); return x * POW(x, y - 1); } __attribute__((constructor)) void constructor() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); } #pragma endregion template class modint{ public: int val = 0; modint(int x = 0) { while (x < 0) x += mod; val = x % mod; } modint(const modint &r) { val = r.val; } // コピーコンストラクタ modint operator -(){ return modint(-val); } // 単項 modint operator +(const modint &r) { return modint(*this) += r; } modint operator -(const modint &r) { return modint(*this) -= r; } modint operator *(const modint &r) { return modint(*this) *= r; } modint operator /(const modint &r) { return modint(*this) /= r; } modint &operator +=(const modint &r) { val += r.val; if (val >= mod) val -= mod; return *this; } modint &operator -=(const modint &r) { if (val < r.val) val += mod; val -= r.val; return *this; } modint &operator *=(const modint &r) { val = val * r.val % mod; return *this; } modint &operator /=(const modint &r) { int a = r.val, b = mod, u = 1, v = 0; while (b) { int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v); } val = val * u % mod; if (val < 0) val += mod; return *this; } bool operator ==(const modint& r) { return this -> val == r.val; } bool operator <(const modint& r) { return this -> val < r.val; } bool operator !=(const modint& r) { return this -> val != r.val; } }; using mint = modint; istream &operator >>(istream &is, mint& x) { int t; is >> t; x = t; return (is); } ostream &operator <<(ostream &os, const mint& x) { return os << x.val; } mint modpow(const mint &a, int n) { if (n == 0) return 1; mint t = modpow(a, n / 2); t = t * t; if (n & 1) t = t * a; return t; } signed main() { int N, M; cin >> N >> M; if (M > N) {cout << 0 << endl; return 0;} if (M == N) {cout << 1 << endl; return 0;} mint zen = 1; zen = modpow(2, N); vector R(1e6 + 1); R[0] = 1; for (int i = 1; i <= 1e6; i++) { R[i] = R[i - 1] * i; } mint yo = 1; for (int i = N; i >= N - M + 2; i--) { yo *= i; } yo /= R[M - 1]; cout << zen - yo << endl; }