题目链接:http://codeforces.com/contest/1097
对于,
状态转移方程:。
#include<bits/stdc++.h>
using namespace std;
#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll unsigned long long
#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<ll,ll>
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=2e4+9;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意:给定n和k,
初始时候是n,一次变换后其等概率的变为其因子之一,
问k次变换的期望值。
题目分析:假如已知E(n),当n变为n*p且p与n互质的时候,
分析计算E(n)最后的权重和概率,假设x出现的概率是q,
那么x变为x*p的概率是1/2^k,保持不变的概率是(2^k-1)/2^k,
相当于乘上E(p),不难发现E(n*p)=E(n)*E(p)。
可知E(n)是积性函数,这样我们质因分解后得到质数p和幂次k,
可以分析E(p^k),转移方程很好列,再加上幂次方最高是64,
k的范围是四次方,所以时间复杂度是合适的,
转移方程:E(k,s)=sigma i=0~k E(i,s-1)。
*/
ll n,k;
ll dp[maxn][80];
int main(){
cin>>n>>k;
ll ans=1LL;
mst(dp,0);
for(ll i=2;i*i<=n;i++) if(n%i==0){
int cnt=0;
while(n%i==0) n/=i,cnt++;
dp[0][0]=1LL;rep(j,1,cnt+1) dp[0][j]=dp[0][j-1]*i%mod;
rep(j,1,cnt+1) dp[0][j]=(dp[0][j]+dp[0][j-1])%mod;
rep(j,1,k+1){
dp[j][0]=dp[j-1][0]%mod;
rep(p,1,cnt+1) dp[j][p]=dp[j-1][p]*powmod(p+1,mod-2)%mod;
if(j!=k) rep(p,1,cnt+1) dp[j][p]=(dp[j][p]+dp[j][p-1])%mod;
}
ans=ans*dp[k][cnt]%mod;
}
/// cout<<ans<<endl;
if(n>1){
ll inv=powmod(2LL,mod-2);
inv=powmod(inv,k);
n=(n-1)%mod;
inv=inv*n%mod;
inv=(inv+1)%mod;
ans=ans*inv%mod;
}
cout<<ans<<endl;
return 0;
}