/** * author: shu8Cream * created: 17.02.2023 21:59:02 **/ #include using namespace std; #define overload3(a,b,c,d,...) d #define rep1(i,n) for (int i=0; i<(n); i++) #define rep2(i,a,n) for (int i=(a); i<(n); i++) #define rep(...) overload3(__VA_ARGS__, rep2, rep1)(__VA_ARGS__) #define rrep1(i,n) for (int i=(n-1); i>=0; i--) #define rrep2(i,a,n) for (int i=(n-1); i>=(a); i--) #define rrep(...) overload3(__VA_ARGS__, rrep2, rrep1)(__VA_ARGS__) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() #define sz(x) int((x).size()) #define pcnt __builtin_popcountll using ll = long long; using P = pair; template using vc = vector; template using vv = vc>; using vi = vc; using vvi = vv; const int INF = 1e9; const ll LINF = 8e18; templateistream& operator>>(istream&i,vc&v){rep(j,sz(v))i>>v[j];return i;} templatestring join(const T&v,const string& d=""){stringstream s;rep(i,sz(v))(i?s<ostream& operator<<(ostream&o,const vc&v){if(sz(v))o<istream& operator>>(istream&i,pair&v){return i>>v.first>>v.second;} templateostream& operator<<(ostream&o,const pair&v){return o< inline bool chmax(T& a, T b) {if(a inline bool chmin(T& a, T b) {if(a>b) { a=b;return true; } return false;} // stable sort template vc argsort(const vc &A) { vc ids(sz(A)); iota(all(ids), 0); sort(all(ids), [&](int i, int j) { return (A[i] == A[j] ? i < j : A[i] < A[j]); }); return ids; } // A[I[0]], A[I[1]], ... template vc rearrange(const vc &A, const vc &I) { vc B(sz(I)); rep(i, sz(I)) B[i] = A[I[i]]; return B; } template vc cumsum(vc &A, int off = 1) { vc B(sz(A) + 1); rep(i, sz(A)) B[i + 1] = B[i] + A[i]; if (off == 0) B.erase(B.begin()); return B; } template string to_string(T s); template string to_string(pair p); string to_string(char c) { return string(1, c); } string to_string(string s) { return s; } string to_string(const char s[]) { return string(s); } template string to_string(T v) { if (v.empty()) return "{}"; string ret = "{"; for (auto x : v) ret += to_string(x) + ","; ret.back() = '}'; return ret; } template string to_string(pair p) { return "{" + to_string(p.first) + ":" + to_string(p.second) + "}"; } void debug_out() { cout << endl; } template void debug_out(Head H, Tail... T) { cout << to_string(H) << " "; debug_out(T...); } #ifdef _DEBUG #define debug(...) debug_out(__VA_ARGS__) #else #define debug(...) #endif struct Sieve{ int n; vector f,primes; Sieve(int n=1):n(n+1),f(n+1){ f[0] = f[1] = -1; for(ll i=2; i<=n; ++i){ if(f[i]) continue; primes.push_back(i); f[i]=i; for(ll j=i*i; j<=n; j+=i){ if(!f[j]) f[j]=i; } } } bool isPrime(int x){ return f[x]==x;} vector factorList(int x){ vector res; while(x!=1){ res.push_back(f[x]); x/=f[x]; } return res; } vector

factor(int x){ vector fl = factorList(x); if(fl.size()==0) return {}; vector

res(1,P(fl[0], 0)); for(int p : fl){ if(res.back().first == p){ res.back().second++; }else{ res.emplace_back(p,1); } } return res; } // mの約数を列挙、返り値は約数の個数 int divisor(ll m, vi& a){ ll cnt = 1; a[0] = 1; while(m!=1){ ll num=1, i=f[m], temp=i; while(m%i==0){m/=i;num++;} rep(k,1,num){ rep(j,0,cnt)a[k*cnt+j]=a[j]*temp; temp*=i; } cnt*=num; } sort(a.begin(),a.begin()+cnt); return cnt; } }sieve(500100); int main() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); ll n; cin >> n; vi a(n); cin >> a; vvi mp(100100); rep(i,n) mp[a[i]].push_back(i); vi dp(n,1); rrep(i,n){ vi div(5000); int cnt = sieve.divisor(a[i],div); rep(j,cnt)for(auto p:mp[div[j]]){ if(p