login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A129141
a(1)=a(2)=1. a(n) = the n-th integer from among those positive integers coprime to (a(n-1)*a(n-2)).
1
1, 1, 3, 5, 8, 13, 15, 16, 31, 19, 11, 13, 15, 29, 28, 39, 61, 29, 19, 21, 40, 97, 57, 37, 41, 26, 59, 61, 29, 31, 33, 53, 56, 81, 121, 59, 40, 97, 99, 65, 92, 119, 111, 86, 139, 93, 73, 76, 105, 227, 109, 52, 115, 151, 72, 169, 185, 81, 113, 89, 61, 63, 110, 307, 177, 100
OFFSET
1,3
EXAMPLE
a(4)*a(5) = 40. So a(6) = the sixth integer from among those positive integers coprime to 40 (i.e., the sixth integer in the sequence 1,3,7,9,11,13,17,19,21,23,...). So a(6) = 13.
MAPLE
a[1]:=1: a[2]:=1: for n from 3 to 80 do A:={}: for i from 1 to 500 do if gcd(i, a[n-1]*a[n-2])=1 then A:=A union {i} else fi od: a[n]:=A[n]: od: seq(a[n], n=1..80); # Emeric Deutsch, Apr 12 2007
CROSSREFS
Cf. A129140.
Sequence in context: A131354 A358767 A092360 * A289013 A357048 A331467
KEYWORD
nonn
AUTHOR
Leroy Quet, Mar 30 2007
EXTENSIONS
More terms from Emeric Deutsch, Apr 12 2007
STATUS
approved