标签：数论 LL long while ans poj1845 include mod

传送门

求A^B的约数和%9901

这道题算是一道比较综合的题吧！

唯一分解定理那些的就不用说了

朴素求约数和：

^B就指数再乘个B好了

答案就是

然后get到一个新技能

A/B mod C = (A mod (B*C))/ B

简单证明：

然后还要注意的一点就是，计算的时候使用龟速乘，避免爆long long

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <math.h>
#define mod 9901
#define N 50000003
#define LL long long
using namespace std;
LL read()
{
    LL f=1,x=0;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}
    return x*f;
}
int notp[50003],prime[50003];
int tot=0;
void init()
{
    for(int i=2;i<=50000;++i)
    {
        if(!notp[i])prime[++tot]=i;
        for(int j=1;j<=tot&&i*prime[j]<=50000;++j)
        {
            notp[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
} 
LL guicheng(LL a,LL x,LL p)
{
    LL ans=0;
    while(x)
    {
        if(x&1)ans=(ans+a)%p;
        a=(a+a)%p;x>>=1;
    }
    return ans;
}
LL guipow(LL a,LL x,LL p)
{
    LL ans=1;
    while(x)
    {
        if(x&1)ans=guicheng(ans,a,p)%p;
        a=guicheng(a,a,p)%p;x>>=1;
    }
    return ans;
}
int main()
{
    init();
    LL A=read(),B=read();
    LL ans=1;
    LL a=A;
    for(int i=1;prime[i]*prime[i]<=a;++i)
    {
        if(a%prime[i]==0)
        {
            LL k=0;
            while(A%prime[i]==0){k++;A/=prime[i];}
            LL p=mod*(prime[i]-1);
            ans=ans*(guipow(prime[i],k*B+1,p)-1)/(prime[i]-1);
            ans%=mod;
        }
    }
    if(A>1)
    {
        LL p=mod*(A-1);
        ans=ans*(guipow(A,B+1,p)-1)/(A-1);
        ans%=mod;
    }
    printf("%lld\n",ans);
} 

标签：数论,LL,long,while,ans,poj1845,include,mod
来源： https://www.cnblogs.com/yyys-/p/11311469.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。