题目链接:https://vjudge.net/contest/559971#problem/D
采用隔板法与逆元进行解决
题意:给定一个数为m,和一个范围的数为n,求解在这个数有多少组为这个范围数个的数的乘积为该数(该组数中可以重复);
思路:是直接对该数进行质因数分解,通过对每个数的质因数的个数进行排列,因为是互不干扰的,所以只要求出各个质因数在n中排序的个数的乘积即可;
例如m=12=2^2*3^1,如果n=3,则需要求两个2在n中的排序的个数以及一个3在n中排序的个数,就是先向这组排列中先放2,再放3,把每种放法的乘积后就是最终答案;
由于这组数中可以有相同的,所以在没有填充的位置补放1即可;
那么该如何放置,如何求解这每个数的放置的方案数?
采用隔板法:隔板法:将n 个球放入m个盒子中,盒子可以为空,那么方案数为C(n+m-1,m-1);
来自这知识它不进我脑子啊写的题解中的隔板法的证明:http://t.csdn.cn/Ta1DS
因为为空的地方采用一进行填充,所以与这个隔板法中的思路进行求解方案数;
还需采用逆元,由于需要取模一个较大的数,而会改变原数的大小,在直接除的话,则会出现小数,以及答案可能不准确,则采用逆元将除以一个除数转化为乘以这个数的逆元;
最终AC代码为下:
#include<bits/stdc++.h>
#define ll long long
const int maxn = 2e5+10;
using namespace std;
map<int,int>a;
ll f[300010],chu[300010];
const ll mod=1e9+7;
vector<int>v;
ll shu(ll x,ll y)//快速幂的计算x^y
{
ll ans=1,b=x;
while(y)
{
if(y&1)
ans=ans*b%mod;
b=b*b%mod;
y>>=1;
}
return ans%mod;
}
ll C(ll x,ll y)
{
return (f[x]*chu[y]%mod*chu[x-y]%mod)%mod;//除以某个数,则为乘以所对应除数的逆元即可
}
void init()//初始化乘积数组以及除数的数组
{
int n=3e5;
f[0]=f[1]=1;
chu[0]=chu[1]=1;//需要对零进行初始化为1,否则会在求解C(a,b)时,a=b出现*0的情况
for(int i=2;i<=n;i++)
{
f[i]=(f[i-1]*i)%mod;
chu[i]=(chu[i-1]*shu(i,mod-2))%mod;//除法采用逆元方法进行求解,防止出现小数
}
}
void solve()
{
int r,l,n,k,i,m;
ll ans=1;
cin>>n>>m;
for(i=2;i<=m;i++)//先求解出这个数的每个质因数的个数进行保存
{
while(m%i==0)
{
if(a[i]==0)
v.push_back(i);
a[i]++;
m/=i;
}
}
for(auto x:v)
{
ans=(ans*C(a[x]+n-1,n-1))%mod;//计算方案数
}
cout<<ans%mod<<endl;
}
int main()
{
int t=1;
//cin>>t;
init();
while(t--) solve();
return 0;
}