CodeForces 396 A.On Number of Decompositions into Multipliers（组合数学+数论）

Description

给出一个正整数$m$，问满足$\prod\limits_{i=1}^na_i=m$的长度为$n$的序列$a_i$个数

Input

两个正整数$n,m(1\le n\le 500,1\le m\le 10^9)$

Output

输出满足条件的序列个数，结果模$10^9+7$

Sample Input

1

15

Sample Output

1

Solution

把$n$质因数分解$n=\prod\limits_{i=1}^kp_i^{a_i}$，那么问题变为：对于每个质因子$p_i$，要把$a_i$个质因子分配到$n$个不同的数字里，一些数字可以没有该因子，由插板法，方案数为$\prod\limits_{i=1}^kC_{n-1+a_i}^{n-1}$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=505;
#define mod 1000000007
int inv[maxn],n;
map<int,int>M;
map<int,int>::iterator it;
void init(int n=500)
{
    inv[1]=1;
    for(int i=2;i<=n;i++)inv[i]=mod-(ll)(mod/i)*inv[mod%i]%mod;
}
int C(int n,int m)
{
    int ans=1;
    for(int i=1;i<=m;i++)ans=(ll)ans*(n-m+i)%mod*inv[i]%mod;
    return ans;
}
void Solve(int x)
{
    for(int i=2;i*i<=x;i++)
        while(x%i==0)
            x/=i,M[i]++;
    if(x>1)M[x]++;
}
int main()
{
    init();
    while(~scanf("%d",&n))
    {
        M.clear();
        for(int i=1;i<=n;i++)
        {
            int a;
            scanf("%d",&a);
            Solve(a);
        }   
        int ans=1;
        for(it=M.begin();it!=M.end();it++)
            ans=(ll)ans*C(n+it->second-1,n-1)%mod;
        printf("%d\n",ans);
    }
    return 0;
}