|
%I #17 Nov 11 2020 08:24:25
%S 3,5,7,10,11,13,14,17,19,22,23,26,28,29,31,34,37,38,41,43,44,46,47,52,
%T 53,58,59,61,62,66,67,68,71,73,74,76,78,79,82,83,86,89,92,94,97,101,
%U 102,103,106,107,109,113,114,116,118,122,124,127,130,131
%N Numbers of the form p*(q-1) for primes p, q with p>q.
%C Sorted A173706. Possible values from the formula for crossing numbers of (p,q) torus knots.
%D Peter R. Cromwell, Knots and Links, Cambridge University Press, November 15, 2004, p.255.
%H Alois P. Heinz, <a href="/A173708/b173708.txt">Table of n, a(n) for n = 1..1000</a>
%F {p*(q-1) : p, q prime and p>q}.
%p a:= proc(n) option remember;
%p local k, p, q;
%p if n=1 then 3
%p else for k from a(n-1)+1 do
%p for p in numtheory[factorset](k) do
%p q:= k/p+1;
%p if isprime (q) and p>q then return k fi
%p od
%p od
%p fi
%p end:
%p seq (a(n), n=1..60); # _Alois P. Heinz_, Nov 26 2010
%t Reap[For[k = 1, k <= 200, k++, s = Solve[p > q && k == p(q-1), {p, q}, Primes]; If[s != {}, Print[{k, p, q} /. s // First]; Sow[k]]]][[2, 1]] (* _Jean-François Alcover_, Nov 11 2020 *)
%Y Cf. A000040, A006093, A037167, A173706.
%K nonn,easy
%O 1,1
%A _Jonathan Vos Post_, Nov 25 2010
%E More terms from _Alois P. Heinz_, Nov 26 2010
|
