[go: up one dir, main page]

login
A098552
3
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,2
COMMENTS
Integer permutation with inverse A098553; A098551(a(n)) = a(A098551(n)) = A098550(n).
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 *)
PROG
(Haskell)
a098552 = a098550 . a098550 -- Reinhard Zumkeller, Nov 21 2014
CROSSREFS
Sequence in context: A057158 A138986 A306294 * A245447 A361947 A037340
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Sep 14 2004
STATUS
approved