CodeForces - 1228C Primes and Multiplication（理解题意+唯一分解定理+快速幂）

题目链接：https://vjudge.net/contest/331381#problem/C

分析:
正确理解题意是关键。
prime(x):x的所有素因子
g(x,p)：使p^k尽可能的大，同时满足x%p ^ k=0，最大的p ^ k就是g(x,p)的值
f(x,y)：g(y,x1) * g(y,x2) * …*g(y,xn).(x1…xn为x的素因子)
给定x和n，求f(x,1) * f(x,2) * …f(x,n)的值
正确转化g函数的功能是关键。
g(x,p)可以转化为p * k,而k可以看成p作为质因子在x中出现的次数。
所以问题可以转化为
对于x的每一个素因子，求在1~n这n个数中，出现的总次数sum。对于每一个素因子prime[i]，求prime[i]*sum,再累乘即可。
代码：

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int mod=1e9+7;
#define LL long long
LL x,n;
LL prime[100001];
int k;
LL quickmi(LL x,LL n)
{

    LL vis=1;
    while(n)
    {
        if(n&1)
            vis=vis*x%mod;
        x=x*x%mod;
        n>>=1;
    }
    return vis;

}
void solve1(LL m)
{
    k=0;
    memset(prime,0,sizeof(prime));
    for(int i=2; i*i<=m; i++)/*唯一分解定理*/
    {
        if(m%i==0)
        {
            prime[k++]=i;
            while(m%i==0)
                m/=i;
        }
    }
    if(m>1)
        prime[k++]=m;
}
int main()
{
    while(~scanf("%lld%lld",&x,&n))
    {
        solve1(x);
        LL ans=1;
        for(int i=0; i<k; i++)
        {
            LL s=n,sum=0;
            while(s)/*求次数*/
            {/*对于每一个prime[i]，是作为素因子出现的，只有prime[i]的倍数才能满足情况*/
                s/=prime[i];
                sum+=s;
            }
            ans=(ans*quickmi(prime[i],sum))%mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}