OFFSET
1,2
COMMENTS
The sequence contains also Sierpiński numbers (i.e., 4521731193704761, 60428287050225649).
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local k;
for k from 0 to ilog2(n) do
if isprime(n - 2^k) then return false fi
od:
true
end proc:
select(filter, [seq((2*i+1)^2, i=0..10^4)]); # Robert Israel, Sep 07 2016
MATHEMATICA
filterQ[n_] := Module[{k}, For[k = 0, k <= Log[2, n], k++, If[PrimeQ[n - 2^k], Return[False]]]; True];
Select[Table[(2i+1)^2, {i, 0, 10^4}], filterQ] (* Jean-François Alcover, Oct 06 2020, after Maple *)
PROG
(Magma) lst:=[]; for s in [1..2015 by 2] do n:=s^2; x:=0; repeat x+:=1; a:=n-2^x; until a lt 1 or IsPrime(a); if a lt 1 then Append(~lst, n); end if; end for; lst;
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Sep 06 2016
STATUS
approved