OFFSET
1,2
COMMENTS
a(n) = (n+1)*A006519(n)-1;
a(2*n-1) = 2*n-1, a(2*n) > 4*n.
LINKS
EXAMPLE
n=20='10100'='101'00' -> a(20)='101'00'11'='1010011'=83.
MATHEMATICA
bra1[n_]:=Module[{idn2=IntegerDigits[n, 2]}, FromDigits[Join[ idn2, Table[1, {IntegerExponent[FromDigits[idn2]]}]], 2]]; Array[bra1, 70] (* Harvey P. Dale, Sep 30 2012 *)
PROG
(Haskell)
a092525 n = f n n where
f x y = if m == 0 then f x' (2 * y + 1) else y
where (x', m) = divMod x 2
-- Reinhard Zumkeller, Oct 06 2012
(Python)
def A092525(n): return (n+1)*(~n&n-1)+n # Chai Wah Wu, Jul 07 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 07 2004
STATUS
approved