|
|
A073932
|
|
Define f(n) = n - largest nontrivial divisor of n or f(n) = n-1 if n is a prime [that is, f(n) = A060681(n)]. Form a triangle in which the n-th row contains terms n, f(n), f(f(n)), ... until a 1 is reached; sequence gives triangle read by rows.
|
|
5
|
|
|
1, 2, 1, 3, 2, 1, 4, 2, 1, 5, 4, 2, 1, 6, 3, 2, 1, 7, 6, 3, 2, 1, 8, 4, 2, 1, 9, 6, 3, 2, 1, 10, 5, 4, 2, 1, 11, 10, 5, 4, 2, 1, 12, 6, 3, 2, 1, 13, 12, 6, 3, 2, 1, 14, 7, 6, 3, 2, 1, 15, 10, 5, 4, 2, 1, 16, 8, 4, 2, 1, 17, 16, 8, 4, 2, 1, 18, 9, 6, 3, 2, 1, 19, 18, 9, 6, 3, 2, 1, 20, 10, 5, 4, 2, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
EXAMPLE
|
Triangle begins:
1;
2, 1;
3, 2, 1;
4, 2, 1;
5, 4, 2, 1;
6, 3, 2, 1;
7, 6, 3, 2, 1;
8, 4, 2, 1;
9, 6, 3, 2, 1;
10, 5, 4, 2, 1;
|
|
MAPLE
|
j := 1:a[1] := 1:for i from 2 to 50 do n := i:j := j+1:a[j] := n:while(n>1) do if isprime(n) then r := n-1: else r := n-n/ifactors(n)[2][1][1]; fi; n := r:j := j+1:a[j] := n: od:od:seq(a[k], k=1..j);
|
|
MATHEMATICA
|
Array[If[# == 1, {1}, NestWhileList[If[PrimeQ@ #, # - 1, # - #/FactorInteger[#][[1, 1]] ] &, #, # > 1 &]] &, 20] // Flatten (* Michael De Vlieger, Apr 15 2020 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,tabf
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|