|
%I #7 Sep 20 2023 11:53:12
%S 1,3,2,6,10,4,7,17,5,25,29,12,20,42,8,9,26,61,69,23,11,31,48,96,13,22,
%T 111,64,14,44,134,15,16,154,36,28,62,18,19,72,109,210,21,34,54,240,
%U 139,89,24,288,39,181,329,27,55,66,137,45,374,30,99,161,236,32
%N Position of A365787(n) in A024619.
%C A365787(n) = A286708(n)/A007947(A286708(n)).
%C Permutation of natural numbers.
%H Michael De Vlieger, <a href="/A365789/b365789.txt">Table of n, a(n) for n = 1..4742</a>
%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%e Let b(n) = A286708(n), rad(n) = A007947(n), and c(n) = A024619(n).
%e a(1) = 1 since b(1)/rad(b(1)) = 36/6 = 6 = c(1).
%e a(2) = 3 since b(2)/rad(b(2)) = 72/6 = 12 = c(3).
%e a(3) = 2 since b(3)/rad(b(3)) = 100/10 = 10 = c(2).
%e a(4) = 6 since b(4)/rad(b(4)) = 108/6 = 18 = c(6).
%e a(5) = 10 since b(5)/rad(b(5)) = 144/6 = 24 = c(10).
%e a(6) = 4 since b(6)/rad(b(6)) = 196/14 = 14 = c(4), etc.
%t nn = 3600;
%t s = Rest@
%t Select[Union@ Flatten@
%t Table[a^2*b^3, {b, nn^(1/3)}, {a, Sqrt[nn/b^3]}],
%t Not @* PrimePowerQ];
%t t = Select[Range[2, nn], Not @* PrimePowerQ];
%t Map[FirstPosition[t, #/(Times @@ FactorInteger[#][[All, 1]])][[1]] &, s]
%Y Cf. A007947, A024619, A120944, A286708, A365786, A365787.
%K nonn
%O 1,2
%A _Michael De Vlieger_, Sep 19 2023
|
