ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

luogu P1593 因子和

2019-10-29 22:03:12  阅读:250  来源: 互联网

标签:return P1593 int luogu 因子 ans ksm include mod


题目背景

感谢@hehe_54321 提供的3组hack数据

题目描述

输入两个正整数aa和bb,求a^ba
b
的因子和。由于结果太大,只要输出它对9901的余数。

输入格式

仅一行,为两个正整数aa和bb (0\leq a,b \leq 500000000≤a,b≤50000000)。


质因数分解

(1+p1^1+p1^2+...+p1^c1)(1+p2^1+p2^2+...+p2^c2)...*(1+pn^1+pn^2+...+pn^cn)

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int N=3e6+10,M=5e7+1,mod=9901;
#define int long long
int a,b,tmp,ans=1;
inline int ksm(int x,int y){
    int ans=1;
    while(y){
        if(y&1)ans=(ans*x)%mod;
        x=(x*x)%mod;
        y=y>>1;
    }
    return ans;
}
inline int sum(int x,int y){
    if(x==0)return 1;
    if(x%2==1)return sum(x/2,y)*(1+ksm(y,x/2+1))%mod;
    return (sum(x/2-1,y)*(1+ksm(y,x/2+1))%mod + ksm(y,x/2))%mod;
}
signed main(){
    cin>>a>>b;
    int lim=sqrt(a);
    for(int i=2;i<=lim;i++)
    if(!(a%i)){
        tmp=0;
        while(!(a%i)){a/=i;tmp++;};
        ans=ans*sum(tmp*b,i)%mod;
    }
    if(a!=1)ans=ans*sum(b,a)%mod;
    cout<<ans<<endl;
}

标签:return,P1593,int,luogu,因子,ans,ksm,include,mod
来源: https://www.cnblogs.com/naruto-mzx/p/11761822.html

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

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有