OFFSET
1,1
COMMENTS
PROG
(Magma)
XOR := func<a, b | Seqint([ (adigs[i] + bdigs[i]) mod 2 : i in [1..n]], 2)
where adigs := Intseq(a, 2, n)
where bdigs := Intseq(b, 2, n)
where n := 1 + Ilog2(Max([a, b, 1]))>;
for i in [1 .. 12] do
xornum:=2^(i)-2;
for j := 2^(i) +1 to 3*2^(i-1) by 2 do
xorcouple:=XOR(j, xornum);
if (IsPrime(j) and IsPrime(xorcouple)) then j;
end if;
end for;
"****";
end for;
CROSSREFS
KEYWORD
nonn
AUTHOR
Brad Clardy, Nov 04 2011
STATUS
approved