OFFSET
1,1
MAPLE
a[1]:=5:i:=2:for k from 6 to 3000 do if isprime(k) and isprime(k+2) then b:=k:c:=k+2: if b mod 4 = 1 then a[i]:=b:i:=i+1 fi: if c mod 4 = 1 then a[i]:=c:i:=i+1 fi:fi od: seq(a[n], n=1..i-1);
MATHEMATICA
Select[Union[Flatten[Select[Partition[Prime[Range[300]], 2, 1], #[[2]]-#[[1]]==2&]]], Mod[ #, 4]==1&] (* Harvey P. Dale, Oct 16 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Miklos Kristof, Sep 27 2006
STATUS
approved