%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