login
A073933
Number of terms in n-th row of triangle in A073932.
8
1, 2, 3, 3, 4, 4, 5, 4, 5, 5, 6, 5, 6, 6, 6, 5, 6, 6, 7, 6, 7, 7, 8, 6, 7, 7, 7, 7, 8, 7, 8, 6, 8, 7, 8, 7, 8, 8, 8, 7, 8, 8, 9, 8, 8, 9, 10, 7, 9, 8, 8, 8, 9, 8, 9, 8, 9, 9, 10, 8, 9, 9, 9, 7, 9, 9, 10, 8, 10, 9, 10, 8, 9, 9, 9, 9, 10, 9, 10, 8, 9, 9, 10, 9, 9, 10, 10, 9, 10, 9, 10, 10, 10, 11, 10, 8
OFFSET
1,2
LINKS
FORMULA
From Antti Karttunen, Aug 23 2017: (Start)
a(1) = 1; for n > 1, a(n) = 1 + a(A060681(n)).
a(n) = 1 + A064097(n).
(End)
MAPLE
a[1] := 1:for i from 2 to 500 do n := i:s := 1:while(n>1) do if isprime(n) then r := n-1: else r := n-n/ifactors(n)[2][1][1]; fi; n := r:s := s+1:od:a[i] := s:od:seq(a[k], k=1..500);
PROG
(Scheme)
(define (A073933 n) (if (= 1 n) n (+ 1 (A073933 (A060681 n)))))
(define (A060681 n) (- n (A032742 n))) ;; See also code under A032742
;; Antti Karttunen, Aug 23 2017
CROSSREFS
One more than A064097.
Sequence in context: A060607 A265690 A061339 * A056792 A292127 A227861
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Aug 19 2002
EXTENSIONS
More terms from Sascha Kurz, Aug 23 2002
Offset corrected from 0 to 1 by Antti Karttunen, Aug 23 2017
STATUS
approved