OFFSET
1,1
LINKS
Nathaniel Johnston, Table of n, a(n) for n = 1..5000
FORMULA
a(n) = A005117(n+1)^6. - Nathaniel Johnston, Jun 21 2011
Sum_{n>=1} 1/a(n) = zeta(6)/zeta(12) - 1 = A269404 - 1. - Amiram Eldar, Oct 13 2020
MAPLE
for n from 2 to 100 do if(numtheory[issqrfree](n))then printf("%d, ", n^6): fi: od: # Nathaniel Johnston, Jun 21 2011
MATHEMATICA
Select[ Range@37^6, Union[Last /@ FactorInteger@# ] == {6} &] (* Robert G. Wilson v *)
Select[Range[2, 37], SquareFreeQ]^6 (* Amiram Eldar, Oct 13 2020 *)
PROG
(PARI) allpwrfact(n, p) = \All prime factors are raised to the power p { local(x, j, ln, y, flag); for(x=4, n, y=Vec(factor(x)); ln = length(y[1]); flag=0; for(j=1, ln, if(y[2][j]==p, flag++); ); if(flag==ln, print1(x", ")); ) }
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Jan 25 2006
EXTENSIONS
More terms from Robert G. Wilson v, Jan 26 2006
STATUS
approved