A336216 Irregular triangle of cycles of purely periodic unitary sigma aliquot sequences with their smallest member as starting number, read by rows.
6, 30, 42, 54, 60, 90, 114, 126, 1140, 1260, 1482, 1878, 1890, 2142, 2178, 2418, 2958, 3522, 3534, 4146, 4158, 3906, 3774, 4434, 4446, 3954, 3966, 3978, 3582, 18018, 22302, 24180, 29580, 35220, 35340, 41460, 41580, 39060, 37740, 44340, 44460, 39540, 39660, 39780, 35820, 32130, 40446
Offset: 1
Examples
The first cycle of size 14 starting at position 16 is: 2418, 2958, 3522, 3534, 4146, 4158, 3906, 3774, 4434, 4446, 3954, 3966, 3978, 3582. Its 7th element is the first number in this sequence smaller than its predecessor. Irregular triangle of cycles: 6 30 42 54 60 90 114 126 1140 1260 1482 1878 1890 2142 2178 2418 2958 3522 3534 4146 4158 3906 3774 4434 4446 3954 3966 3978 3582 18018 22302 ...
Programs
-
Mathematica
a063919[1] = 1; a063919[n_] := Total[Select[Divisors[n], GCD[#, n/#]==1&]] - n/;n>1 (* Jean-François Alcover *) aliquotSequence[n_] := NestWhileList[a063919, n, UnsameQ, All] a336216[n_] := Module[{list={}, listS={}, i, seq, seqS}, For[i=2, i<=n, i++, seq=aliquotSequence[i]; If[First[seq]==Last[seq], seqS=Sort[Most[seq]]; If[!MemberQ[listS, seqS], AppendTo[listS, seqS]; AppendTo[list, Most[seq]]]]]; list] (* list of cycles *) Flatten[a336216[35000]] (* data - first 11 rows of triangle *)
Comments