poj1845（唯一分解定理，等比数列求和，约数个数公式）

最新推荐文章于 2021-03-24 01:30:30 发布

林伏案

最新推荐文章于 2021-03-24 01:30:30 发布

阅读量764

点赞数

分类专栏： =====数论===== 素数与整除问题

本文链接：https://blog.csdn.net/qq_29169749/article/details/52810119

版权

=====数论===== 同时被 2 个专栏收录

22 篇文章 0 订阅

订阅专栏

素数与整除问题

11 篇文章 0 订阅

订阅专栏

该博客探讨如何利用唯一分解定理和等比数列求和公式来计算A^B%9901的值。通过质因数分解A，然后使用约数个数公式和模运算性质，结合二分法求解等比数列和，最终得出答案。示例代码中包含实现这一过程的函数。

摘要由CSDN通过智能技术生成

/*
translation:
	输入两个数A，B。求A^B%9901的值
solution:
	唯一分解定理，
	首先将a进行质因数分解，存储在factor二维数组当中。其形式可以描述为：
	p1^q1 * p2^q2 *...* pn^qn。而后要求所有因子的和，可以用到公式：
	约数个数=(1+p1+p1^2+...+p1^q1)...(1+pn+pn^2+pn^qn)。根据这个
	公式，再利用模运算的性质就可以算出答案了。其中计算过程中需要用到等比数列的求和方式。
note:
	*一个数的约数个数的计算公式
	*二分方式求解等比数列的和
date:
	2016.10.13
*/
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>

using namespace std;
typedef long long ll;

const int maxn = 10000 + 1;
const int M = 9901;

ll factor[maxn][2];	//1存放指数,0存放底数
ll a, b;	int id;

void getFac(ll x)
{
	id = 0;
	for(ll i = 2; i*i <= x; i++) {
		if(x % i) continue;
		factor[++id][0] = i;
		factor[id][1] = 0;
		while(x % i == 0) {
			x /= i;
			factor[id][1]++;
		}
	}

	if(x != 1) {
		factor[++id][0] = x;
		factor[id][1] = 1;
	}
}

ll powMod(ll a, ll n)
{
	if(n == 0)	return 1;
	ll x = powMod(a, n/2);
	ll ans = x * x % M;
	if(n % 2 == 1)	ans = ans * a % M;
	return ans;
}

ll getModSum(ll p, ll n)
{
	if(n == 0)    return 1;
    if(p == 0)    return 0;
    if(n & 1)     return (getModSum(p, n/2) * (1 + powMod(p, n/2 + 1))) % M;
    else	      return (getModSum(p, n/2-1) * (1 + powMod( p, n/2 + 1 )) + powMod(p, n/2)) % M;
}

int main()
{
	while(cin >> a >> b) {
		getFac(a);
		ll ans = 1;
		for(int i = 1; i <= id; i++) {
			ans = (ans * getModSum(factor[i][0], factor[i][1]*b)) % M;
		}
		cout << ans << endl;
	}
}