#include using namespace std; const int MOD = 1e9 + 7; long long mod_pow(long long a,long long b,long long mod) { if( b == 0 ) return 1; long long res = mod_pow(a*a%mod,b/2,mod); if( b&1 ) res = (res*a)%mod; return res; } long long inv(long long a,long long mod) { return mod_pow(a, mod - 2, mod); } long long nCr(int n, int r, int p) { long long ret = 1; while(true){ if(r == 0)break; int N = n % p; int R = r % p; if(N < R)return 0; for(int i = 0;i < R;i++){ ret = ret * (N-i) % p; } long imul = 1; for(int i = 0;i < R;i++){ imul = imul * (i+1) % p; } ret = ret * inv(imul, p) % p; n /= p; r /= p; } return ret; } int main() { long long N, K; cin >> N >> K; // 包除原理ぽよねー long long ret = 0; for(int i = 1; i <= K; i++) { long long sign = pow(-1, K - i); (ret += sign * nCr(K, i, MOD) * mod_pow(i, N, MOD) + MOD) %= MOD; } cout << ret << endl; }