#include #define range(i, n) for(int (i) = 0; (i) < (n); (i)++) #define reversed_range(i, n) for (int (i) = (n) - 1; (i) >= 0; (i)--) using namespace std; template using vec = vector; using lint = long long; using ulint = unsigned long long; struct UnionFind { private: vector parent; vector num_components; public: UnionFind(int n) { num_components.resize(n, 1); for (int i = 0; i < n; i++) { parent.emplace_back(i); } } int root(int x) { if (parent.at(x) == x) return x; parent.at(x) = root(parent.at(x)); return parent.at(x); } bool same(int x, int y) { return root(x) == root(y); } int num_connected(int x) { return num_components.at(x); } void unite(int x, int y) { int root_x = root(x); int root_y = root(y); if (root_x == root_y) return; int rank_x = num_components.at(root_x); int rank_y = num_components.at(root_y); if (rank_x < rank_y) { parent.at(root_x) = root_y; num_components.at(root_y) += rank_x; } else { parent.at(root_y) = root_x; num_components.at(root_x) += rank_y; } } }; template ostream& operator <<(ostream& os, vec v) { os << "["; for (int i = 0; i < v.size() - 1; i++) { os << v.at(i) << ", "; } return os << v.at(v.size() - 1) << "]"; } int gcd(int a, int b) { if (b == 0) return a; return gcd(b, a % b); } vec sieve(int n) { vec is_prime(n + 1, true); is_prime.at(0) = false; is_prime.at(1) = false; for (int i = 2; i <= n; i++) { if (!is_prime.at(i)) continue; for (int j = 2 * i; j <= n; j += i) { is_prime.at(j) = false; } } return is_prime; } int main() { int N, P; cin >> N >> P; vec is_prime = sieve(N); UnionFind uf(N + 1); for (int i = 2; i <= N; i++) { if (is_prime.at(i)) { for (int j = i; j <= N; j += i) { uf.unite(i, j); } } } cout << uf.num_connected(uf.root(P)) << "\n"; }