// 2022.10.21 // #define _GLIBCXX_DEBUG // // #define __USE_MATH_DEFINES // const double PI = acos(-1.0) #include using namespace std; // (*) snipet // type name using ll = long long; using pii = pair; using pil = pair; using pli = pair; using pll = pair; using vi = vector; using vvi = vector; using vvvi = vector; using vl = vector; using vvl = vector; using vvvl = vector; using vpii = vector; using vpil = vector; using vpli = vector; using vpll = vector; using vb = vector; using vvb = vector; using vvvb = vector; using vd = vector; using vvd = vector; using vvvd = vector; template using pq = priority_queue>; template using pql = priority_queue, greater>; // rep #define rep0(goal) for (int cnt = 0; cnt < int(goal); ++cnt) #define rep(cnt, goal) for (int cnt = 0; cnt < int(goal); ++cnt) #define rep2(cnt, start, goal) for (int cnt = int(start); cnt < int(goal); ++cnt) #define rep3(cnt, start, goal) for (int cnt = int(start); cnt > int(goal); --cnt) // all #define all(ctn) begin(ctn), end(ctn) // chmax, chmin template bool chmax(T &x, const T y) { if (x < y) { x = y; return true; } else { return false; }} template bool chmin(T &x, const T y) { if (x > y) { x = y; return true; } else { return false; }} // etc. void Yn(const bool &exp) { if (exp) cout << "Yes" << endl; else cout << "No" << endl;} void set_prec(const int &dig) { cout << fixed << setprecision(dig); cerr << fixed << setprecision(dig);} // template T INF() { return numeric_limits::max();} // template T MIN_INF() { return numeric_limits::min();} #include using namespace atcoder; // using mint = atcoder::modint1000000007; using mint = atcoder::modint998244353; // using mint = atcoder::modint; void pr(const mint &MINT) { cerr << MINT.val(); } // mint using pmm = pair; using pim = pair; using pmi = pair; using plm = pair; using pml = pair; using vm = vector; using vvm = vector; using vvvm = vector; // debug template void pr(const T &obj) { cerr << obj; } // single template void pr(const T &first, const Ts &...rest) { pr(first); cerr << ", "; pr(rest...); } // multi(plural) template void pr(const pair &PAIR) { cerr << "("; pr(PAIR.first); cerr << ", "; pr(PAIR.second); cerr << ")"; } // pair template void pr(const tuple &trp) { cerr << "("; pr(get<0>(trp)); cerr << ", "; pr(get<1>(trp)); cerr << ", "; pr(get<2>(trp)); cerr << ")"; } // tuple(3) template void pr(const vector &vec) { cerr << "{"; for (T obj : vec) { pr(obj); cerr << ", "; } cerr << "}"; } // vector template void pr(const vector> &vv) { rep(row, vv.size()) { cerr << endl; cerr << "[" << row << "]: "; pr(vv[row]); }} // vector(multi-D) template void pr(const set &SET) { cerr << "{"; for (T obj : SET) { pr(obj); cerr << ", "; } cerr << "}"; } // set template void pr(const multiset &MS) { cerr << "{"; for (T obj : MS) { pr(obj); cerr << ", "; } cerr << "}"; } // multiset template void pr(const map &MAP) { cerr << "{"; for (pair p : MAP) { pr(p.first); cerr << ": "; pr(p.second); cerr << ", "; } cerr << "}";} // map template void pr(const queue &que) { queue q = que; cerr << "{"; while (!q.empty()) { pr(q.front()); q.pop(); cerr << ", ";} cerr << "}";} // queue template void pr(const deque &deq) { deque d = deq; cerr << "{"; while (!d.empty()) { pr(d.front()); d.pop_front(); cerr << ", ";} cerr << "}";} // deque template void pr(const stack &stc) { stack s = stc; vector v; while (!s.empty()) { v.push_back(s.top()); s.pop();} reverse(all(v)); cerr << "{"; for (T obj : v) cerr << obj << ", "; cerr << "}";} // stack template void pr(const priority_queue &pq) { priority_queue p = pq; cerr << "{"; while (!p.empty()) { pr(p.top()); p.pop(); cerr << ", "; } cerr << "}";} // priority_queue template void pr(const priority_queue, greater> &pq) { priority_queue, greater> p = pq; cerr << "{"; while (!p.empty()) { pr(p.top()); p.pop(); cerr << ", ";} cerr << "}";} // priority_queue(from less) #define db(obj) cerr << #obj << ": "; pr(obj); cerr << " "; #define dl(obj) db(obj); cerr << endl; #define dm(...) cerr << "(" << #__VA_ARGS__ << "): ("; pr(__VA_ARGS__); cerr << ") "; #define dml(...) dm(__VA_ARGS__); cerr << endl; // N元集合からM元集合への全射の数 // N -> n, M -> k // global ll N; int M; vl fac; // fac[i] = i! vl inv; // inv[i] = 1/i! mod m ll m = 1e9 + 7; // 998244353; // MOD void input() { cin >> N >> M; } ll my_pow_mod(ll x, ll n, ll m) { assert(0 <= n); assert(1 <= m); ll res = 1; ll y = x % m; if (x < 0) x += m; while (n) { if (n % 2) { res = (res * y) % m; } n /= 2; y = (y * y) % m; } return res; } ll my_inv_mod(ll x, ll m) { assert(1 <= m); assert(x % m != 0); // m: 素数 return my_pow_mod(x, m - 2, m); // assert(my_gcd(x, m) == 1); // return my_pow_mod(x, phi(m) - 1, m); } // factorial, inv_factorial vl my_facts(int N) { assert(0 <= N); vl res(N + 1); res[0] = 1; rep(i, N) res[i + 1] = (res[i] * (i + 1)) % m; return res; } vl my_inv_facts(int N) { assert(0 <= N); vl res(N + 1); res[0] = 1; rep(i, N) res[i + 1] = (res[i] * my_inv_mod(i + 1, m)) % m; return res; } void init() { fac = my_facts(M + 1); inv = my_inv_facts(M + 1); } ll surjection(ll n, int k) { // O(k*logn) assert(0 <= n); assert(0 <= k); // if (n < k) return 0; ll res = 0; rep(i, k + 1) { ll temp = ((fac[k] * inv[k - i] % m) * inv[i] % m); temp = (temp * my_pow_mod(i, n, m)) % m; if ((k - i) % 2 == 0) res = (res + temp) % m; else res = (res - temp + m) % m; } return res; } void solve() { cout << surjection(N, M) << endl; }; int main() { input(); init(); solve(); }