MOD = 998244353
N,M = map(int, input().split())
A = [1]*(M+1)
for i in range(2,M+1):
    A[i] = (A[i-1]*i)%MOD
B = [1]*(M+1)
B[M] = pow(A[M], MOD-2, MOD)
for i in range(M-1,0,-1):
    B[i] = (B[i+1]*(i+1))%MOD
q = M//N
r = M%N
ans = A[M]
for _ in range(N-r):
    ans = (ans*B[q])%MOD
for _ in range(r):
    ans = (ans*B[q+1])%MOD
print(ans)