#include #include const long long MOD = 1000000007; int main() { int N, M; std::cin >> N >> M; auto pow = [](long long a, int n) { long long res = 1; for (; n; n >>= 1) { if (n & 1) res = res * a % MOD; a = a * a % MOD; } return res; }; std::vector table(N + 1), inv_table(N + 1); table[0] = 1; for (int i = 1; i <= N; i++) table[i] = table[i - 1] * i % MOD; inv_table[N] = pow(table[N], MOD - 2); for (int i = N; i > 0; i--) inv_table[i - 1] = inv_table[i] * i % MOD; auto comb = [&table, &inv_table](int n, int r) { if (n < r) return 0LL; long long tmp = inv_table[n - r] * inv_table[r] % MOD; return table[n] * tmp % MOD; }; long long ans = pow(M, N); for (int i = M - 1; i > 0; i--) { if ((M + i) % 2) { ans -= comb(M, i) * pow(i, N) % MOD; ans += ans < 0 ? MOD : 0; } else { ans += comb(M, i) * pow(i, N) % MOD; ans %= MOD; } } std::cout << ans << "\n"; }