def fact(num):
	i=2
	ans=dict()
	while i*i <=num:
		while num%i==0:
			num=num//i
			if i in ans:
				ans[i]+=1
			else:
				ans[i]=1
		i+=1
	if num!=1:
		ans[num]=1
	return ans

mod = 10**9+7
n,m = map(int,input().split())

primes = fact(m)
ans = 1

for _,value in primes.items():
    dp = [[0 for _ in range(n+1)] for _ in range(value + 1)]
    dp[0][0] = 1

    for i in range(n):
        s = 0
        for j in range(value+1):
            s += dp[j][i]
            s %= mod
        for j in range(value+1):
            dp[j][i+1] += s
            dp[j][i+1] %= mod
            s -= dp[value-j][i]
    
    cnt = 0
    for j in range(value+1):
        cnt += dp[j][n]
        cnt += mod
    ans *= cnt
    ans %= mod

print(ans)