n,k = map(int,input().split())
mod = 10**9+7

### for bigger prime 
fact = [1,1]
finv = [1,1]
inv = [0,1]
 
for i in range(2,2*10**6+5):
    fact.append((fact[-1]*i)%mod)
    inv.append((inv[mod%i]*(mod-mod//i))%mod)
    finv.append((finv[-1]*inv[-1])%mod)
 
def nCr(n,r,mod):
    if r > n:
        return 0
    else: 
        return fact[n]*finv[r]%mod*finv[n-r]%mod


ans = 1

for i in range(2,int(n**0.5)+1):
    if n%i:
        continue
    count = 0
    while n%i == 0:
        n //= i
        count += 1
    ans *= nCr(count+k,k,mod)
    ans %= mod
if n != 1:
    ans *= k+1
    ans %= mod
print(ans)