#include <bits/stdc++.h>
using namespace std;
using ll=long long;
using pll=pair<ll,ll>;
using tll=tuple<ll,ll,ll>;
const ll INF=(1ll<<60);
#define rep(i,n) for (ll i=0;i<(ll)(n);i++)
#define all(v) v.begin(),v.end()
template<class T> void chmin(T &a,T b){
    if(a>b){
        a=b;
    }
}
template<class T> void chmax(T &a,T b){
    if(a<b){
        a=b;
    }
}
ll modpow(ll a,ll b,ll mod){
    if(mod==1) return 0;
    if(b==0) return 1;
    ll p=a,ret=1;
    for(ll i=0;i<=62;i++){
        if(b&(1ll<<i)){
            ret*=p;
            ret%=mod;
        }
        p=(p%mod)*(p%mod);
        p%=mod;
    }
    return ret;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    const ll mod=1e9+7;
    ll n,p;
    cin >> n >> p;
    vector<ll> f(n+1);
    f[0]=1;
    f[1]=1;
    for(ll i=2;i<=n;i++) f[i]=f[i-1]*i%mod;
    ll cnt=0;
    for(ll i=1;i<=n;i++){
        ll x=i;
        while(x%p==0){
            x/=p;
            cnt++;
            cnt%=mod;
        }
    }
    if(n<p) cout << "0\n";
    else cout << cnt*modpow(f[n],f[n],mod)%mod << '\n';
}