Factorization(组合数学,隔板法应用)

题目链接

题意:

找到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+x1n1
对于每个质因子都这样处理后相乘。

#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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值