#include "bits/stdc++.h"
#include <cassert>
using namespace std;
using ll = long long;
using ull = unsigned long long;
using ld = long double;
using pi = pair<int, int>;
using pl = pair<ll, ll>;
using pd = pair<double, double>;
using pld = pair<ld, ld>;
const int M = 9901, MAX_SQRT_NR = 7071;  // 7071 = floor(sqrt(5e7)).


int main()
{
    bitset<MAX_SQRT_NR + 1> is_prime = move(bitset<MAX_SQRT_NR + 1>{}.set());
    vector<int> primes{2};
    auto eratosthenes = [&]() -> void
    {
        is_prime[0] = is_prime[1] = false;
        for (int j = 4; j <= MAX_SQRT_NR; j += 2)
        {
            is_prime[j] = false;
        }

        int i;
        for (i = 3; i * i <= MAX_SQRT_NR; i += 2)
        {
            if (is_prime[i])
            {
                primes.emplace_back(i);
                for (int j = i * i; j <= MAX_SQRT_NR; j += i)
                {
                    is_prime[j] = false;
                }
            }
        }

        for (i += !(i & 1); i <= MAX_SQRT_NR; i += 2)
        {
            if (is_prime[i])
            {
                primes.emplace_back(i);
            }
        }
    };
    
    auto mod = [](const int a) -> int
    {
        return (a % M + M) % M;
    };
    
    vector<int> mod_inverse(M, 1);
    auto precalc_mod_inverse = [&]() -> void
    {
        for (int i = 2; i < M; ++i)
        {
            mod_inverse[i] = mod(-(M / i) * mod_inverse[M % i]);
        }
    };
    
    eratosthenes();
    precalc_mod_inverse();
    
    ifstream cin("sumdiv.in");
    ofstream cout("sumdiv.out");

    int A, B;
    cin >> A >> B;

    vector<pi> power_prime_fac;
    for (const int prime : primes)
    {
        if (prime * prime > A)
        {
            break;
        }

        int exp{};
        while (A % prime == 0)
        {
            A /= prime;
            ++exp;
        }
        if (exp)
        {
            power_prime_fac.emplace_back(prime, exp * B);
        }
    }
    if (A != 1)
    {
        power_prime_fac.emplace_back(A, B);
    }

    auto fast_exp = [&](int a, int b) -> int
    {
        int res{1};
        while (b)
        {
            if (b & 1)
            {
                res = mod(res * a);
            }
            a = mod(a * a);
            b /= 2;
        }
        return res;
    };

    int res{1};
    for (const auto &prime_pair : power_prime_fac)
    {
       res = mod(res * mod((fast_exp(prime_pair.first, prime_pair.second + 1) - 1) * mod_inverse[(prime_pair.first - 1) % M]));
    }

    cout << res;
}