基础数论（Sumdiv）（快速幂+算术基本定理+约数之和+费马小定理+模运算+逆元）

Sumdiv（好题，知识点多）

这道题涉及到的知识点有点多，而且对定理的理解也要求的比较严格，这一个题让我卡了一下午，哈哈
首先看题目
有两个自然数a和b(a,b≤50000000)

求a的b次方的所有约数之和模9901

输入格式

一行，包含由空格分隔的两个自然数a和b

输出格式

一行，a的b次方的约数和模9901

样例输入

2 3

样例输出

15

样例解释

8的约数是1,2,4,8， 它们的总和是15

15模9901是15

首先先复习一下那些定理们：
算术基本定理： 算术基本定理可表述为：任何一个大于1的自然数 N,如果N不为质数，那么N可以唯一分解成有限个质数的乘积N=P1a1P2a2P3a3…Pnan，这里P1<P2<P3…<Pn均为质数，其中指数ai是正整数。
约数之和: 对于一个大于1正整数n可以分解质因数：n=p1^a1*p2a2*p3^a3*…*pkak,
则由约数个数定理可知n的正约数有(a₁+1)(a₂+1)(a₃+1)…(ak+1)个，
那么n的(a₁+1)(a₂+1)(a₃+1)…(ak+1)个正约数的和为
f(n)=(p1^0+p11+p1^2+…p1a1)(p2^0+p21+p2^2+…p2a2)…(pk^0+pk1+pk^2+…pkak）

费马小定理： 如果 p是一个质数，而整数a不是p的倍数，则有a^（p-1）≡1（mod p）

模运算： 加法减法乘法的模运算都满足分配律，唯独除法不满足

逆元： 在一个群中任取x，若x * y % mod = e,则y是x的逆元，除一个数等于乘以这个数的逆元

思路

首先把 a 给展开，用算数基本定理，然后用约数之和这个结论整出个等比数列来，等比数列的分母如果不是9901的倍数那就求一下逆元再乘上，如果是9901的倍数则（pk^0 + pk1 + pk2 + … + pks）% 9901 = (1 + 1 + 1 + … + 1) = s + 1,因为每项用二项式定理展开再mod的话除了那个1都成0了，所以就只剩下那个1了，所以就特判一下就行了，这个特判没了的话会一直wa，wa的心态崩了都，具体看代码把

#include<iostream>
#include<cstdio>
using namespace std;
typedef unsigned long long LL;
typedef pair<int,int> PII;
const int N = 10005,mod = 9901;
LL ksm(LL x,LL n){
	LL ans = 1;
	x %= mod;
	while(n){
		if(n & 1) ans = ans * x % mod;
		n >>= 1;
		x = x * x % mod;
	}
	return ans;
}
int main(){
	LL a,b,res = 1;
	cin>>a>>b;
	for(int i=2;i*i<=a;i++){
		if(a % i == 0){
			LL cnt = 0;
			while(a % i == 0){
				a /= i;
				cnt ++;
			}
			cnt *= b;
			LL s = i - 1;
			if(s % mod == 0){//特判
				res = res * (cnt+1) % mod;
				continue;
			}
			s = ksm(s,mod-2);
			res = res * (ksm(i,cnt+1)-1 + mod) % mod * s % mod;
		}
	}
	if(a > 1){
		LL s = a - 1;
		if(s % mod == 0){
			res = res * (b + 1) % mod;
		}
		else{
			s = ksm(s,mod-2);
			res = res * (ksm(a,b+1)-1 + mod) % mod * s % mod;
		}
	}
	cout<<res<<endl;
	return 0;
}