OFFSET
1,1
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
864 is an element of this sequence because 864 = 2^5 * 3^3, so the sum of prime factors is 2 + 2 + 2 + 2 + 2 + 3 + 3 + 3 = 19 which is prime, hence relatively prime to 864. That is the same sum of prime factors as 640 = 2^7 * 5, hence 640 is also a member of this sequence. The sum of prime factors need not be prime for this membership, for example, 2432 = 2^7 * 19 has sum of prime factors 2 + 2 + 2 + 2 + 2 + 2 + 2 + 19 = 33 = 3 * 11, which is composite, yet relatively prime to 2432.
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
Corrected and extended by Ray Chandler, Jul 20 2005
STATUS
approved