#pragma GCC optimize("Ofast") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long int ll; typedef unsigned long long ull; mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()); ll myRand(ll B) { return (ull)rng() % B; } constexpr ll mod = 1e9+7; ll mod_pow(ll a,ll b){ a%=mod; if(b==0)return 1; if(b==1)return a; ll res=mod_pow(a,b/2)%mod; res*=res; res%=mod; if(b%2)res*=a; return res%mod; } struct perm{ private: int sz; vector p,invp; public: perm(int n){ sz=n+1; p.resize(sz),invp.resize(sz); p[0]=1; for(int i=1;i<=sz-1;i++){ p[i]=p[i-1]*i%mod; } invp[sz-1]=mod_pow(p[sz-1],mod-2); for(int i=sz-2;i>=0;i--){ invp[i]=invp[i+1]*(i+1)%mod; } } ll comb(ll x,ll y){ if(x> n >> k; // binom(n-1+k,k) auto comb=[&](ll n,ll k)->ll{ ll res = 1; for(ll i=1;i<=k;i++){ res *= (n+1-i); res %= mod; res *= mod_pow(i,mod-2); res %= mod; } return res; }; ll res = 1; for(ll i=2;i*i<=n;i++){ if(n%i==0){ int cnt=0; while(n%i==0){ n/=i; cnt++; } ll sum = 1; for(int j=1;j<=cnt;j++){ sum += comb(k-1+j, j); sum %= mod; } res *= sum; res %= mod; } } if(n > 1){ res *= (1+comb(k, 1)); res %= mod; } cout << res << endl; }