OFFSET
1,1
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
576 = 2^6 * 3^2 is an element of this sequence because its sum of prime factors is 2 + 2 + 2 + 2 + 2 + 2 + 3 + 3 = 18 = 2 * 3^2 which is a factor of 576 and not relatively prime to 576.
PROG
(PARI) list(lim)=my(v=List()); forprime(p=2, lim\128, forprime(q=2, min(p, lim\64\p), my(pq=p*q); forprime(r=2, min(lim\pq\32, q), my(pqr=pq*r); forprime(s=2, min(lim\pqr\16, r), my(pqrs=pqr*s); forprime(t=2, min(lim\pqrs\8, s), my(pqrst=pqrs*t); forprime(u=2, min(lim\pqrst\4, t), my(pqrstu=pqrst*u); forprime(w=2, min(lim\pqrstu\2, u), my(pqrstuw=pqrstu*w, n); forprime(x=2, min(lim\pqrstuw, w), n=pqrstuw*x; if(gcd(n, p+q+r+s+t+u+w+x)>1, listput(v, n)))))))))); Set(v) \\ Charles R Greathouse IV, Feb 01 2017
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Jul 18 2005
EXTENSIONS
Extended by Ray Chandler, Jul 20 2005
STATUS
approved