POJ 1845 Sumdiv【数学】

JA_yichao

于 2022-02-26 10:03:45 发布

阅读量229

点赞数

分类专栏：题解数学文章标签： c++ 题解 POJ 数学等比数列

本文链接：https://blog.csdn.net/Jackma_mayichao/article/details/123145721

版权

模幂运算因数分解大整数计算算法效率提升

关键词由CSDN通过智能技术生成

题解同时被 2 个专栏收录

414 篇文章 2 订阅

订阅专栏

数学

40 篇文章 1 订阅

订阅专栏

题目

在这里插入图片描述

思路

在这里插入图片描述

代码

#include<iostream>
#include<cstdio>
using namespace std;
const long long mod=9901;
long long a,b,now,ans=1;
long long ksm(long long x,long long y)
{
	long long js=1;
	while(y)
	 {
	 	if(y&1)
	 	  js=(js*x)%mod;
	 	x=(x*x)%mod;
	 	y>>=1;
	 }
	return js;
}
long long solve(long long p,long long c)
{
	if(!c)
	  return 1;
	if(c&1)
	  return solve(p,(c-1)/2)*(ksm(p,(c+1)/2)+1)%mod;
	else
	  return (solve(p,c-1)+ksm(p,c))%mod;
}
int main()
{
	scanf("%lld%lld",&a,&b);
	for(long long f=2; f*f<=a; f++)
	 {
	 	while(a%f==0)
	 	  a/=f,now++;
	 	if(now)
	 	  ans=ans*solve(f,now*b)%mod;
	    now=0;
	 }
	if(a!=1)
	  ans=ans*solve(a,b)%mod;
	printf("%lld",ans);
	return 0;
}