A339107 a(1) = 3; for n>1, a(n) = the smallest positive number not occurring earlier such that a(n-1)*a(n) is divisible by a(n-1)+a(n), where a(n) is not a prime or 4.
3, 6, 12, 24, 8, 56, 42, 21, 28, 70, 30, 15, 10, 40, 60, 20, 80, 48, 16, 112, 84, 14, 35, 140, 105, 120, 72, 9, 18, 36, 45, 90, 135, 108, 54, 27, 216, 270, 180, 144, 240, 160, 96, 32, 224, 168, 126, 63, 378, 189, 252, 315, 210, 231, 132, 44, 77, 462, 22, 99, 198, 165, 110, 374, 204, 68, 221, 2652
Offset: 1
Keywords
Examples
a(2) = 6. The only candidate for a(2) that satisfies 3*a(2) being divisible by 3+a(2) is a(2) = 6. a(3) = 12. The candidates for a(3) given a(2) = 6 are 3,6,12,30, all of which satisfy 6*a(3) being divisible by 6+a(3). 3 and 6 have already appeared so the next smallest candidate is chosen, being 12. a(4) = 24. The candidates for a(4) given a(3) = 12 are 4,6,12,24,36,60,132, all of which satisfy 12*a(4) being divisible by 12+a(4). 4 is not allowed as there would be no candidates for a(5), and 6 and 12 have already appeared, so the next smallest candidate is chosen, being 24. a(5) = 8. There are 10 candidates for a(5) given a(4) = 24, the smallest that has not appeared is 8.
Links
- Robert Israel, Table of n, a(n) for n = 1..10000
- Scott R. Shannon, Line graph of the first 100000 terms.
Programs
-
Maple
R:= 3: a:= 3: S:= {4,3}: for i from 2 to 100 do Cands:= remove (t -> t < 1 or isprime(t), map(`-`,numtheory:-divisors(a^2),a) minus S); a:= min(Cands); R:= R, a; S:= S union {a}; od: R; # Robert Israel, Mar 23 2023
Comments