from functools import reduce
n, m = map(int, input().split())
mo = 998244353
f = lambda m, n: reduce(lambda x, y: x * y % mo, ((pow(2, m, mo) - pow(2, i, mo)) for i in range(n)))
fact = lambda n: reduce(lambda x, y: x * y % mo, range(1, n + 1))
rev = lambda x: pow(x, mo - 2, mo)
print(pow(2, n, mo), f(m, n) * rev(fact(n)) % mo, f(m, n) * rev(f(n, n)) % mo)