A130700 a(0)=1; a(n) is the n-th integer from among those positive integers which are coprime to (a(n-1)*n).

1, 1, 3, 4, 7, 6, 17, 8, 15, 16, 23, 12, 35, 19, 33, 31, 33, 29, 55, 27, 73, 37, 49, 27, 71, 31, 57, 43, 67, 30, 109, 32, 63, 64, 71, 51, 113, 38, 79, 62, 103, 42, 145, 56, 111, 88, 105, 106, 145, 73, 127

Offset

0,3

Links

Table of n, a(n) for n=0..50.

Example

a(8)*9 = 135. Those positive integers which are coprime to 135 form the sequence: 1,2,4,7,8,11,13,14,16,17,19,22,... The 9th of these integers is 16, so a(9) = 16.

Maple

a[0]:=1: for n to 50 do CP:={}: for j while nops(CP)<=n do if gcd(n*a[n-1], j) =1 then CP:=`union`(CP, {j}) else end if end do: a[n]:=CP[n] end do: seq(a[n], n=0..50); # Emeric Deutsch, Aug 28 2007

Crossrefs

Cf. A069213.

Sequence in context: A244974 A077580 A069213 * A117134 A334127 A095001

Adjacent sequences: A130697 A130698 A130699 * A130701 A130702 A130703

Keyword

nonn

Author

Leroy Quet, Aug 12 2007

Extensions

More terms from Emeric Deutsch, Aug 28 2007

Status

approved