|
%I #10 Jun 09 2017 03:37:30
%S 2,7,17,37,52,77,107,122,157,217,247,262,292,317,367,437,497,527,542,
%T 587,632,637,737,752,857,877,887,997,1057,1087,1102,1177,1267,1277,
%U 1312,1372,1477,1507,1522,1717,1732,1757,1777,1822,1837,1927,1997,2047,2077,2117
%N If a, b are in the sequence, so is ab+3.
%H Robert Israel, <a href="/A009302/b009302.txt">Table of n, a(n) for n = 1..10000</a>
%p N:= 10000: # to get all terms <= N
%p S:= {2}:
%p Agenda:= {2};
%p while Agenda <> {} do
%p x:= Agenda[1];
%p T:= select(`<=`, map(s -> s*x+3, S) minus S, N);
%p S:= S union T;
%p Agenda:= Agenda[2..-1] union T;
%p od:
%p sort(convert(S,list)); # _Robert Israel_, Jun 08 2017
%K nonn
%O 1,1
%A _David W. Wilson_
|
