{-# LANGUAGE BangPatterns #-}
{-# LANGUAGE BinaryLiterals #-}
{-# LANGUAGE MultiWayIf #-}
{-# LANGUAGE NumDecimals #-}
{-# LANGUAGE ScopedTypeVariables #-}
{-# LANGUAGE TupleSections #-}

import Control.Applicative
import Control.Arrow
import Control.Monad
import Control.Monad.ST

import Data.Char
import Data.List
import Data.Maybe

import qualified Data.ByteString.Char8 as B

import Data.Array.ST.Safe
import Data.STRef

import Debug.Trace
import Text.Printf

readInt = readLn :: IO Int
readInts = map ( fst . fromJust . B.readInt ) . B.words <$> B.getLine
readNInts = readInt >>= flip replicateM readInt

which a b f = if f then a else b
mp [ a, b ] = ( a, b )

modBase = 1000000007 :: Integer

a *% b = ( a * b ) `mod` modBase

primeFactorization 1 = []
primeFactorization n = f 2 n
	where
		f p n
			| n == p = [p]
			| n `mod` p == 0 = p : f p ( n `div` p )
			| otherwise = f ( p + 1 ) n

nCr 0 _ = 1
nCr _ 0 = 1
nCr n r = nCr ( n - 1 ) ( r - 1 ) * n `div` r

nHr n r = nCr ( n + r - 1 ) r

main = do
	[ n, m ] <- map fromIntegral <$> readInts
	let degrees = map ( fromIntegral . length ) $ group $ primeFactorization m
	print $ foldl (*%) 1 $ map ( n `nHr` ) degrees