const MOD: usize = 998244353; fn power(base: usize, times: usize) -> usize { if times == 0 { return 1usize; } if times == 1 { return base; } let temp = power(base, times/2); temp * temp % MOD * power(base, times%2) % MOD } fn main() { let mut nm = String::new(); std::io::stdin().read_line(&mut nm).ok(); let nm: Vec = nm.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let n = nm[0]; let m = nm[1]; let mut powers = vec![0usize; m+1]; let p = n % (MOD-1); for i in 1..=m { powers[i] = power(i, p); } let nmod = n % MOD; let mut result = 0usize; let summary = (1+m) * nmod % MOD * power(2, MOD-2) % MOD; for i in 1..=m { result += i * (MOD + powers[i] - powers[i-1]) % MOD; result += MOD - i * (MOD + powers[m-i+1] - powers[m-i]) % MOD; result %= MOD; } println!("{}", result * summary % MOD); }