AtCoder Beginner Contest 110 D - Factorization（隔板法）

题目链接：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;
}