|
%I #16 Jan 21 2016 09:00:13
%S 5,10,2,4,6,3,9,12,8,14,7,21,15,18,16,20,22,11,33,24,26,13,39,27,30,
%T 25,35,28,32,34,17,51,36,38,19,57,42,40,44,46,23,69,45,48,50,52,54,56,
%U 49,63,60,55,65,70,58,29,87,66,62,31,93,72,64,68,74,37,111,75,78,76,80,82,41
%N EKG sequence started at 5 instead of 2.
%H Peter Kagey, <a href="/A169841/b169841.txt">Table of n, a(n) for n = 1..10000</a>
%p N:= 1000: # for terms before the first term > N
%p V:= Vector(N):
%p A[1]:= 5: V[5]:= 1:
%p for n from 2 do
%p S:= {seq(seq(k*p,k=1..N/p),p=numtheory:-factorset(A[n-1]))};
%p for s in sort(convert(S,list)) do
%p if V[s] = 0 then
%p A[n]:= s;
%p break
%p fi
%p od;
%p if V[s] = 1 then break fi;
%p V[s]:= 1;
%p od:
%p seq(A[i],i=1..n-1); # _Robert Israel_, Jan 18 2016
%Y For other initial terms, see A064413, A169837, A169839, A169843, A169845, A169847, A169849, A169851, A169853, A169855.
%Y Cf. A169842.
%K nonn
%O 1,1
%A _T. D. Noe_ and _N. J. A. Sloane_, Jun 02 2010
|
