/**
 * code generated by JHelper
 * More info: https://github.com/AlexeyDmitriev/JHelper
 * @author champon
 */

#include <iostream>
#include <fstream>

#include <bits/stdc++.h>
#define ok cout << "ok" << endl;
#define sp << " " <<
using namespace std;
typedef long long ll;
const ll mod = 1e9+7;
const ll inf = 9e18;
template<typename T1, typename T2> ll ceill(T1 a, T2 b){ return (a+b-1)/b; }
void ans(bool b, ostream& out=cout){ out << (b ? "Yes" : "No") << endl; }
void ans2(bool b, ostream& out=cout){ out << (b ? "YES" : "NO") << endl; }
template<typename T> inline bool chmax(T &a, T b){ if(a<b){a=b; return true;} return false; }
template<typename T> inline bool chmin(T &a, T b){ if(b<a){a=b; return true;} return false; }
template<typename A, size_t N, typename T> void Fill(A (&array)[N], const T &val){fill((T*)array, (T*)(array+N), val);}
int dx[] = {1, 0, -1, 0, 1, -1, 1, -1};
int dy[] = {0, 1, 0, -1, 1, -1, -1, 1};

#ifndef INCLUDE_DEBUG_HPP
#define INCLUDE_DEBUG_HPP

#include <vector>
#include <map>
using namespace std;

class debug
{
public:
    typedef long long ll;
    template<typename T> void vd(vector<T> vec);
    template<typename T> void vd2(vector<vector<T>> vec);
    template<typename T> void vpd(vector<pair<int, int>> vec);
    template<typename T> void ad(T *arr, ll n);
    template<typename T, size_t N> void ad2(T (&arr)[N], ll n);
    template<typename T> void md(T mp);
};

#endif //INCLUDE_DEBUG_HPP

template<typename T> void debug::vd(vector<T> vec)
{
    for(auto vv : vec) cout << vv << " ";
    cout << endl;
}

template<typename T> void debug::vd2(vector<vector<T>> vec)
{
    for(auto vv : vec) for(auto vvv : vv) cout << vvv << " ";
    cout << endl;
}

template<typename T> void debug::vpd(vector<pair<int, int>> vec) {
    for(auto const& vv : vec) {
        cout << vv.first << " " << vv.second << endl;
    }
    cout << endl;
}

template<typename T> void debug::ad(T *arr, ll n)
{
    for(int i=0; i<=n; i++) cout << arr[i] << " ";
    cout << endl;
}

template<typename T, size_t N> void debug::ad2(T (&arr)[N], ll n)
{
    for(int i=0; i<=N; i++) for(int j=0; j<=n; j++) cout << arr[i][j] << " ";
    cout << endl;
}

template<typename T> void debug::md(T mp)
{
    cout << "key : value" << endl;
    for(auto const& m : mp) cout << m.first << " : " << m.second << endl;
    cout << endl;
}
debug debug;

#ifndef INCLUDE_FACTOR_HPP
#define INCLUDE_FACTOR_HPP



using namespace std;

class prime {
public:
    typedef long long ll;
    template<typename T> map<ll, ll> factor(T val);
    template<typename T> bool isPrime(T val);
};

#endif //INCLUDE_FACTOR_HPP

template<typename T> map<ll, ll> prime::factor(T val)
{
    vector<T> f;
    for(ll i=2; i*i<=val; i++){
        while(val%i == 0) {
            f.push_back(i);
            val /= i;
            if(val == 0) break;
        }
    }
    if(val > 1) f.push_back(val);
    map<ll, ll> res;
    for(auto const& e : f){
        res[e]++;
    }
    return res;
}

template<typename T> bool prime::isPrime(T val)
{
    for(ll i=2; i*i<=val; i++){
        if(val%i == 0) return false;
    }
    return true;
}
#ifndef DEFINE_COMBINATION_HPP
#define DEFINE_COMBINATION_HPP


using namespace std;

template<typename T>
class Comb
{
    typedef long long ll;
    ll mod = 1e9+7;
  public:
    vector<vector<T>> comb;
    vector<T> fact;
    vector<T> finv;
    vector<T> inv;
    Comb(){};
    Comb(int n) { build(n); }

    void build(int n);
    T result(int n, int k); // return nCk
};

#endif //DEFINE_COMBINATION_HPP

template<typename T>
void Comb<T>::build(int n)
{
    inv.assign(n+1, 0);
    fact.assign(n+1, 0);
    finv.assign(n+1, 0);
    fact[0] = fact[1] = finv[0] = finv[1] = 1;
    inv[1] = 1;
    for(int i=2; i<=n; i++){
        fact[i] = fact[i-1]*i % mod;
        inv[i] = mod - inv[mod%i]*(mod/i)%mod;
        finv[i] = finv[i-1]*inv[i]%mod;
    }
}

template<typename T>
T Comb<T>::result(int n, int k)
{
    if(n < k) return 0;
    if(n < 0 || k < 0) return 0;
    return fact[n]*(finv[k]*finv[n-k]%mod)%mod;
}

// verify
// https://atcoder.jp/contests/abc034/tasks/abc034_c


class DFactorization {
public:
	void solve(std::istream& in, std::ostream& out)
	{
        ll n, m;
        in >> n >> m;
        prime pr;
        map<ll, ll> fact = pr.factor(m);

        ll res = 1;
        Comb<ll> cmb(1e6);
        for(auto const& e : fact) {
            //cout << e.first sp e.second sp cmb.result(e.second+n-1, e.second) << endl;
            res *= cmb.result(e.second + n - 1, e.second);
            res %= mod;
        }

        out << res << endl;
	}
};


int main() {
	DFactorization solver;
	std::istream& in(std::cin);
	std::ostream& out(std::cout);
	solver.solve(in, out);
	return 0;
}


Measured memory_bytes: 23680