#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; class range {private: struct I{int x;int operator*(){return x;}bool operator!=(I& lhs){return x ostream& operator<<(ostream& os, const pair& p){ return os << "{" << p.first << ", " << p.second << "}"; } template ostream& operator<<(ostream& os, const vector& obj) { os << "{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const set& obj) { os << "set{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template ostream& operator<<(ostream& os, const map& obj) { os << "map{"; for (const auto& e : obj) os << e << ", "; return os << "}"; } template void take(vector& vec, int n) { vec.resize(n); for (int i = 0; i < n; ++i) cin >> vec[i]; } #ifdef ONLINE_JUDGE #define dump(expr) ; #else #define dump(expr) { cerr << "\033[33m#L" << __LINE__ << ": " << expr << "\033[39m" << endl; } #endif using EdgeCost = ll; struct E { int fr, to, eid; EdgeCost cost; E(int fr_, int to_, int eid_, EdgeCost cost_) : fr(fr_), to(to_), eid(eid_), cost(cost_) {} friend ostream &operator<<(ostream &os, const E &e) { os << "(" << e.fr << " -> " << e.to << ")"; return os; } }; using V = vector; struct G { int n, e_count; vector adj; vector rev; G(int n_) : n(n_), e_count(0), adj(n), rev(n) {} void add_directed_edge(int fr, int to, EdgeCost cost = (EdgeCost)1) { int eid = e_count; adj[fr].push_back(E(fr, to, eid, cost)); rev[to].push_back(E(fr, to, eid, cost)); e_count++; } void add_undirected_edge(int fr, int to, EdgeCost cost = (EdgeCost)1) { int eid = e_count; adj[fr].push_back(E(fr, to, eid, cost)); rev[to].push_back(E(fr, to, eid, cost)); adj[to].push_back(E(to, fr, eid, cost)); rev[fr].push_back(E(to, fr, eid, cost)); e_count++; } }; EdgeCost shortest_path(G& g, int src, int dst) { vector vis(g.n, 0); priority_queue> q; q.push(make_pair((EdgeCost)0, src)); while (!q.empty()) { EdgeCost curr_cost = -q.top().first; int idx = q.top().second; q.pop(); if (vis[idx]) continue; vis[idx] = 1; if (idx == dst) { return curr_cost; } for (auto it : g.adj[idx]) { if (vis[it.to]) continue; q.push({-(curr_cost + it.cost), it.to}); } } return (EdgeCost)-1; // impossible } namespace solver { template struct In2 { T1 a; T2 b; friend std::istream& operator>>(std::istream& is, In2& obj) { T1 t1; T2 t2; is >> t1 >> t2; obj = {t1, t2}; return is; } }; template struct In3 { T1 a; T2 b; T3 c; friend std::istream& operator>>(std::istream& is, In3& obj) { T1 t1; T2 t2; T3 t3; is >> t1 >> t2 >> t3; obj = {t1, t2, t3}; return is; } }; ll n, m, p; vector vs; void read() { cin >> n >> m >> p; take(vs, n); } template void vec_unique(vector& v) { ranges::sort(v); auto result = ranges::unique(v); v.erase(result.begin(), result.end()); } template ll vec_index(vector& v, T val) { auto it = ranges::lower_bound(v, val); return it - v.begin(); } using RetType = ll; RetType run() { vector kinds; for (ll d = 1; d * d <= m; ++d) { kinds.push_back(d); kinds.push_back(m / d); } vec_unique(kinds); vec_unique(vs); ll n_kinds = kinds.size(); G g(n_kinds + 1); for (int i : range(n)) { ll a = vs[i]; ll cost = 1; while (a % p == 0) a /= p, cost++; if (a == 1) continue; for (auto upper : kinds) { g.add_directed_edge(vec_index(kinds, upper), vec_index(kinds, upper / a), cost); } } ll maxval = *ranges::max_element(vs); for (auto upper : kinds) { if (upper < maxval) g.add_directed_edge(vec_index(kinds, upper), n_kinds, 1); } return shortest_path(g, n_kinds - 1, n_kinds); } } // namespace template void run(F f) { if constexpr (std::is_same_v) f(); else cout << f() << endl; } int main(int argc, char** argv) { cerr << fixed << setprecision(12); cout << fixed << setprecision(12); int testcase = 1; if (argc > 1) testcase = atoi(argv[1]); while (testcase--) { solver::read(); } run(solver::run); }