OFFSET
1,1
COMMENTS
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures
FORMULA
Sum_{n>=1} 1/a(n) = P(4)*P(6) - P(10) = A085964 * A085966 - P(10) = 0.000320..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020
MATHEMATICA
f[n_]:=Sort[Last/@FactorInteger[n]]=={4, 6}; Select[Range[50000000], f] (*and*) lst={}; Do[Do[If[n!=m, AppendTo[lst, Prime[n]^6*Prime[m]^4]], {n, 50}], {m, 50}]; Take[Union@lst, 50]
PROG
(PARI) list(lim)=my(v=List(), t); forprime(p=2, (lim\16)^(1/6), t=p^6; forprime(q=2, (lim\t)^(1/4), if(p==q, next); listput(v, t*q^4))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 24 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, May 10 2011
STATUS
approved