OFFSET
1,1
COMMENTS
Previous name was: Primes in A160668 in order of appearance.
PROG
(UBASIC)
10 'recipseq, Enoch Haga, May 22 2009
20 N=3:print N:C=2
30 A=3:S=sqrt(N)
40 B=N/A
50 if A*B=int(N) then 70
60 A=A+2:if A<S then 40
70 if N=prmdiv(N) then print N; :else 130
80 if alen(N)=1 then print 10^1-N; :P=prmdiv(10^1-N):goto 120
90 if alen(N)=2 then print 10^2-N; :P=prmdiv(10^2-N):goto 120
100 if alen(N)=3 then print 10^3-N; :P=prmdiv(10^3-N):goto 120
110 if alen(N)=4 then print 10^4-N; :P=prmdiv(10^4-N)
120 print P; C:C=C+1:stop
130 N=N+2:S=sqrt(N):goto 40
CROSSREFS
KEYWORD
dead
AUTHOR
Enoch Haga, May 23 2009
STATUS
approved