#include #include #define SELECTER(_1,_2,_3,SELECT,...) SELECT #define rep1(i,n) for(int i=0;i<(int)n;++i) #define rep2(i,a,n) for(int i=(int)a;i<(int)n;++i) #define rep(...) SELECTER(__VA_ARGS__,rep2,rep1)(__VA_ARGS__) #define RSELECTER(_1, _2, _3, RSELECT, ...) RSELECT #define rrep1(i,n) for(int i=(int)n-1;i>=0;--i) #define rrep2(i,a,n) for(int i=(int)n-1;i>=(int)a;--i) #define rrep(...) RSELECTER(__VA_ARGS__, rrep2, rrep1)(__VA_ARGS__) #define all(a) a.begin(),a.end() #define rall(a) a.rbegin(),a.rend() #define fi first #define se second #define PrintR LogOutput #ifdef _DEBUG #define Log(...) LogOutput(__VA_ARGS__) #else #define Log(...) #endif using namespace std; using namespace atcoder; using ll=long long; using pii=pair; using pll=pair; using pdd=pair; using veci=vector; using vecpii=vector>; using vecll=vector; using vecpll=vector>; using vecpdd=vector>; using vecs=vector; using vecb=vector; using vecd=vector; template inline bool chmax(T& a,T b){if(a inline bool chmin(T& a,T b){if(a>b) {a=b;return true;} return false;} constexpr ll INF=numeric_limits::max() / 4; constexpr ll MOD=998244353; using tp=tuple; using tpll=tuple; using mint=modint998244353; using mint10=modint1000000007; template::value && !is_same::value, nullptr_t> = nullptr> T LogOutput(T a, bool b = true){ if(b) cout << a << endl; return a; } template ::value || is_same::value, nullptr_t> = nullptr> auto LogOutput(T a, bool b = true) -> decltype(a.val()){ if(b) cout << a.val() << endl; return a.val(); } template void LogOutput(Args&&... args){ stringstream ss; ((ss << LogOutput(args, false) << ' '), ...); cout << ss.str().substr(0, ss.str().length() - 1) << endl; } template void LogOutput(vector& data){ for(int i = 0 ; i < data.size() ; ++i){ cout << LogOutput(data[i], false) << " \n"[i == data.size() - 1]; } } template void LogOutput(vector>& data){ for(auto d : data) LogOutput(d); } vecpll prime(int n){ vecpll res; int m = n; for(int i = 2; i * i <= m ; ++i){ int cnt = 0; while(n % i == 0) cnt++, n /= i; if(cnt) res.emplace_back(i, cnt); } if(n != 1) res.emplace_back(n, 1); return res; } ll f(ll k){ modint::set_mod(1000000006); return modint(2).pow(k / 2).val(); } void solve(){ int n; ll k;cin>>n>>k; auto t = prime(n); if(k < 60){ rep(j, k){ map mp; rep(i, t.size()){ auto u = prime(t[i].fi + 1); rep(p, u.size()){ mp[u[p].fi] += u[p].se * t[i].se; } } t.clear(); for(auto [l, r] : mp) t.emplace_back(l, r); } mint10 ans = 1; rep(i, t.size()){ ans *= mint10(t[i].fi).pow(t[i].se); } cout< mp; rep(i, t.size()){ auto u = prime(t[i].fi + 1); rep(p, u.size()){ mp[u[p].fi] += u[p].se * t[i].se; } } t.clear(); for(auto [l, r] : mp) t.emplace_back(l, r); } k -= 60; mint10 ans = 1; rep(i, t.size()){ if(k % 2){ ans *= mint10(t[i].fi == 2 ? 3 : 2).pow(f(k) * t[i].se); }else{ ans *= mint10(t[i].fi).pow(f(k) * t[i].se); } } cout<