login
A113469
a(0)=1. a(n) = number of positive integers <= n that are not coprime to a(n-1).
1
1, 0, 2, 1, 0, 5, 1, 0, 8, 4, 5, 2, 6, 8, 7, 2, 8, 8, 9, 6, 13, 1, 0, 23, 1, 0, 26, 14, 16, 14, 17, 1, 0, 33, 13, 2, 18, 24, 25, 7, 5, 8, 21, 18, 29, 1, 0, 47, 1, 0, 50, 30, 38, 27, 18, 36, 37, 1, 0, 59, 1, 0, 62, 32, 32, 32, 33, 26, 37, 1, 0, 71, 1, 0, 74, 38, 40, 46, 41, 1, 0
OFFSET
0,3
EXAMPLE
a(8) = 8. So a(9) = the number of positive integers <= 9 that are not coprime to 8 (i.e., that are divisible by 2). So a(9) = 4.
MAPLE
A113469 := proc(n) option remember ; local a, i, aprev ; if n= 0 then RETURN(1) ; fi ; a :=0 ; aprev := A113469(n-1) ; for i from 1 to n do if gcd(i, aprev) <> 1 or aprev=0 then a := a+1 ; fi ; od ; RETURN(a) ; end: for n from 0 to 80 do printf("%d, ", A113469(n)) ; od ; # R. J. Mathar, Jun 07 2007
CROSSREFS
Sequence in context: A322119 A363731 A112334 * A060137 A350262 A243821
KEYWORD
nonn
AUTHOR
Leroy Quet, May 10 2007
EXTENSIONS
More terms from R. J. Mathar, Jun 07 2007
STATUS
approved