OFFSET
1,1
COMMENTS
The 'Px+1 map' is defined as follows: if there exists p = smallest prime < P which divides x then x = x/p, otherwise x = P*x + 1.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..11325 (antidiagonals 1..150, flattened).
EXAMPLE
Array begins:
[ 1] 2, 5, 11, 23, 47, 95, 191, 383, 767, ... = A153893
[ 2] 3, 10, 5, 16, 8, 4, 2, 1, 4, ... = A033478
[ 3] 5, 26, 13, 66, 33, 11, 56, 28, 14, ... = A057688
[ 4] 7, 50, 25, 5, 1, 8, 4, 2, 1, ... = A368113
[ 5] 11, 122, 61, 672, 336, 168, 84, 42, 21, ... = A368114
[ 6] 13, 170, 85, 17, 222, 111, 37, 482, 241, ... = A057684
[ 7] 17, 290, 145, 29, 494, 247, 19, 324, 162, ... = A368115
[ 8] 19, 362, 181, 3440, 1720, 860, 430, 215, 43, ... = A057685
[ 9] 23, 530, 265, 53, 1220, 610, 305, 61, 1404, ... = A057686
[10] 29, 842, 421, 12210, 6105, 2035, 407, 37, 1074, ... = A057687
... | | |
MATHEMATICA
Px1[p_, n_]:=Catch[For[i=1, i<PrimePi[p], i++, If[Divisible[n, Prime[i]], Throw[n/Prime[i]]]]; p*n+1];
A368085list[dmax_]:=With[{a=Reverse[Table[NestList[Px1[Prime[n], #]&, Prime[n], dmax-n], {n, dmax}]]}, Array[Diagonal[a, #]&, dmax, 1-dmax]];
A368085list[15] (* Generates 15 antidiagonals *)
CROSSREFS
KEYWORD
AUTHOR
Paolo Xausa, Dec 11 2023
STATUS
approved