A098552 - OEIS

%I #14 Nov 22 2014 03:19:35

%S 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,

%T 51,23,85,21,20,123,33,75,55,64,65,88,17,69,91,58,27,63,78,42,81,49,

%U 54,50,13,40,102,48,72,44,111,46,41,31,32,18,56,92,93,94,77,74,52,112,62

%N A098550(A098550(n)).

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

%H Reinhard Zumkeller, <a href="/A098552/b098552.txt">Table of n, a(n) for n = 1..10000</a>

%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>

%t 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 *)

%o (Haskell)

%o a098552 = a098550 . a098550  -- _Reinhard Zumkeller_, Nov 21 2014

%Y  Cf. A098550, A098553.

%K nonn

%O 1,2

%A _Reinhard Zumkeller_, Sep 14 2004