#include #include // url #define rep(i,a,b) for(int i=a;i=a;i--) #define all(x) (x).begin(),(x).end() #define pb(x) push_back(x); templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b; // const ll mod=998244353; using mint = static_modint<1000000007>; const ll mod=1e9+7; int dx[4]={1,0,-1,0}; int dy[4]={0,1,0,-1}; const string zton="0123456789"; const string atoz="abcdefghijklmnopqrstuvwxyz"; const string ATOZ="ABCDEFGHIJKLMNOPQRSTUVWXYZ"; const ll inf=(1ll<<60); // const int inf=(1<<30); lld dist(lld x1,lld x2,lld y1,lld y2){ lld res=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2); res=sqrt(abs(res)); return res; } lld arg(lld x,lld y){ const lld eps=1e-8; lld res=0; if(abs(x)+abs(y)<=eps)return 0.0; else if(abs(x)<=eps){ if(y>=0.0)return (M_PI/2); else return (M_PI/2+M_PI); } else if(abs(y)<=eps){ if(x>=0.0)return 0.0; else return M_PI; } res=atan2(abs(y),abs(x)); if(x<=0&&y>=0)res=(M_PI-res); else if(x<=0&&y<=0)res+=(M_PI); else if(x>=0&&y<=0)res=(M_PI*2-res); return res; } ll gcd(ll a,ll b){ if(a==0||b==0)return a+b; ll r; r=a%b; if(r==0){ return b; } else{ return gcd(b,r); } } std::ostream &operator<<(std::ostream &dest, __int128_t value) { std::ostream::sentry s(dest); if (s) { __uint128_t tmp = value < 0 ? -value : value; char buffer[128]; char *d = std::end(buffer); do { --d; *d = "0123456789"[tmp % 10]; tmp /= 10; } while (tmp != 0); if (value < 0) { --d; *d = '-'; } int len = std::end(buffer) - d; if (dest.rdbuf()->sputn(d, len) != len) { dest.setstate(std::ios_base::badbit); } } return dest; } __int128 parse128(string &s) { __int128 ret = 0; for (int i = 0; i < s.length(); i++){ if ('0' <= s[i] && s[i] <= '9'){ ret = 10 * ret + s[i] - '0'; } } return ret; } using l3 = __int128; typedef pair P; map memo; mint f(ll N,ll K){ if(memo.count({N,K})){ return memo[{N,K}]; } if(K==0){ return N; } if(N==3){ mint res=2; if(K%2==0){ res++; } else{ res*=2; } ll k=K/2; ll ex=pow_mod(2,k,mod-1); return res.pow(ex); } map cnt; ll N2=N; for(ll i=2;i*i<=N2;i++){ while(N%i==0){ cnt[i]++; N/=i; } } if(N!=1)cnt[N]++; mint res=1; for(auto p:cnt){ // cout << p.first << endl; mint mul=f(p.first+1,K-1); res*=(mul.pow(p.second)); // cout << p.first << " " << mul.val() << endl; // cout << p.first << " " << p.second << endl; } // cout << N << endl; return res; } mint g(ll N,ll K){ rep(i,0,K){ map cnt; for(ll i=2;i*i<=N;i++){ while(N%i==0){ cnt[i]++; N/=i; } } if(N!=1)cnt[N]++; mint res=1; for(auto p:cnt){ mint mul=(p.first+1); mul=mul.pow(p.second); res*=mul; } N=res.val(); // cout << i << " " << N << endl; } return N; } int main(void){ ios::sync_with_stdio(false); cin.tie(nullptr); ll N,K;cin >> N >> K; mint ans=f(N,K); cout << ans.val() << endl; // rep(k,2,11){ // cout << f(N,k).val() << " " << g(N,k).val() << endl; // } // cout << f(N,9).val() << " " << g(N,9).val() << endl; // map cnt; // for(ll i=2;i*i<=K;i++){ // while(K%i==0){ // cnt[i]++; // K/=i; // } // } // if(K!=1)cnt[K]++; }