login
EKG sequence started at 5 instead of 2.
13

%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