#include <bits/stdc++.h>
using namespace std;
long long m, x;
vector<pair<long long, int> > v;
vector<pair<long long, int> > u;
map<long long, long long> phi;
vector<long long> dfdfd;
long long product(long long p, long long q, long long mod) {
  if (q == 0) return 0;
  if (q & 1) return (product(p, q - 1, mod) + p) % mod;
  return (product(p, q >> 1, mod) << 1) % mod;
}
long long modpow(long long x, long long y, long long mod) {
  if (y == 0) {
    return 1;
  }
  if (y & 1) return product(modpow(x, y - 1, mod), x, mod);
  long long m = modpow(x, y >> 1, mod);
  return product(m, m, mod);
}
void solve(int i, long long val, long long totient) {
  if (i == v.size()) {
    phi[val] = totient;
    dfdfd.emplace_back(val);
    return;
  }
  long long cur = 1;
  for (int x = 0; x <= v[i].second; x++) {
    long long mult = cur == 1 ? 1 : (cur / v[i].first) * (v[i].first - 1);
    solve(i + 1, val * cur, totient * mult);
    cur *= v[i].first;
  }
}
long long ans;
vector<long long> can;
vector<long long> modulo;
void solver(int x, long long val, int ptr) {
  if (x == v.size()) {
    if (val == 1) return;
    ans += phi[val] / can[ptr];
    return;
  }
  for (int i = 0; i <= v[x].second; i++) {
    while (val > 1 && modulo[ptr] % val != 1) {
      ++ptr;
    }
    solver(x + 1, val, ptr);
    val *= v[x].first;
  }
}
int main(int argc, char const *argv[]) {
  cin >> m >> x;
  long long tmp = m;
  for (long long i = 2; i * i <= m; i++) {
    int cnt = 0;
    while (tmp % i == 0) {
      tmp /= i;
      ++cnt;
    }
    if (cnt) {
      v.emplace_back(i, cnt);
    }
  }
  if (tmp > 1) v.emplace_back(tmp, 1);
  solve(0, 1, 1);
  tmp = phi[m];
  for (long long i = 1; i * i <= tmp; i++) {
    if (tmp % i == 0) {
      can.emplace_back(i);
      if ((i * i) != tmp) {
        can.emplace_back(tmp / i);
      }
    }
  }
  sort(can.begin(), can.end());
  for (auto i : can) {
    modulo.emplace_back(modpow(x, i, m));
  }
  ans = 1;
  solver(0, 1, 0);
  printf("%lld\n", ans);
  return 0;
}
