CF396A On Number of Decompositions into Multipliers（数论+组合数学）

最新推荐文章于 2018-09-18 20:23:37 发布

Icefox_zhx

最新推荐文章于 2018-09-18 20:23:37 发布

阅读量327

点赞数

分类专栏： codeforces ---------数学--------- 组合数学乘法逆元数论

本文链接：https://blog.csdn.net/Icefox_zhx/article/details/79464770

版权

---------数学--------- 同时被 3 个专栏收录

209 篇文章 0 订阅

订阅专栏

codeforces

122 篇文章 0 订阅

订阅专栏

组合数学

42 篇文章 0 订阅

订阅专栏

将m分解质因数，
$m=\prod_{i=1}^t p_i^{k_i}$
则 $ans=\prod_{i=1}^t C(n+k_i-1,n-1)$ ，相当于对于每一个 $p_i$ ，我要把 $k_i$ 分配到n个位置上去，可以为0，问方案数。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 510
#define mod 1000000007
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,inv[16000],mul[16000],imul[16000];
ll ans=1;map<int,int>cnt;
inline int C(int x,int y){
    if(!y) return 1;
    return (ll)mul[x]*imul[x-y]%mod*imul[y]%mod;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i){
        int x=read(),xx=x;
        for(int j=2;j*j<=xx;++j){
            while(x%j==0) x/=j,cnt[j]++;
        }if(x>1) cnt[x]++;
    }inv[1]=mul[1]=imul[1]=1;
    for(int i=2;i<=15500;++i) inv[i]=(ll)inv[mod%i]*(mod-mod/i)%mod;
    for(int i=2;i<=15500;++i) mul[i]=(ll)mul[i-1]*i%mod;
    for(int i=2;i<=15500;++i) imul[i]=(ll)imul[i-1]*inv[i]%mod; 
    for(map<int,int>::iterator it=cnt.begin();it!=cnt.end();++it)
        ans=(ll)ans*C(it->second+n-1,n-1)%mod;
    printf("%d\n",ans);
    return 0;
}