cp's OEIS Frontend

A365901 Irregular triangle read by rows giving trajectory beginning with A002182(n) under recursion of x -> f(x) until reaching 1, where f(x) = x/rad(x), rad(x) = A007947(x).

%I A365901 #9 Oct 12 2023 17:50:17
%S A365901 1,2,1,4,2,1,6,1,12,2,1,24,4,2,1,36,6,1,48,8,4,2,1,60,2,1,120,4,2,1,
%T A365901 180,6,1,240,8,4,2,1,360,12,2,1,720,24,4,2,1,840,4,2,1,1260,6,1,1680,
%U A365901 8,4,2,1,2520,12,2,1,5040,24,4,2,1,7560,36,6,1
%N A365901 Irregular triangle read by rows giving trajectory beginning with A002182(n) under recursion of x -> f(x) until reaching 1, where f(x) = x/rad(x), rad(x) = A007947(x).
%C A365901 Let h(n) = A002182(n).
%C A365901 Since highly composite numbers h(n) are products of primorials (i.e., in A025487), the squarefree kernel is always a primorial (i.e., in A002110), and the trajectory always reaches 1.
%H A365901 Michael De Vlieger, <a href="/A365901/b365901.txt">Table of n, a(n) for n = 1..11397</a> (rows n = 1..1200, flattened)
%H A365901 Michael De Vlieger, <a href="/A365901/a365901.png">Log log scatterplot of log_10(a(n))</a>, n = 1..1053678, showing 2^16 rows.
%F A365901 Length of row n is A112779(n)+1.
%F A365901 T(n,2) = A301413(n), n > 1.
%e A365901 Row 1 = {1} since h(1) = 1, already 1.
%e A365901 Row 2 = {2, 1} since h(2) = 2, 2/rad(2) = 2/2 = 1, reaching 1.
%e A365901 Row 3 = {4, 2, 1} since h(3) = 4, 4/rad(4) = 4/2 = 2, and we follow the trajectory of 2 thereafter.
%e A365901 Row 6 = {24, 4, 2, 1} since h(6) = 24, 24/rad(24) = 24/6 = 4, and we follow the trajectory of 4 thereafter.
%e A365901 Row 14 = {720, 24, 4, 2, 1} since h(14) = 720, 720/rad(720) = 720/30 = 24, which appends row 6 thereafter.
%e A365901 .
%e A365901 First rows of this sequence:
%e A365901 row 1:   1
%e A365901     2:   2, 1
%e A365901     3:   4, 2, 1
%e A365901     4:   6, 1
%e A365901     5:  12, 2, 1
%e A365901     6:  24, 4, 2, 1
%e A365901     7:  36, 6, 1
%e A365901     8:  48, 8, 4, 2, 1
%e A365901     9:  60, 2, 1
%e A365901    10: 120, 4, 2, 1
%e A365901    11: 180, 6, 1
%e A365901    12: 240, 8, 4, 2, 1
%e A365901    ...
%t A365901 a2182 = Import["https://oeis.org/A002182/b002182.txt", "Data"][[All, -1]]; Array[NestWhileList[#/(Times @@ FactorInteger[#][[All, 1]]) &, a2182[[#]], # > 1 &] &, 20] // TableForm
%Y A365901 Cf. A002110, A002182, A007947, A112779, A301413, A365900.
%K A365901 nonn,tabf
%O A365901 1,2
%A A365901 _Michael De Vlieger_, Oct 11 2023