#include #include using namespace std; using namespace atcoder; typedef int64_t lint; #define rep(i, n) for(int i=0; i; using vvi = vector>; template inline void vin(vector& v) { rep(i, v.size()) cin >> v.at(i); } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } template inline void drop(T x) { cout << x << endl; exit(0); } template void vout(vector v) { rep(i, v.size()) { cout << v.at(i) << ' '; } cout << endl; } constexpr lint LINF = LLONG_MAX/2; using mint = modint1000000007; using vm = vector; vector divisor(lint n) { vector v; for (lint i = 1; i*i <= n; i++) { if (n % i == 0) { v.push_back(i); if (i*i != n) v.push_back(n/i); } } sort(v.begin(), v.end()); return v; } template map coord_comp(set s) { int a = 0; map m; for (auto u : s) m[u] = ++a; return m; } int main() { lint N, M; cin >> N >> M; mint a=0, b=0, c=0; lint x, y, z; vi v = divisor(M); set s; z = v.size(); rep(i, z) s.insert(v[i]); map m = coord_comp(s); vvi w(z); rep(i, z) { rep(j, z) { if (v[i] > LINF/v[j]) continue; if (s.count(v[i]*v[j])) { w[i].pb(j); } } } std::vector dp(N, vm(z)); rep(i, z) { dp[0][i] = 1; } repx(i, 1, N) { rep(j, z) { rep(k, w[j].size()) { dp[i][w[j][k]] += dp[i-1][j]; } } } rep(i, z) a += dp[N-1][i]; std::cout << a.val() << '\n'; }