#include using namespace std; long long gcd(long long a,long long b){ long long m=min(a,b); long long M=max(a,b); long long temp; while(m>0){ temp=m; m=M%m; M=temp; } return M; } int main(){ int N,mod=1e9+7; long long K,A; map dp; cin>>N>>K; dp[K]=1; for(int i=0;i>A; map next=dp; for(auto j=dp.begin();j!=dp.end();++j){ long long m=gcd(j->first,A); next[j->first/m]=(next[j->first/m]+j->second)%mod; } dp=next; } cout<<(dp[1]-(K==1))<