A374291 Squares of powerful numbers.

1, 16, 64, 81, 256, 625, 729, 1024, 1296, 2401, 4096, 5184, 6561, 10000, 11664, 14641, 15625, 16384, 20736, 28561, 38416, 40000, 46656, 50625, 59049, 65536, 82944, 83521, 104976, 117649, 130321, 153664, 160000, 186624, 194481, 234256, 250000, 262144, 279841, 331776

Offset

1,2

Comments

First differs from A340588 at n = 12.

4-full (or 3-full) squares.

Numbers whose exponents in their prime factorization are all even numbers >= 4.

This sequence is closed under multiplication.

The sequence {A000290(n)*A078615(A000290(n)), n>=1} is a permutation of this sequence, and the sequence {a(n)/A078615(a(n)), n>=1} is a permutation of {A000290(n), n>=1}.

The sequence {A335988(n)*A007947(A335988(n)), n>=1} is a permutation of this sequence, and the sequence {a(n)/A007947(a(n)), n>=1} is a permutation of A335988.

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Index entries for sequences related to powerful numbers.

Formula

a(n) = A000290(A001694(n)) = A001694(n)^2.

Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/(p^2*(p^2-1))) = zeta(4)*zeta(6)/zeta(12) = 15015/(1382*Pi^2) = 1.10082313486953808844... .

Sum_{n>=1} 1/a(n)^s = Product_{p prime} (1 + 1/(p^(2*s)*(p^(2*s)-1))) = zeta(4*s)*zeta(6*s)/zeta(12*s), for s > 1/4.

Mathematica

powQ[n_] := n==1 || AllTrue[FactorInteger[n][[;; , 2]], # > 1 &]; Select[Range[600], powQ]^2

Prog

(PARI) is(k) = issquare(k) && ispowerful(sqrtint(k));

Crossrefs

Intersection of A000290 and A036967 (or A036966).

Intersection of A000290 and A337050.

Subsequence of A322449.

Cf. A000290, A001694, A007947, A078615, A335988, A340588.

Sequence in context: A370787 A322449 A117453 * A340588 A352475 A305242

Adjacent sequences: A374288 A374289 A374290 * A374292 A374293 A374294

Keyword

nonn,easy

Author

Amiram Eldar, Jul 02 2024

Status

approved