CF1228C Primes and Multiplication（快速幂+埃氏筛法+因式分解）

最新推荐文章于 2024-02-19 17:51:32 发布

原创最新推荐文章于 2024-02-19 17:51:32 发布 · 257 阅读

1 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

6 篇文章

订阅专栏

本文详细解析了如何求解x的素因子在1-n中所有数的贡献问题，通过快速幂和埃氏筛法实现，提供了完整的代码示例，适用于算法竞赛和数学问题求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接:https://vjudge.net/contest/373416#problem/I

思路：求出x的素因子在1-n中所有数的贡献。

Input
10 2
Output
2
Input
20190929 1605
Output
363165664
Input
947 987654321987654321
Output
593574252

代码如下：

#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
#define ll long long
vector<ll>v;
ll book[1000010],su[1000010];
const int mod = 1e9 + 7;
int r=0;
ll qq(ll x,ll n)//快速幂
{
    ll res=1;
    while(n)
    {
        if(n%2)
            res=(res*x)%mod;
        x=(x*x)%mod;
        n/=2;
    }
    return res;
}
void prime()//埃氏筛法
{
    book[2]=0;
    for(ll i=2; i<=1000000; i++)
    {
        if(!book[i])
        {
            su[r++]=i;
            for(ll j=i*i; j<=1000000; j+=i)
                book[j]=1;
        }
    }
}
int main()
{
    prime();
    ll x,n;
    scanf("%lld%lld",&x,&n);
    for(int i=0; i<r; i++)
    {
        if(x%su[i]==0)
        {
            v.push_back(su[i]);//x的因子数
            while(x%su[i]==0)
                x/=su[i];
        }
    }
    if(x>1)
        v.push_back(x);
    ll ans=1;
    for(int i=0; i<v.size(); i++)
    {
        ll p=0;
        ll y=n;
        while(y)//在n中有多少个素数
        {
            p+=y/v[i];
            y/=v[i];
        }
        ll xx=1;
        xx=qq(v[i],p);
        ans=(ans*(xx%mod))%mod;
    }
    printf("%lld\n",ans % mod);
    return 0;
}