hdu 4390

链接：http://acm.hdu.edu.cn/showproblem.php?pid=4390

给你n个数(n<=20),b1,b2,...bn,问你有多少个a1,a2,..an使得ta₁*a₂*...*a_n=b₁*b₂*…*b_n ，并且每个ai要大于1

数据范围1000也可以做

一看就是组合问题，首先，第一步，肯定是先分解质因数，然后得到一个数组a,a[i]表示第i种素数的个数，至于这个素数是什么已经不重要了

然后就相当于依次将每种素数放入n个容器中，记录下每种素数各有多少放法,假设有tot种素数 f(i) 表示第i种素数的放法，那么如果n个数中某些数可以为1，即n个容器某些位置可以不放数，总的方案就是f(1)*f(2)*...*f(tot)，然后再减去所有出现某些容器不放数的情况，可以用容斥原理做，即当前求的总方案数减去有一个容器放0的方案数，加上有两个位置放0的方案数，减去至少有三个位置放0的方案数。。。。。

这样得出的就是答案了。

将m个数放进n个容器中，就相当于有n-1块隔板，总的方案就是C(n-1,n+m-1)

#include<cstdio>
#include<cmath>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int mod = 1000000007;
int C[1000][1000];
vector<int> p;
void init(){
    for(int i=0;i<1000;i++){
        C[i][0]=1;
        C[i][i]=1;
        for(int j=1;j<i;j++){
            C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
        }
    }
}
void make(int num){
    for(int i=2;i*i<=num;i++){
        if(num%i==0){
            while(num%i==0){
                num/=i;
                p.push_back(i);
            }
            if(num==1) break;
        }
    }
    if(num>1) p.push_back(num);
}
int n;
int f(int n,int m){
    return C[n+m-1][n-1];
}
void solve(){
    sort(p.begin(),p.end());

    int tot=0;
    int a[100];
    int sz=p.size();
    for(int i=0,j;i<sz;i=j)  {
        j=i;
        while(j<sz && p[j]==p[i])    {
            j++;
        }
        a[++tot]=j-i;
    }
    __int64 ans=1;
    for(int i=1;i<=tot;i++)     ans*=f(n,a[i]),ans%=mod;
    for(int i=1;i<=n;i++)
    {
        if(i&1)
        {
            __int64 tmp=1;
            tmp*=C[n][i]; tmp%=mod;
            for(int j=1;j<=tot;j++)
            {
                tmp*=f(n-i,a[j]);
                tmp%=mod;
            }
            ans-=tmp;
            ans=(ans%mod+mod)%mod;
        }
        else 
        {
            __int64 tmp=1;
            tmp*=C[n][i];tmp%=mod;
            for(int j=1;j<=tot;j++)
            {
                tmp*=f(n-i,a[j]);
                tmp%=mod;
            }
            ans+=tmp;
            ans=(ans%mod+mod)%mod;
        }
    }
    printf("%I64d\n",ans);
}
int num[25];
int main()
{
    init();
    while(scanf("%d",&n)!=EOF)
    {
        p.clear();
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&num[i]);
            make(num[i]);
        }
        solve();
    }
    return 0;
}