#include using namespace std; typedef signed long long ll; #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define FORR2(x,y,arr) for(auto& [x,y]:arr) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) template bool chmax(T &a, const T &b) { if(a bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;} //------------------------------------------------------- int N; int A[202020]; vector>> V; template class UF { public: vector par,rank,cnt,G[um]; UF() {par=rank=vector(um,0); cnt=vector(um,1); for(int i=0;irank[y]) return par[x]=y; rank[x]+=rank[x]==rank[y]; return par[y]=x; } }; UF<2<<20> uf; const int prime_max = 1010101; vector prime; int NP,divp[prime_max]; void cprime() { if(NP) return; for(int i=2;i=i&&j enumpr(ll n) { map V; while(divp[n]>1) { V[divp[n]]++; n/=divp[n]; } if(n>1) V[n]++; return V; } void solve() { int i,j,k,l,r,x,y; string s; cprime(); cin>>N; FOR(i,N) { cin>>x; auto m=enumpr(x); FORR2(a,b,m) { V.push_back({a*b,{i,202020+a}}); } } sort(ALL(V)); reverse(ALL(V)); ll ret=0; FORR2(c,p,V) { if(uf[p.first]==uf[p.second]) { ret+=c; } else { uf(p.first,p.second); } } cout<