OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
From Amiram Eldar, Oct 13 2020: (Start)
a(n) = A005117(n+1)^7.
Sum_{n>=1} 1/a(n) = zeta(7)/zeta(14) - 1. (End)
MATHEMATICA
Select[ Range@34^7, Union[Last /@ FactorInteger@# ] == {7} &] (* Robert G. Wilson v *)
Select[Range[2, 34], SquareFreeQ]^7 (* 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