#include using namespace std; using ll = long long; bool is_prime(int n) { if (n % 2 == 0) return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0) return false; } return n != 1; } vector divisor(int n) { vector res; for (int i = 1; i * i <= n; i++) { if (n % i == 0) { res.emplace_back(i); if (i != n / i) res.emplace_back(n / i); } } return res; } map prime_factor(int n) { map res; for (int i = 2; i * i <= n; i++) { while (n % i == 0) { ++res[i]; n /= i; } } if (n != 1) res[n] = 1; return res; } vector sieve(int n) { vector res; vector is_prime(n + 1, true); is_prime[0] = false; is_prime[1] = false; for (int i = 2; i <= n; i++) { if (!is_prime[i]) continue; res.emplace_back(i); for (int j = 2 * i; j <= n; j += i) is_prime[j] = false; } return res; } struct UnionFind { //各要素が属する集合の代表(根)を管理する //もし、要素xが根であればdata[x]は負の値を取り、-data[x]はxが属する集合の大きさに等しい vector data; UnionFind(int sz) : data(sz, -1) {} bool unite(int x, int y) { x = find(x); y = find(y); bool is_union = (x != y); if (is_union) { if (data[x] > data[y]) swap(x, y); data[x] += data[y]; data[y] = x; } return is_union; } int find(int x) { if (data[x] < 0) { //要素xが根である return x; } else { data[x] = find(data[x]); //data[x]がxの属する集合の根でない場合、根になるよう更新される return data[x]; } } bool same(int x, int y) { return find(x) == find(y); } int size(int x) { return -data[find(x)]; } }; int solve(int n, int p, vector& ok) { UnionFind uf(n); for (int i = 1; i <= n; i++) { for (int j = 1; j < i; j++) { if (__gcd(i, j) != 1) uf.unite(i-1, j-1); } } for (int i = 1; i <= n; i++) { if (uf.same(i-1,p-1)) { ok[i] = true; } } return uf.size(p-1); } int main() { cin.tie(0); ios::sync_with_stdio(false); int n, p; cin >> n >> p; if (p == 1) { cout << 1 << endl; return 0; } vector s = sieve(n); map pf = prime_factor(p); ll hoge = pf.begin()->first; if (hoge != 2 && hoge * 2 <= n) { hoge = 2; } vector ok(n + 1, false); for (auto& p : pf) ok[p.first] = true; ok[hoge] = true; for (int x : s) { if (!ok[x] && x * hoge > n) continue; ok[x] = true; for (int i = 2 * x; i <= n; i += x) { ok[i] = true; } } //vector ok2(n, false); //solve(n, p, ok2); int ans = 0; for (int i = 1; i <= n; i++) { ans += ok[i]; /* if (ok[i] != ok2[i]) { cerr << i << endl; }*/ } cout << ans << endl; return 0; }