%I #11 Jan 07 2024 09:44:53
%S 1,2,6,3,15,10,70,35,35,14,154,231,3003,858,1430,715,12155,24310,
%T 461890,46189,969969,176358,4056234,676039,676039,104006,104006,7429,
%U 215441,6463230,200360130,100180065,367326905,43214930,60500902,30250451,1119266687,117817546,4594884294,11487210735
%N a(1) = 1; for n > 1, a(n) = A055231(a(n-1) * n), where A055231(k) is the powerfree part of k.
%C Terms are squarefree. - _Michael De Vlieger_, Jan 07 2024.
%e a(4) = 3 as a(3)*4 = 6*4 = 24 = 2*2*2*3, and the powerfree part of 24 is 3.
%t f[x_] := Times @@ Map[#1^(#2 Boole[#2 == 1]) & @@ # &, FactorInteger[x]]; a[n_] := f[n* a[n - 1]]; a[1] = 1; Array[a, 120] (* _Michael De Vlieger_, Jan 07 2024 *)
%Y Cf. A368825 (addition), A055231, A005117, A124010.
%K nonn
%O 1,2
%A _Scott R. Shannon_, Jan 07 2024