import sys
input = sys.stdin.readline
import math 

N,K=map(int,input().split())
mod=10**9+7


x=N
L=int(math.sqrt(x))

FACT=dict()

for i in range(2,L+2):
    while x%i==0:
        FACT[i]=FACT.get(i,0)+1
        x=x//i

if x!=1:
    FACT[x]=FACT.get(x,0)+1

def calc(x):
    DP=[0]*(x+1)
    DP[-1]=1

    for i in range(K):
        NDP=[0]*(x+1)

        for j in range(x+1):
            if DP[j]!=0:
                for k in range(j+1):
                    NDP[k]+=DP[j]
                    NDP[k]%=mod
        DP=NDP

    return sum(DP)

ANS=1
for f in FACT:
    ANS*=calc(FACT[f])
    ANS%=mod

print(ANS%mod)