A098552 A098550(A098550(n)).

1, 2, 3, 4, 5, 14, 7, 16, 9, 8, 26, 12, 34, 10, 15, 6, 39, 22, 28, 11, 30, 25, 35, 24, 45, 95, 51, 23, 85, 21, 20, 123, 33, 75, 55, 64, 65, 88, 17, 69, 91, 58, 27, 63, 78, 42, 81, 49, 54, 50, 13, 40, 102, 48, 72, 44, 111, 46, 41, 31, 32, 18, 56, 92, 93, 94, 77, 74, 52, 112, 62

Offset: 1

Author

Reinhard Zumkeller, Sep 14 2004

Keywords

• nonn

Comments

Integer permutation with inverse A098553; A098551(a(n)) = a(A098551(n)) = A098550(n).

Links

• Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
• Index entries for sequences that are permutations of the natural numbers

Crossrefs

Cf. A098550, A098553.

Programs

Haskell

a098552 = a098550 . a098550  -- Reinhard Zumkeller, Nov 21 2014

Mathematica

f[lst_List] := Block[{k = 4}, While[ GCD[ lst[[-2]], k] == 1 || GCD[ lst[[-1]], k] > 1 || MemberQ[lst, k], k++]; Append[lst, k]]; lst[[#]] & /@ Take[lst = Nest[f, {1, 2, 3}, 120], 71] (* Robert G. Wilson v, Nov 21 2014 *)