hello 2019 D

一开始sb了考虑总的因子疯狂T，做题太少了。。。没意识到会有辣么多因子。。。

神仙说1e9以内的最多的就有800个因子的了。。。

然后我们可以考虑质因子

我觉得已经说得很明白了。。。

唔逆元好像exgcd比费马小定理预处理快？

en...

#include <bits/stdc++.h>
#define mk(a,b) make_pair(a,b)
#define pii pair<int,int>
using namespace std;
typedef long long ll;
const ll mod = 1e9+7;
const int N = 3e7+1;
ll qpow(ll a,ll x){
    ll res = 1;
    while (x){
        if(x&1)
            res=(res*a)%mod;
        a=(a*a)%mod;
        x/=2;
    }
    return res;
}
void exgcd(ll a,ll b,ll& d,ll& x,ll& y) {
    if (!b) {
        d = a;
        x = 1;
        y = 0;
    } else {
        exgcd(b, a % b, d, y, x);
        y -= x * (a / b);
    }
}
ll inv(ll a, ll p) {
    ll d, x, y;
    exgcd(a, p, d, x, y);
    return d == 1 ? (x + p) % p : -1;
}

ll n,k;
ll dp[55][55][10005];
map<ll,int> mp;
map<ll,int> m;
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>k;
    int id = 0;
    for(int i=2;1ll*i*i<=n;i++){
        if(n%i==0) {
            m[i] = id;
            while (n%i==0){
                mp[i]++;
                n/=i;
            }
            dp[id][0][0]=1;
            for(int j=1;j<=mp[i];j++)
                dp[id][j][0]=dp[id][j-1][0]*i%mod;
            id++;
        }
    }
    if(n!=1){
        m[n]=id;
        dp[id][1][0]=n;
        dp[id][0][0]=1;
        mp[n]++;
        id++;
    }
    for(int i=1;i<=k;i++){
        for(auto p:mp){
            int ind = m[p.first];
            ll num = p.first;
            for(int j=0;j<=p.second;j++){//每个次幂
                for(int k=0;k<=j;k++){
                    dp[ind][j][i]+=dp[ind][k][i-1];
                }
                dp[ind][j][i]%=mod;
                dp[ind][j][i]=dp[ind][j][i]*inv(j+1,mod)%mod;
            }
        }
    }
    ll ans = 1;
    for(auto p:mp){
        ans=(ans*dp[m[p.first]][p.second][k]%mod);
    }
    cout<<ans<<endl;
}