#include #define REP(i, n) for(int i = 0;i < n;i++) #define VSORT(v) sort(v.begin(), v.end()) #define VRSORT(v) sort(v.rbegin(), v.rend()) #define ll long long using namespace std; typedef pair P; typedef pair LP; typedef pair PP; typedef pair LPP; //typedef vectorvec; //typedef vectorvec; //typedef vector mat; typedef vector> Graph; const int dx[8] = {1, 0, -1, 0, 1, -1, -1, 1}; const int dy[8] = {0, 1, 0, -1, 1, 1, -1, -1}; const int INF = 1000000000; const ll LINF = 1000000000000000000;//1e18 const ll MOD = 1000000007; const double PI = acos(-1.0); const double EPS = 1e-10; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template inline void add(T &a, T b){a = ((a+b) % MOD + MOD) % MOD;}; vector > prime_factorize(long long n) { vector > res; for (long long p = 2; p * p <= n; ++p) { if (n % p != 0) continue; int num = 0; while (n % p == 0) { ++num; n /= p; } res.push_back(make_pair(p, num));//p^num } if (n != 1) res.push_back(make_pair(n, 1)); return res; } const int MAX = 1110000; long long fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void comb_init() { fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (int i = 2; i < MAX; i++){ fac[i] = fac[i - 1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD; finv[i] = finv[i - 1] * inv[i] % MOD; } } // 二項係数計算 long long comb(int n, int k){ if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD; } ll mod_pow(ll x, ll n, ll mod){ ll res = 1; while(n > 0){ if(n & 1) res = res * x % mod;//最下位ビットが立っているときにx^(2^i)をかける x = x * x % mod; //xを順次二乗していく n >>= 1; } return res; } int main(){ cin.tie(0); ios::sync_with_stdio(false); ll n, m; cin >> n >> m; comb_init(); ll ans = 0; for(int k=0;k<=m;k++){ ll tmp = comb(m, k) * mod_pow(m - k, n, MOD); tmp %= MOD; if(k % 2 == 0) ans += tmp; else{ ans -= tmp; ans += MOD; } ans %= MOD; } cout << ans << endl; }