#include <iostream>
#include <fstream>
#include <cmath>
#include <cassert>

using namespace std;

ifstream in("sumdiv.in");
ofstream out("sumdiv.out");

int const radmax = 8000;
int const nmax = 1001;
int const modulo = 9901;

int ciur[radmax]; //ciur[i] == true <=> composed number
int nfactor, power[nmax];
long long factor[nmax];
int a, b;

void computeciur() {
  ciur[0] = ciur[1] = 1;
  ciur[2] = ciur[3] = 0;
  for (int i = 4; i < radmax; i += 2) {
    ciur[i] = 1;
  }
  for (int i = 5; i < radmax; i += 2) {
    if (ciur[i] == 0) {
      long long j = ((long long) i) * i;
      while (j < radmax) {
        ciur[j] = 1;
        j += i;
      }
    }
  }
}

void decompose(int n) {
  nfactor = 0;
  long lim = sqrt(n);

  if (n % 2 == 0) {
    factor[nfactor] = 2;
    int npow = 0;
    while ((n & 1) == 0) {
      npow++;
      n = n >> 1;
    }
    power[nfactor] = npow;
    nfactor++;
  }

  for (int i = 3; i <= lim; i += 2) {
    if (ciur[i] == 0 && n % i == 0) {
      factor[nfactor] = i % modulo;
      int npow = 0;
      while (n % i == 0) {
        npow++;
        n = n / i;
      }
      power[nfactor] = npow;
      nfactor++;
    }
  }
  if (1 < n) {
    factor[nfactor] = n % modulo;
    power[nfactor] = 1;
    nfactor++;
  }
//debug
//  for (int i = 0; i < nfactor; i++) {
//    cout << factor[i] << "^" << power[i] << endl;
//  }
}

int effpower(int a, int b) {
//  cout << a << " " << b << " -> ";
  if (b == 0)
    return 1;
  else if (b == 1)
    return a;
  else if ((b & 1) == 0) {
    int result = effpower(a, b >> 1);
//    cout << a << " " << b << " " << result << endl;
    return (result * result) % modulo;
  } else {
    int result = effpower(a, b >> 1);
    return (((result * result) % modulo) * a) % modulo;
  }
}

long long numdiv() {
  long long result = 1;
  for (int i = 0; i < nfactor; i++) {
    result *= (power[i] + 1);
  }
  return result;
}

int sumdiv() {
  int sd = 1;
  for (int i = 0; i < nfactor; i++) {
    int modif1 = effpower(factor[i], power[i] + 1) + modulo - 1;
    int modif2 = effpower(factor[i] - 1, modulo - 2);
    sd = (((sd * modif1) % modulo) * modif2) % modulo;
//    cout << effpower(factor[i], power[i] + 1) << " " << modif1 << " " << modif2 << " " << sd << endl;
  }
  return sd;
}

int main() {
//  cout << effpower(2, 4) << endl;
  computeciur();

  in >> a >> b;
  decompose(a);
  for (int i = 0; i < nfactor; i++) {
    power[i] *= b;
  }
  out << sumdiv() << endl;
  return 0;
}