cp's OEIS Frontend

A169841 EKG sequence started at 5 instead of 2.

%I A169841 #16 Jan 21 2016 09:00:13
%S A169841 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 A169841 25,35,28,32,34,17,51,36,38,19,57,42,40,44,46,23,69,45,48,50,52,54,56,
%U A169841 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 A169841 EKG sequence started at 5 instead of 2.
%H A169841 Peter Kagey, <a href="/A169841/b169841.txt">Table of n, a(n) for n = 1..10000</a>
%p A169841 N:= 1000: # for terms before the first term > N
%p A169841 V:= Vector(N):
%p A169841 A[1]:= 5: V[5]:= 1:
%p A169841 for n from 2 do
%p A169841   S:= {seq(seq(k*p,k=1..N/p),p=numtheory:-factorset(A[n-1]))};
%p A169841   for s in sort(convert(S,list)) do
%p A169841     if V[s] = 0 then
%p A169841       A[n]:= s;
%p A169841       break
%p A169841     fi
%p A169841   od;
%p A169841   if V[s] = 1 then break fi;
%p A169841   V[s]:= 1;
%p A169841 od:
%p A169841 seq(A[i],i=1..n-1); # _Robert Israel_, Jan 18 2016
%Y A169841 For other initial terms, see A064413, A169837, A169839, A169843, A169845, A169847, A169849, A169851, A169853, A169855.
%Y A169841 Cf. A169842.
%K A169841 nonn
%O A169841 1,1
%A A169841 _T. D. Noe_ and _N. J. A. Sloane_, Jun 02 2010