#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #ifdef LOCAL # include "debug_print.hpp" # define debug(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__) #else # define debug(...) (static_cast(0)) #endif using namespace std; #define rep1(a) for(int i = 0; i < a; i++) #define rep2(i, a) for(int i = 0; i < a; i++) #define rep3(i, a, b) for(int i = a; i < b; i++) #define rep4(i, a, b, c) for(int i = a; i < b; i += c) #define overload4(a, b, c, d, e, ...) e #define rep(...) overload4(__VA_ARGS__, rep4, rep3, rep2, rep1)(__VA_ARGS__) #define rrep1(a) for(int i = (a)-1; i >= 0; i--) #define rrep2(i, a) for(int i = (a)-1; i >= 0; i--) #define rrep3(i, a, b) for(int i = (b)-1; i >= a; i--) #define rrep4(i, a, b, c) for(int i = (b)-1; i >= a; i -= c) #define rrep(...) overload4(__VA_ARGS__, rrep4, rrep3, rrep2, rrep1)(__VA_ARGS__) #define ALL(v) v.begin(), v.end() #define RALL(v) v.rbegin(), v.rend() #define UNIQUE(v) v.erase( unique(v.begin(), v.end()), v.end() ); #define pb push_back using ll = long long; using D = double; using LD = long double; using P = pair; using vi = vector; using vl = vector; template using vc = vector; template using vvc = vector>; template using vvvc = vector>; template using vvvvc = vector>; template using vvvvvc = vector>; #define vv(type, name, h, ...) \ vector> name(h, vector(__VA_ARGS__)) #define vvv(type, name, h, w, ...) \ vector>> name( \ h, vector>(w, vector(__VA_ARGS__))) #define vvvv(type, name, a, b, c, ...) \ vector>>> name( \ a, vector>>( \ b, vector>(c, vector(__VA_ARGS__)))) template using PQ = priority_queue>; template using minPQ = priority_queue, greater>; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b ostream &operator<<(ostream &os, const pair &p) { os << p.first << " " << p.second; return os; } template istream &operator>>(istream &is, pair &p) { is >> p.first >> p.second; return is; } template ostream &operator<<(ostream &os, const vector &v) { int s = (int)v.size(); for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i]; return os; } template istream &operator>>(istream &is, vector &v) { for (auto &x : v) is >> x; return is; } void in() {} template void in(T &t, U &...u) { cin >> t; in(u...); } void out() { cout << "\n"; } template void out(const T &t, const U &...u) { cout << t; if (sizeof...(u)) cout << sep; out(u...); } void outr() {} template void outr(const T &t, const U &...u) { cout << t; outr(u...); } struct eratosthenes{ vector isprime; vector minfactor; vector mobius; eratosthenes(int n) : isprime(n+1,true), minfactor(n+1, -1), mobius(n+1,1){ isprime[1] = false; minfactor[1] = 1; for(int i = 2; i <= n; i++){ if(!isprime[i]) continue; minfactor[i] = i; mobius[i] = -1; for(int j = i+i; j <= n; j += i){ isprime[j] = false; if(minfactor[j] == -1) minfactor[j] = i; if((j/i)%i == 0) mobius[j] = 0; else mobius[j] *= -1; } } } vector> factorize(int n){ vector> res; while(n > 1){ int p = minfactor[n]; int e = 0; while(minfactor[n] == p){ n /= p; e++; } res.push_back({p,e}); } return res; } vector divisor(int n){ vector res; res.pb(1); auto v = factorize(n); for(auto x : v){ int s = res.size(); rep(i,s){ int m = 1; rep(j,x.second){ m *= x.first; res.push_back(res[i]*m); } } } return res; } template void fzt(vector &f){ int n = f.size(); for(int p = 2; p < n; p++){ if(!isprime[p]) continue; for(int q = (n-1)/p; q > 0; q--){ f[q] += f[p*q]; } } } template void fmt(vector &f){ int n = f.size(); for(int p = 2; p < n; p++){ if(!isprime[p]) continue; for(int q = 1; q*p < n; q++){ f[q] -= f[p*q]; } } } template vector conv(const vector &f, const vector &g){ int n = max(f.size(), g.size()); vector nf(n), ng(n), h(n); rep(i,f.size()) nf[i] = f[i]; rep(i,g.size()) ng[i] = g[i]; fzt(nf); fzt(ng); rep(i,n) h[i] = nf[i]*ng[i]; fmt(h); return h; } }; #include "atcoder/dsu.hpp" int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); int n; in(n); eratosthenes si(200001); int c1 = 0; atcoder::dsu d(n+200001); rep(i,n){ int a; in(a); if(a==1) c1++; else{ for(auto [p,_] : si.factorize(a)){ d.merge(i,n+p); } } } int cg = 0; int m = 200000; for(auto &g : d.groups()){ if(g.size()==1) continue; cg++; for(int e : g) if(e >= n) chmin(m, e-n); } int ans = (cg-1+c1)*m; out(ans); }