#include #include using namespace std; using namespace atcoder; #define reps(i,s,n) for(int i = s; i < n; i++) #define rep(i,n) reps(i,0,n) #define Rreps(i,n,e) for(int i = n - 1; i >= e; --i) #define Rrep(i,n) Rreps(i,n,0) #define ALL(a) a.begin(), a.end() using ll = long long; using vec = vector; using mat = vector; ll N,M,H,W,Q,K,A,B; string S; using P = pair; using tp = tuple; const ll INF = (1LL<<60); template bool chmin(T &a, const T b){ if(a > b) {a = b; return true;} else return false; } template bool chmax(T &a, const T b){ if(a < b) {a = b; return true;} else return false; } template void my_printv(std::vector v,bool endline = true){ if(!v.empty()){ for(std::size_t i{}; i; const ll MAX_N = ll(4e+6) + 10; vm fact(MAX_N, mint(1)), fact_inv(MAX_N, mint(1)), n_inv(MAX_N, mint(1)); void makefact(){ mint tmp; reps(i,2,MAX_N) fact[i] = fact[i-1] * mint::raw(i); fact_inv[MAX_N - 1] = fact[MAX_N - 1].inv(); Rreps(i, MAX_N - 1, 1){ fact_inv[i] = fact_inv[i + 1] * mint::raw(i + 1); n_inv[i + 1] = fact[i] * fact_inv[i + 1]; } } mint nCm(ll n, ll m){ return fact[n] * fact_inv[n-m] * fact_inv[m]; } mint nCm_inv(ll n, ll m){ return fact[n-m] * fact[m] * fact_inv[n]; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); makefact(); cin>>N>>K; mint res(1); auto calc = [&](int k){ mint sum(0); rep(i, k + 1) { sum += nCm(K + i - 1, i); } res *= sum; }; for(ll i = 2; i * i <= N; ++i){ if(N % i == 0){ int num = 0; while(N % i == 0) { ++num; N /= i; } calc(num); } } if(N != 1) calc(1); cout<