OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
MAPLE
a:= proc(n) local i, k, m, r;
k:= 2*n+1;
m:= k &^ k mod (2^(1+ilog2(k)));
r:= 0;
for i from 0 while (m>0 or k>0) do
r:= r +2^i* irem(m, 2, 'm') *irem(k, 2, 'k')
od; r
end:
seq(a(n), n=0..100); # Alois P. Heinz, Jun 21 2012
MATHEMATICA
Table[BitAnd[n, n^n], {n, 1, 141, 2}] (* Harvey P. Dale, Nov 26 2014 *)
PROG
(Python)
print([k**k & k for k in range(1, 222, 2)])
CROSSREFS
KEYWORD
AUTHOR
Alex Ratushnyak, Jun 14 2012
STATUS
approved