A147661 - OEIS

%I #13 Nov 01 2022 07:13:16

%S 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,

%T 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,

%U 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

%N a(n) = squarefree part of n^n.

%H Alois P. Heinz, <a href="/A147661/b147661.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A007913(A000312(n)). - _Michel Marcus_, Nov 01 2022

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

%o (PARI) a(n) = core(n^n); \\ _Michel Marcus_, Nov 01 2022

%Y Cf. A000312, A007913, A055204, A147662.

%K nonn

%O 1,3

%A _Artur Jasinski_, Nov 09 2008