|
%I #1 Sep 22 2004 03:00:00
%S 1,4,2,6,3,9,5,11,7,14,8,16,10,18,12,20,13,24,15,26,21,28,17,30,19,22,
%T 27,32,23,36,25,38,33,40,29,43,31,34,39,45,35,49,37,47,41,51,42,55,44,
%U 53,46,57,48,61,50,63,52,59,54,67,56,65,58,69,60,75,62,71,64,77,66,73
%N Permutation of natural numbers inverse to A097465.
%e a(3)=2 because A097465(2)=3.
%p A:=proc(n) option remember; local t, S; S:=({$1..1000} minus {seq(A(i),i=1..n-1)}) minus {A(n-1)-1,A(n-1)+1}; t:=min(S[]); while igcd(A(n-1),t)>1 do S:=S minus {t}; t:=min(S[]) od; t end: A(1):=1: L:=[seq(A(n),n=1..200)]: a:=n-> if member(n,L,'p') then p else 0 fi: seq(a(n),n=1..195);
%Y Cf. A097465.
%K nonn
%O 1,2
%A Alec Mihailovs (alec(AT)mihailovs.com), Aug 23 2004
|
