约数之和【质因数分解+等比数列求和】

题目：

假设现在有两个自然数A和B，S是ABAB的所有约数之和。

请你求出S mod 9901的值是多少。

输入格式

在一行中输入用空格隔开的两个整数A和B。

输出格式

输出一个整数，代表S mod 9901的值。

数据范围

0≤A,B≤5×1070≤A,B≤5×107

输入样例：

2 3

输出样例：

15

注意: A和B不会同时为0。

 

解题报告：由于A，B，都非常地大，咱们不可能直接去求解A^B，所以应该使用一下中学阶段学的幂次的性质，A=（a*b*c)

那么A^B=(a*b*c)^B=a^B*b^B*c^B。然后咱们分别去求解每个素因子的约数和，这里涉及了一下基性函数，每个素因子的求和的乘积=所有因子的加和。在求解每个素因子的和的时候，使用了一下等比数列求和，单独写了一个递归函数去求解的。

ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

const int mod=9901;

int ksc(int a,int b)
{
	int res=0;
	while(b)
	{
		if(b&1)
			res=(res+a)%mod;
		a=(a+a)%mod;
		b>>=1;	
	}
	return res;	
}
int ksm(int a,int b)
{
	a%=mod;
	int res=1;
	while(b)
	{
		if(b&1)
			res=ksc(res,a)%mod;
		a=ksc(a,a)%mod;
		b>>=1;
	}
	return res;
}

int sum(int p,int k)
{
	if(k==0)
		return 1;
	if(k%2==0)
	 	return (p%mod*sum(p,k-1)+1)%mod;
	return ((1+ksm(p,k/2+1))*sum(p,k>>1))%mod;
}
int main()
{
	int a,b;
	std::ios::sync_with_stdio(false);
	cin>>a>>b;
	int res=1;
	for(int i=2;i<=a;i++)
	{
		int s=0;
		while(a%i==0)
		{
			s++;
			a/=i;
		}
		if(s>0)
			res=(res*sum(i,s*b))%mod;
	}
	if(a==0)
		res=0;
	cout<<res<<endl;
}