#include "bits/stdc++.h"

#include <unordered_set>
#include <unordered_map>

//#include <iostream>//POJ
//#include <vector>
//#include <string>
//#include <iomanip>
//#include <math.h>
//#include <algorithm>
//#include <cstring>
using namespace std;
#define l0 (lint)0
#define setc cin.tie(0);cout.tie(0);ios::sync_with_stdio(0)
#define dd(n) cout<<fixed<<setprecision(n)
//#define repi(i,a,b) for(int i=(a), i##_len=(b); i<i##_len; ++i)
#define repi(i,a,b) for(int i=(a);i<b;++i)
#define rep(i,n) repi(i,0,n)
#define repi_(i,a,b) for(int i=(a), i##_len=(b); i<=i##_len; ++i)
#define rep_(i,n) repi_(i,0,n)
#define repir(i,a,b) for(int i=(a)-1, i##_first=(b); i>=i##_first; --i)
#define repr(i,n) repir(i,n,0)
#define repir_(i,a,b) for(int i=(a), i##_first=(b); i>=i##_first; --i)
#define repr_(i,n) repir_(i,n,0)
#define prep(i,n) for(int i=0;i*i<n;++i)
#define prepi(i,a,b) for(lint i=a;i*i<=b;++i)
#define all(x) (x).begin(),(x).end()
#define vsort(v) sort((v).begin(), (v).end())
#define vcsort(v,c) sort((v).begin(),(v).end(),c)
#define vrev(v) reverse((v).begin(), (v).end())
#define fnd(v,x) ((v).find(x)!=string::npos)
#define mod(a,b) (a<0?a%b+abs(b):a%b)
#define cl(a,b) ((a+b-1)/b)

//temlyate<typename T> T vpop(vector<T> *v) {
//	T d = v->back();
//	v->pop_back();
//	return d;
//}
//char vpop(string *v) {
//	char d = v->back();
//	v->pop_back();
//	return d;
//}
#define siz(v) ((int)(v).size())
#define ers(v, n) (v).erase((v).begin() + n)
#define vmin(v) *min_element(v.begin(), v.end())
#define vmax(v) *max_element(v.begin(), v.end())
#define contain(q) !q.empty()
#define cont(q) !q.empty()
//#define qpop(q, a, b) a=q.back().first;b=q.back().second;q.pop()
//#define pqpop(p, a, b) a=q.top().first;b=q.top().second;q.pop()
#define el "\n"
#define sp " "
#define pi 3.14159265358979
#define co(x) cout<<x<<el
#define coc(c, a, b) co(((c)?a:b))
#define cim(x) {cin>>x;--x;}
#define cim2(a,b) {cin>>a>>b;--a;--b;}
#define cosp(x) cout<<(x)<<' '
#define YES(c) coc(c,"YES", "NO")
#define Yes(c) coc(c,"Yes", "No")
#define yes(c) coc(c,"yes", "no")
#define POSSIBLE(c) coc(c, "POSSIBLE", "IMPOSSIBLE")
#define Possible(c) coc(c, "Possible", "Impossible")
#define possible(c) coc(c, "possible", "impossible")
#define inf INT_MAX
#define linf INT64_MAX*0.99
#define wildcard(T) numeric_limits<T>::min()
#define noway() {co(-1);gc();return 0;}

#define pb push_back
#define pq priority_queue
#define np next_permutation

#ifdef tqktmp_2
#define gc() getchar(); getchar()
#else
#define gc() 1
#endif
const unsigned int bf0 = (1 << 0);
const unsigned int bf1 = (1 << 1);
const unsigned int bf2 = (1 << 2);
const unsigned int bf3 = (1 << 3);
const unsigned int bf4 = (1 << 4);
const unsigned int bf5 = (1 << 5);
const unsigned int bf6 = (1 << 6);
const unsigned int bf7 = (1 << 7);

//#define lint lint
//typedef lint lint;
typedef int_fast64_t lint;
typedef vector<int> IV; typedef vector<string> SV;
typedef vector<lint> LIV;
typedef vector<vector<int> > IVV;
typedef pair<int,int> P; typedef pair<lint,lint> LP;
typedef vector<P> PV; typedef vector<LP> LPV;

const int
dx8[8] ={0,1,1,1,0,-1,-1,-1},
dy8[8] ={1,1,0,-1,-1,-1,0,1},
dx9[9] ={0,1,1,1,0,0,-1,-1,-1},
dy9[9] ={1,1,0,-1,0,-1,-1,0,1},
dx6[2][6]={{0,1,0,-1,-1,-1},{1,1,1,0,-1,0}},
dy6[6]={-1,0,1,1,0,-1},
dx5[5] ={0,1,0,0,-1},
dy5[5] ={1,0,0,-1,0},
dx4[4] ={0,1,0,-1},
dy4[4] ={1,0,-1,0};
//temlyate<class... A> void cim_(A... args) {
//	for (A *i : initializer_list<A*>{ args... }) {
//		cin>>*i;--*i;
//	}
//	return;
//}
lint input[5];

template<class T>inline bool maxi(T &a,const T &b){ if(a < b){ a = b; return 1; } return 0; }
template<class T>inline bool mini(T &a,const T &b){ if(b < a){ a = b; return 1; } return 0; }

//kokomade<C-m><C-h>

#define MAX_N 200020
const int md = 1e9+7;

struct AB{ lint a,b; };//a,b
struct ABC{ lint a,b,c; };//a,b,c
struct C{ lint x,y; };//coords
struct SEC{ lint l,r; };//section

#define MAXN 100010

lint powmod_h;
lint powmod(lint a,int p){
	if(!p)return 1;
	if(p%2)return powmod(a,p-1)* a % md;
	else{
		powmod_h = powmod(a,p / 2);
		return powmod_h*powmod_h % md;
	}
}

lint ncr_ret,ncr_div;
lint ncr(lint n,lint r){
	ncr_ret=1;ncr_div=1;
	if(r>n/2)return ncr(n,n-r);
	rep(i,r){
		ncr_ret=ncr_ret*(n-i)%md;
		ncr_div=ncr_div*(i+1)%md;
	}ncr_ret*=powmod(ncr_div,md-2);
	return ncr_ret%md;
}

int main(){
	setc;

	lint n,m;cin>>n>>m;
	lint ans=1;
	lint plf=m,pp;
	prepi(i,2,m){
		if(!(plf%i)){
			pp=1;plf/=i;
			while(plf%i==0){
				++pp;plf/=i;
			}
			//ここで、因数i^ppがわかる i is prime
			//ma[i]=pp;//処理
			ans*=ncr(pp+n-1,n-1);
			ans%=md;
		}
	}
	if(plf!=1){
		//ここで、因数plf^1がわかる plf is prime
		ans*=ncr(n,n-1);
		ans%=md;
	}co(ans);
	gc();
}

Measured memory_bytes: 256