1 /**
2 大意: 计算 a^b 的所有因子的和， 和mod 9901
3 思路； 将a 进行质因子分解，那么所有因子的和为
4 （2^0+ 2^1 + 2^2 +....+ 2^a1）*(3^0 + 3^1+..+ 3^a2)*.....
5 注意： 求模n下a的逆，需要 gcd（a，n） =1，，所以此处不适合
6 用二分求 0-n 次幂的和。。。 又遇到了一次。。。
7 注意中间过程的取模。，，，不然会溢出。。
8 **/
9
10 #include <iostream>
11 #include <cmath>
12 #include <algorithm>
13 #include <cstring>
14 using namespace std;
15 long long ddiv[50000];
16 long long res[50000];
17 int cnt;
18 const int mod = 9901;
19 void get_dr(long long n){
20 memset(res,0,sizeof(res));
21 long long m = (long long )sqrt(n+0.5);
22 cnt =0;
23 for(int i=2;i<=m;i++)if(n%i==0){
24 ddiv[cnt] = i;
25 while(n%i==0){
26 res[cnt]++;
27 n = n/i;
28 }
29 cnt++;
30 }
31 if(n>1){
32 ddiv[cnt] = n;
33 res[cnt]++;
34 cnt++;
35 }
36 }
37
38 long long pow_mod(long long a,long long n){
39 if(n==0)
40 return 1;
41 long long c =1;
42 a =a %mod;
43 while(n){
44 if(n&1)
45 c =c*a%mod;
46 a =a*a%mod;
47 n>>=1;
48 }
49 return c;
50 }
51
52 long long sum (long long p,long long n){
53 if(n==0)
54 return 1;
55 if(n%2){ // 因为是从0开始的。。所以此处才是偶数的情况
56 return (sum(p,n/2)*(1+pow_mod(p,n/2+1)))%mod;
57 }
58 else // 奇数的情况
59 return (sum(p,n/2-1)*(1+pow_mod(p,n/2+1))+pow_mod(p,n/2))%mod;
60 }
61
62 int main()
63 {
64 long long a,b;
65 while(cin>>a>>b){
66 if(b==0){
67 cout<<1<<endl;
68 }else if(a==0){
69 cout<<0<<endl;
70 }else{
71 get_dr(a);
72 //for(int i=0;i<cnt;i++)
73 // cout<<ddiv[i]<<endl;
74 for(int i=0;i<cnt;i++){
75 res[i] *= b;
76 }
77 //for(int i=0;i<cnt;i++)
78 // cout<<res[i]<<endl;
79 long long ans =1;
80 for(int i=0;i<cnt;i++){
81 ans = (ans*(sum(ddiv[i],res[i]))%mod)%mod;
82 ans =ans%mod;
83 }
84 ans = ans%mod;
85 cout<<ans<<endl;
86 }
87 }
88 return 0;
89 }
