【POJ】1845 Sumdiv

SSL_GYX

于 2022-02-26 10:29:56 发布

阅读量233

点赞数

CC 4.0 BY-SA版权

分类专栏：数论/数学文章标签： c++

本文链接：https://blog.csdn.net/SSL_guyixin/article/details/123146094

数论/数学专栏收录该内容

14 篇文章

订阅专栏

本文详细介绍了Sumdiv问题的解决思路，通过分解质因数的方法来计算aba^bab的因数个数，并利用等比数列求和公式进行求和。代码部分展示了使用C++实现的算法，包括ksm函数和sum函数，最终计算出答案。该文章适合对算法和数论感兴趣的编程爱好者阅读。

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

Sumdiv

Link

在这里插入图片描述

解题思路

显然， $a^b$ 的分解质因数就是 $a$ 的分解质因数数量每个乘以 $b$ 。
求和直接用等比数列求和公式即可。

code

#include<iostream>
#include<cstdio>
#include<cmath>
#define int long long
using namespace std;

const int mod=9901;

int a,b,ans;
int x[10000],y[10000],tot;

int ksm(int p,int q)
{
	int t=p,s=1;
	while(q)
	{
		if(q&1) s=(s*t)%mod;
		t=(t*t)%mod;
		q>>=1;
	}
	return s;
}

int sum(int p,int q)
{
	if(q==0) return 1;
	if(q%2) return sum(p,q/2)*(1+ksm(p,q/2+1))%mod;
	else return (sum(p,q-1)+ksm(p,q))%mod;
}

signed main()
{
	cin>>a>>b;
	for(int i=2;i*i<=a;i++)
	{
		if(a%i!=0) continue;
		x[++tot]=i;
		while(a%i==0) y[tot]++,a/=i;
	}
	if(a!=1) x[++tot]=a,y[tot]=1;
	int ans=1;
	for(int i=1;i<=tot;i++)
		ans=(ans*sum(x[i],y[i]*b))%mod;
	cout<<ans<<endl;
}