OFFSET
1,1
COMMENTS
Every positive square (A000290 without 0) is the product of a unique subset of these numbers. The lexicographically earliest (when ordered) minimal set of generators for the positive squares as a group under A059897(.,.); the intersection of A050376 and A000290. - Peter Munn, Aug 25 2019
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Group.
Eric Weisstein's World of Mathematics, Prime Power.
Wikipedia, Generating set of a group.
FORMULA
a(n) ~ n^2 log^2 n. - Charles R Greathouse IV, Oct 19 2015
Sum_{n>=1} 1/a(n) = Sum_{k>=1} P(2^k) = 0.53331724743088069672..., where P is the prime zeta function. - Amiram Eldar, Nov 26 2020
EXAMPLE
3^(2^2) = 81, therefore 81 is a term.
PROG
(PARI) lst(lim)=my(v=List(apply(n->n^2, primes(primepi(sqrtint(lim))))), t); forprime(p=2, (lim+.5)^(1/4), t=p^2; while((t=t^2)<=lim, listput(v, t))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Apr 10 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, May 11 2003
STATUS
approved