/*
    给出A,B求A^B内的所有约数之和   mod 9901
    约数之和，有公式 (p1^0 +  + p1^n1)**(pk^0 +  + pk^nk)
    用等比数列求和公式 (p^(n+1)-1)/(p-1)
    这里有1/(p-1)，需要求(p-1)关于mod的逆元
    但是，
    求逆元是有条件的 ax = 1 mod m    存在逆元a^-1 需要 gcd(a,m)  = 1

    如果不满足，可用公式 a/b%m = a%(bm)/b
    这个公式是没条件的，当然b|a
    不过当心bm在ipow那里，最大有bm*bm会溢出
    不过这题，刚好对于(p-1)%m==0的那些数据，不会导致溢出
    
    另外一种做法就是二分了
    对每个1 +p+..+p^n 
    n为奇数 : (1++p^n/2) + p^(n/2+1)(1++p^n/2)
    n为偶数 : (1++p^n/2) + p^n/2(p++p^n/2)
    注意维护p
*/
#include<iostream>
#include<cstring>
#include<map>
#include<algorithm>
#include<stack>
#include<queue>
#include<cmath>
#include<string>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<set>
#include<list>
#include<numeric>
#include<cassert>
#include<ctime>

using namespace std;

const long long MOD = 9901;

long long ipow(long long a, long long n, long long MOD){
    if(n == 0){
        return 1;
    }
    if(n == 1){
        return a % MOD;
    }
    long long ans = ipow(a, n/2, MOD);
    ans = ans * ans % MOD;
    if(n & 1 ){
        ans = ans * a % MOD;
    }
    return ans;
}

//ax + ny = d
long long ExtendEuclid(long long a, long long n, long long &x , long long &y){//ax = b mod n
    if(n == 0){
        x = 1;
        y = 0;
        return a;
    }else{
        //  nx' + a%ny'
        //= nx' + (a - a/n*n)y'
        //= ay' + n(x'-a/ny')
        //=ax + ny
        long long _x, _y, d = ExtendEuclid(n, a % n, _x, _y);
        x = _y;
        y = _x - a/n*_y;
        return d;
    }
}

long long inv(long long a , long long n){//普通逆元求法   ax+ny = 1
    long long x,y;
    ExtendEuclid(a,n,x,y);
    return x % n;
}

long long inv2(long long a, long long n){//x关于素数p的逆元为 x^(p-2) mod p    注意 gcd(x,p) = 1!!!
    return ipow(a,n-2,n);
}
 

long long solve(long long p , long long n){//p^0 +  + p^n
    long long mod = MOD;
    //ax = 1 mod m 。a有逆元的条件是gcd(a,n) = 1 !!!! 否则无解，不能用ExtendEuclid
    //a/b%m = a%(bm)/b  这道题刚好mod = bm，在上面做乘法时不会溢出
    if((p-1) % mod == 0){
        mod *= p-1;
        return (ipow(p,n+1,mod)-1) / (p-1);
    }else{
        //虽然公式a/b%m = a%(bm)/b 在这里也可以用的，但是某些数据，导致ipow那里溢出
        //用逆元就没事
        return (ipow(p,n+1, MOD)-1) * inv(p-1, MOD) % MOD;
    }
}

long long cal(long long a, long long n, long long &p){//a^0+..+a^n  , p = a^n
    if(n == 0){
        p = 1;
        return 1;
    }
    if(n == 1){
        p = a % MOD;
        return (1 + a) % MOD;
    }
    long long ans = cal(a, n/2, p);
    //注意别忘记维护p了
    if(n&1){
        ans = ans * (1 + p*a) % MOD;
        p = p * p * a % MOD;
        return ans;
    }else {
        ans = (ans + (ans- 1) * p % MOD) % MOD;
         p = p * p % MOD;
        return ans;
    }
}

int main()
{
#ifndef ONLINE_JUDGE
    //freopen("in","r",stdin);
#endif

    for(long long a , b ; ~scanf("%lld%lld", &a, &b); ){
        if(a == 0){
            printf("0\n");
            continue;
        }
        long long ans = 1;
        for(long long p = 2 ; p*p <= a ; p++){
            if(a % p == 0){
                long long cnt = 0;
                while(a % p == 0){
                    a /=p;
                    cnt++;
                }
                cnt *= b;
                ans = ans * solve(p ,cnt) % MOD;
            }
        }
        if(a > 1){
            ans = ans * solve(a, b) % MOD;
        }
        printf("%lld\n", (ans+MOD)%MOD);
    }
    return 0;
}