from math import isqrt
n,m = map(int,input().split())
mod = 998244353; s = set()
for i in range(1,isqrt(m)+1):
    if m%i==0: s |= {i,m//i}
s = sorted(s); d = {}
for v in s:
    d[v] = 0
    for u in s: d[v] += v%u==0
    d[v] = pow(d[v],n,mod)
    for u in s:
        if u<v and v%u==0: d[v] -= d[u]
print(d[m]%mod)