A147661 a(n) = squarefree part of n^n.

1, 1, 3, 1, 5, 1, 7, 1, 1, 1, 11, 1, 13, 1, 15, 1, 17, 1, 19, 1, 21, 1, 23, 1, 1, 1, 3, 1, 29, 1, 31, 1, 33, 1, 35, 1, 37, 1, 39, 1, 41, 1, 43, 1, 5, 1, 47, 1, 1, 1, 51, 1, 53, 1, 55, 1, 57, 1, 59, 1, 61, 1, 7, 1, 65, 1, 67, 1, 69, 1, 71, 1, 73, 1, 3, 1, 77, 1, 79, 1, 1, 1, 83, 1, 85, 1, 87, 1, 89

Offset: 1

Author

Artur Jasinski, Nov 09 2008

Keywords

• nonn

Links

• Alois P. Heinz, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A000312, A007913, A055204, A147662.

Programs

Mathematica

data = Table[Sqrt[n^n], {n, 1, 200}]; sp = data /. Sqrt[] -> 1; sfp = data/sp /. Sqrt[x] -> x

PARI

a(n) = core(n^n); \\ Michel Marcus, Nov 01 2022

Formula

a(n) = A007913(A000312(n)). - Michel Marcus, Nov 01 2022