login
Primes in A000960.
1

%I #29 May 11 2026 07:33:57

%S 3,7,13,19,79,109,181,223,307,349,613,709,1009,1093,1321,1471,1483,

%T 1579,1693,1933,2161,2269,3019,3529,4603,4783,7069,7219,7369,7741,

%U 7963,8419,8923,9511,9679,11731,11827,13411,13711,13903,14503,15559,16273,16927,17203

%N Primes in A000960.

%H Robert Israel, <a href="/A393107/b393107.txt">Table of n, a(n) for n = 1..1500</a>

%p S:= [$1..10^5]: nS:= 10^5:

%p for k from 2 while k <= nS do

%p R:= [seq(i*k,i=1..floor(nS/k))];

%p S:= subsop(seq(r=NULL,r=R),S);

%p nS:= nS-floor(nS/k);

%p od:

%p select(isprime,S); # _Robert Israel_, Mar 01 2026

%Y Intersection of A000960 and A000040.

%K nonn

%O 1,1

%A _Aminos Abijad_, Feb 01 2026