其他分享
首页 > 其他分享> > luogu P1593 因子和

luogu P1593 因子和

作者:互联网

题目背景

感谢@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