OFFSET
1,14
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
FORMULA
a(n) is half the Hamming distance between the binary expansion of prime(n) and its reversal.
EXAMPLE
a(5) = 1 since prime(5) = 11 = 1011_2 becomes a palindrome if we change the third bit to 0.
MATHEMATICA
A144477[n_]:=With[{p=IntegerDigits[Prime[n], 2]}, HammingDistance[p, Reverse[p]]/2]; Array[A144477, 100] (* Paolo Xausa, Nov 13 2023 *)
PROG
(PARI)
HD(p)=
{
v=binary(p); H=0; j=#v;
for(k=1, #v, H+=abs(v[k]-v[j]); j--);
return(H)
};
for(n=1, 100, p=prime(n); an=HD(p)/2; print1(an, ", "))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Washington Bomfim, Jan 15 2011, following a suggestion from Joerg Arndt.
EXTENSIONS
Edited by N. J. A. Sloane, Apr 23 2020 at the suggestion of Harvey P. Dale
STATUS
approved