题意:
找到n个数的成绩等于m的总排列方案数(mod 1e9+7)。
分析:
由乘积我们可以想到将m质因数分解得到质因数的种类数和每个质因数的个数。
只要n个数种包含这全部的质因子就行。便转换成了排列问题。
当我们通过质因数分解得到每个质因子个数之后,对于每个质因子x有cnt个:
将这cnt个x放到n个位置上,某个位置可能为0个x,
由此想到隔板法,但是隔板法是处理至少为1的情况,
这时候我们假设先给n个位置上都先放置一个数x,
然后就变成了n+x个数放到n个位置上,每个位置至少放一个的总方案数;
可以写成:
C
n
+
x
−
1
n
−
1
C_{n+x-1}^{n-1}
Cn+x−1n−1
对于每个质因子都这样处理后相乘。
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define sqr(x) (x)*(x)
const int mod = 1e9+7;
const int maxn = 1e6+10;
double dp[maxn];
ll f[maxn];
ll qpow(ll a, ll b){
ll res = 1;
while(b){
if(b&1) res = res * a % mod;
b >>= 1;
a = a * a % mod;
}
return res;
}
void init(){
f[0] = 1;
for(int i=1; i<=1e6+10; i++){
f[i] = f[i-1] * i % mod;
}
}
ll cal(ll n, ll m){
if(n<m) return 0;
return 1ll*f[n]*qpow(f[m], mod-2)%mod*qpow(f[n-m], mod-2)%mod;
}
int main(){
init();
ll n, m;
while(cin >> n >> m){
ll res = 1;
for( ll i=2; i*i<=m; i++){
if(m % i == 0){
ll cnt = 0;
while(m % i == 0){
m /= i;
cnt++;
}
res = (res * cal(n+cnt-1, n-1)%mod + mod) % mod;
}
}
if(m != 1) res = res * n % mod;
cout << res << endl;
}
return 0;
}