CF Round 232 C. On Number of Decompositions into Multipliers 素因子分解,组合数

题意:给出长度为n的序列a, 其中m=a[1]*a[2]..*a[n]. 

n<=1e5,a[i]<=1e9 问将m拆分成n个数的乘积的方案数?

将每个a[i]素因子分解. 从而得到m的每个素因子幂.

初始n个位置都是1. 第i个素因子幂为val,现在将第i个素因子放入n个位置中的方案为C(val+n-1,n-1).总方案累乘即可.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5,mod=1e9+7;
int n,x,vis[N],pr[N],pn=0;
map<int,ll> mp;
ll f[N];
void init()
{
	for(int i=2;i<N;i++){
		if(vis[i])	continue;
		pr[++pn]=i;
		for(int j=i+i;j<N;j+=i)
			vis[j]=true;
	}
	f[0]=1;
	for(ll i=1;i<N;i++)	f[i]=(f[i-1]*i)%mod;
}
ll powmod(ll x,ll n){ll s=1;for(;n;n>>=1,x=(x*x)%mod) if(n&1)s=(s*x)%mod; return s;}
ll C(ll n,ll m)
{
	if(n<m)	return 0;
	ll a=f[n],b=(f[n-m]*f[m])%mod;
	return (a*powmod(b,mod-2))%mod;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	init();
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>x;
		for(int j=1;j<=pn&&x!=1;j++)
			while(x%pr[j]==0)
				mp[pr[j]]++,x/=pr[j];
		if(x!=1)	mp[x]++;
	}
	ll res=1;
	for(auto it=mp.begin();it!=mp.end();it++){
		ll val=it->second;
		ll cnt=C(val+n-1,n-1);
		res=(res*cnt)%mod;
	//	cout<<(it->first)<<' '<<val<<' '<<cnt<<'\n';
	}
	cout<<res<<'\n';
	return 0;
}