#include<stdio.h>
long long int modpow(long long int a, long long int n, long long int p)
{
	long long int res = 1;
	while (n > 0)
	{
		if (n % 2 > 0)
			res = res * a % p;
		a = a * a % p;
		n /= 2;
	}
	return res;
}
long long int conb(long long int n, long long int r, long long int p)
{
	long long int res = 1, i;
	for (i = 1; i <= r; i++)
	{
		res = res * ((n - i + 1) % p) % p;
		res = res * modpow(i, p - 2, p) % p;
	}
	return res;
}
int main()
{
	long long int n, k;
	scanf("%lld %lld", &n, &k);
	long long int a[100003], b[100003], aa, bb;
	long long int i, j, l;
	aa = bb = 0;
	for (i = 1; i * i <= n; i++)
	{
		if (n % i == 0)
		{
			a[aa] = i;
			aa++;
			if (i * i < n)
			{
				b[bb] = n / i;
				bb++;
			}
		}
	}
	while (bb > 0)
	{
		bb--;
		a[aa] = b[bb];
		aa++;
	}
	long long int count = 0;
	for (i = 2, j = n; i * i <= j; i++)
	{
		while (j % i == 0)
		{
			count++;
			j /= i;
		}
	}
	if (j > 1)
		count++;
	long long int p = 1000000007;
	long long int dp[10000][1000];
	for (i = 1; i < aa; i++)
	{
		dp[i][0] = 1;
		for (j = 1; j < count; j++)
			dp[i][j] = 0;
		for (j = 1; j < i; j++)
			if (a[i] % a[j] == 0)
				for (l = 0; l < count - 1; l++)
					dp[i][l + 1] = (dp[i][l + 1] + dp[j][l]) % p;
	}
	if (count > k)
		count = k;
	long long int ans = 1;
	for (i = 1; i < aa; i++)
		for (j = 0; j < count; j++)
			ans = (ans + conb(k, j + 1, p) * dp[i][j] % p) % p;
	printf("%lld\n", ans);
	return 0;
}