OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..505
EXAMPLE
a(1)=2 since 2^2 is in A003714.
MAPLE
fibbin:= proc(n) Bits:-Xor(n, 2*n) = 3*n end proc:
select(t -> isprime(t) and fibbin(t^2), [2, seq(i, i=3..10^5, 2)]); # Robert Israel, Feb 15 2016
PROG
(PARI) msb(n)=my(k=1); while(k<=n, k<<=1); k>>1;
lista(nn) = {forprime(p=1, nn, n=p^2; k=bitand(n, n<<1); if(k, n=bitor(n, msb(k)-1), print1(p, ", ")); ); } \\ Michel Marcus, Feb 15 2016
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Vladimir Shevelev, Sep 20 2008
EXTENSIONS
More terms from Michel Marcus, Feb 15 2016
STATUS
approved