#include using namespace std; struct iofast_t { iofast_t() { ios::sync_with_stdio(false); cin.tie(nullptr); } } iofast; struct uns_t {} uns; template auto vec(Element init, Head arg, Args ...args) { if constexpr (sizeof...(Args) == 0) return std::vector(arg, init); else return std::vector(arg, vec(init, args...)); } template auto vec(uns_t, Head arg, Args ...args) { return vec(Element(), arg, args...); } template > T &chmin(T &l, T r, Compare &&f = less()) { return l = min(l, r, f); } template > T &chmax(T &l, T r, Compare &&f = less()) { return l = max(l, r, f); } #include using mint = atcoder::modint1000000007; int main() { int64_t n, k; cin >> n >> k; auto pf = vec>(uns, 0); for (int i = 2; i <= n / i; ++i) { if (n % i != 0) { continue; } int count = 0; while (n % i == 0) { ++count; n /= i; } pf.push_back({ i, count }); } if (n != 1) { pf.push_back({ n, 1 }); } auto comb = [](int64_t n, int64_t r) { mint ans = 1; for (int i = 0; i < r; ++i) { ans *= n - i; } while (r) { ans /= r--; } return ans; }; mint ans = 1; for (auto [x, y] : pf) { mint sub = 0; for (int i = 0; i <= y; ++i) { sub += comb(k + i - 1, i); } ans *= sub; } cout << ans.val() << endl; }