CodeForces - 1228C Primes and Multiplication（题意转化+快速幂）

最新推荐文章于 2024-05-26 21:30:00 发布

zaiyang遇见

最新推荐文章于 2024-05-26 21:30:00 发布

阅读量193

点赞数

分类专栏： # 数论

本文链接：https://blog.csdn.net/lylzsx20172018/article/details/106200439

版权

数论专栏收录该内容

61 篇文章 0 订阅

订阅专栏

题目链接：https://vjudge.net/contest/373416#problem/I
在这里插入图片描述
Input

10 2

Output

Input

20190929 1605

Output

363165664

Input

947 987654321987654321

Output

593574252

翻译：
给定x和n，计算f(x,1) $*$ f(x,2) $*$ …f(x,n)
f(x,y)=g(y, $x_1$ ) $*$ g(y, $x_2$ ) $*$ …g(y, $x_n$ )（ $x_1$ , $x_2$ … $x_n$ 为x的素因子）

g(x,p)=max(p^k1,p^k2,p^k3,p^kn)（x能被p^k整除）

例如：
x=10,n=2
f(10,1) $*$ f(10,2)=g(1,2) $*$ g(1,5) $*$ g(2,2) $*$ g(2,5)

分析：

f(x,1) $*$ f(x,2) $*$ …f(x,n)=
g(1, $x_1$ ) $*$ g(2, $x_1$ ) $*$ g(3, $x_1$ ) $*$ …g(n, $x_1$ ) $*$
g(1, $x_2$ ) $*$ g(2, $x_2$ ) $*$ g(3, $x_2$ ) $*$ …g(n, $x_2$ ) $*$
…
g(1, $x_n$ ) $*$ g(2, $x_n$ ) $*$ g(3, $x_n$ ) $*$ …g(n, $x_n$ )

1~n个数中，能够整除 $x_1$ 的数字一共有n/ $x_1$ 个
整除 $x_1^2$ 的数字一共有n/ $x_1^2$ 个
以此类推。。。。

可以得出在这个n个数中，有n/ $x_1$ 项贡献 $x_1$ ;
有n/ $x_2^2$ 项贡献 $x_2^2$
以此类推。。。。。

考虑：
能够整除 $x_1^2$ 的数字必定能够整除 $x_1$ ,即整除 $x^k$ 的数字必定能够整除 $x^t$ （0<t<k）。
需要去重：
n/ $x_1$ 项贡献的 $x_1$ 全部乘上:quick_mi( $x_1$ ,n/ $x_1$ )(快速幂)
计算n/ $x_1^2$ 项的贡献时，只需要乘以quick_mi( $x_1$ ,n/ $x_1^2$ )

例如：
2 4 6 8 10 ( $x_1$ =2)
第一次乘以 $2^5$
第二次：4和8能整除 $x_2^2$ ,第一次，4和8已经对答案贡献了2，只需要乘以 $2^2$ 即可。

完整代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int mod=1e9+7;
const int N=1e6+10;
LL x,n,prime[N];
int tot;
void solve(LL x)///求x的所有素因子
{
    for(int i=2; i*i<=x; i++)
    {
        if(x%i==0)
        {
            prime[tot++]=i;
            while(x%i==0)
                x/=i;
        }
    }
    if(x>1)
        prime[tot++]=x;

}
LL quick_mi(LL a,LL b)
{
    LL res=1;
    while(b)
    {
        if(b&1)
            res=(res*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return res;
}
LL g(LL n,LL m)
{
    LL t=0;
    for(LL i=m; i<=n; i)
    {
        t+=(n/i);
        n/=i;///1~n中能整除prime[i],prime[i]^2,,,prime[i]^n项的和
    }
    return quick_mi(m,t);
}
int main()
{
    scanf("%lld%lld",&x,&n);
    solve(x);
    LL res=1;
    for(int i=0; i<tot; i++)
    {
        LL sum;
        sum=g(n,prime[i]);
        res=(res*sum)%mod;
    }
    printf("%lld\n",res);
    return 0;
}