Number Sequence HDU - 4390 容斥

题意：

首先给出n个数， 然后你可以随机再找n个数， 但必须要让你找的这n个数的乘积等于给出的n个数的乘积， 然后问你可以找出这n个数的方案数并模1e9+7；

思路：

首先我们先把给出的n个数乘起来变成一个数， 再之后根据唯一分解定理我们可以把这个数分解成一些质数， 那么问题不就变成了把这些质数随机组合成n个数的方案数了么， 但要注意此题是有顺序关系的， 比如 4 6 和 6 4 这两个是两种方案；

我们可以按每个分解出来的数来讨论， 假设a[i]是分解出来的某一个数， p[i]是这个数被分解出来的个数， 那么我们分这个数的时候不就相当于将p[i]个一个的数放到n个位置中去么？ 我们很容易就会想到一个叫做隔板法的方法， 我们可以把p[i]扩展成p[i] + n个a[i]， 然后这p[i] + n个a[i] 中产生了p[i] + n - 1个空隙， 我们不就可以在这p[i] + n - 1个空隙里边选出n - 1个空隙， 这样每个区域选出来的数再减去1（反正都是一样的数嘛）不就相当于将p[i]个一样的a[i] 放入了n个位置中了么？

当然， 我们很快就发现了一个问题， 这样放之后， 有些位置没有放数啊！ 那我们怎么办呢？ 容斥。 我们首先求出一个总的答案， 然后再分别求出至少一个不放的， 至少两个不放的，至少三个不放的等等， 然后根据容斥原理不就求出答案了么。

代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>

#define MOD 1000000000 + 7
const int maxed = 1000000 + 10;
typedef long long ll;

int n;
int cnt_pri, cnt, pri[maxed], p[maxed], sum[maxed];
ll c[400 + 2][400 + 2];
bool is_pri[maxed];

int main()
{
    ll slove(int x);
    c[0][0] = 1;
    for (int i = 1; i <= 400; ++i) {
        c[i][0] = 1;
        for (int j = 1; j <= i; ++j)
            c[i][j] = (c[i - 1][j - 1] + c[i - 1][j]) % (MOD);
    }
    for (int i = 2; i < maxed; ++i)
        if (!is_pri[i])
            for (int j = i + i; j < maxed; j += i)
                is_pri[j] = true;
    for (int i = 2; i < maxed; ++i)
        if (!is_pri[i])
            pri[++cnt_pri] = i;
    while (scanf("%d", &n) != EOF) {
        memset(sum, 0, sizeof(sum));
        int a;
        for (int i = 1; i <= n; ++i) {
            scanf("%d", &a);
            for (int j = cnt_pri; j >= 1; --j)
                while (a && !(a % pri[j])) {
                    sum[pri[j]] += 1;
                    a /= pri[j];
                }
        }
        cnt = 0;
        for (int i = 2; i < maxed; ++i)
            if (sum[i])
                p[++cnt] = i;
        ll answer = 0;
        for (int i = 0; i < n; ++i) {
            if (i % 2)
                answer = (answer - c[n][n - i] * slove(n - i) % (MOD) + (MOD)) % (MOD);
            else
                answer = (answer + c[n][n - i] * slove(n - i) % (MOD)) % (MOD);
            //std::cout << answer << std::endl;
        }
        printf("%lld\n", answer);
    }
    return 0;
}

ll slove(int x)
{
    ll ans = 1;
    for (int i = 1; i <= cnt; ++i)
        ans = ans * c[sum[p[i]] + x - 1][x - 1] % (MOD);
    return ans;
}