use std::io::*; const MOD: u64 = 1_000_000_007; mod enumeration { pub struct Enumeration { fact: Vec, finv: Vec, modulo: u64, } impl Enumeration { pub fn new(n: usize, modulo: u64) -> Enumeration { let mut fact: Vec = vec![1; n]; let mut inv: Vec = vec![1; n]; let mut finv: Vec = vec![1; n]; for i in 2..n { fact[i] = fact[i - 1] * (i as u64) % modulo; inv[i] = modulo - (modulo / (i as u64)) * inv[(modulo as usize % i)] % modulo; finv[i] = inv[i] * finv[i - 1] % modulo; } Enumeration { fact: fact, finv: finv, modulo: modulo, } } pub fn c(&self, n: usize, k: usize) -> u64 { if n < k { return 0; } self.fact[n] * self.finv[n - k] % self.modulo * self.finv[k] % self.modulo } } } fn mod_pow(mut x: u64, mut e: u64, m: u64) -> u64 { let mut res = 1; while e > 0 { if e & 1 == 1 { res = res * x % m; } x = x * x % m; e >>= 1; } res } fn main() { let mut s: String = String::new(); std::io::stdin().read_to_string(&mut s).ok(); let mut itr = s.trim().split_whitespace(); let n: u64 = itr.next().unwrap().parse().unwrap(); let m: u64 = itr.next().unwrap().parse().unwrap(); let en = enumeration::Enumeration::new(100010, MOD); let mut ans = mod_pow(m, n, MOD); for i in (1..m).rev() { let c = en.c(m as usize, i as usize); if (m - i) % 2 == 1 { ans = (ans + MOD - c * mod_pow(i, n, MOD) % MOD) % MOD; } else { ans = (ans + c * mod_pow(i, n, MOD) % MOD) % MOD; } } println!("{}", ans); }