Codeforces 615D Multipliers

题目大意

给出一个数$n$的质因数分解：$n=p_1\times p_2\times \cdots \times p_m$，
求$n$的所有因数的乘积，模$10^9+7$。

时间限制

2s

数据范围

$m\le 2\times 10^5$
$p_i\le 2\times 10^5$，且保证$p_i$是质数

题解

$n$非常大，不可能计算出来，更加不可能枚举$n$的因数，
换一种思路，就是考虑每个质因数对答案的贡献。
我们不妨将$n$表示成$p_i^{k_i}$这种形式。
即$n={\Pi }_{i=1}^m{p}_i^{k_i}$
显然$n$的因数个数就有$S={\Pi }_{i=1}^m{k}_i+1$

现在考虑一个质因数$p_i$对答案的贡献，
它自身可以组成的因数有：$p_i^0,p_i^1,\cdots ,p_i^{k_i}$
它们出现的次数是$\frac{S}{k_i+1}$
因此，总的贡献就是$({\Pi }_{j=0}^{k_i}{p}_i^j{)}^{\frac{S}{k_i+1}}$
根据欧拉公式，我们知道指数上的模数为$10^9+6$，显然它不一定存在逆元，所以在这里需要特殊维护一下。

Code

//#pragma GCC optimize (2)
//#pragma G++ optimize (2)
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <vector>
#include <queue>
#define G getchar
#define ll long long
using namespace std;

int read()
{
    char ch;
    for(ch = G();(ch < '0' || ch > '9') && ch != '-';ch = G());
    int n = 0 , w;
    if (ch == '-')
    {
        w = -1;
        ch = G();
    } else w = 1;
    for(;'0' <= ch && ch <= '9';ch = G())n = (n<<1)+(n<<3)+ch-48;
    return n * w;
}

const int N = 200005;
const int mo = 1000000007;

int n , t;
int a[N] , nxt[N];
ll ans , s[N];

ll ksm(ll x , ll y)
{
    ll s = 1;
    for (; y ; )
    {
        if (y & 1) s = s * x % mo;
        x = x * x % mo;
        y = y >> 1;
    }
    return s;
}

int main()
{
    //freopen("b.in","r",stdin);
    //freopen("b.out","w",stdout);
    n = read();
    ans = 1;
    for (int i = 0 ; i < n ; i++)
    {
        t = read();
        a[t]++;
    }
    s[N - 1] = 1;
    for (int i = 1 ; i < N ; i++)
        if (a[i])
        {
            if (s[a[i] + 1] == 0)
            {
                int x = N - 1;
                for (; nxt[x] ;)
                {
                    s[x] = s[x] * (a[i] + 1) % (mo - 1);
                    x = nxt[x];
                }
                nxt[x] = a[i] + 1;
                if (x != N - 1) s[a[i] + 1] = s[x] * x % (mo - 1); else s[a[i] + 1] = s[x];
                s[x] = s[x] * (a[i] + 1) % (mo - 1);
            }
            else
            {
                int x = N - 1;
                for (; x ;)
                {
                    s[x] = s[x] * (a[i] + 1) % (mo - 1);
                    x = nxt[x];
                }
            }
        }

    for (int i = 1 ; i < N ; i++)
        if (a[i]) ans = ans * ksm(i , (ll)a[i] * (a[i] + 1) / 2 % (mo - 1) * s[a[i] + 1] % (mo - 1)) % mo;
    printf("%lld\n", ans);
    return 0;
}