OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
MAPLE
a:= n-> Bits[Or](n, n^2):
seq(a(n), n=1..100); # Alois P. Heinz, Mar 29 2018
MATHEMATICA
(* program should not be used to extend this sequence *)
d = 15; Mb = Array[ 2^(d - #)&, d ]; Map[ Plus@@(Mb*#)&, Array[ IntegerDigits[ #, 2, d ] + IntegerDigits[ #^2, 2, d ]&, 180 ]/. (2->1) ] (* Mocquard and Bouchon *)
Table[BitOr[n, n^2], {n, 100}] (* Vladimir Joseph Stephan Orlovsky, Jul 19 2011 *)
PROG
(Haskell)
import Data.Bits ((.|.))
a007745 n = n .|. (n ^ 2) -- Reinhard Zumkeller, Apr 28 2014
(PARI) a(n) = bitor(n, n^2); \\ Michel Marcus, Mar 30 2018
(Magma) [BitwiseOr(n, n^2): n in [1..50]]; // Bruno Berselli, Mar 30 2018
CROSSREFS
KEYWORD
base,nonn,easy
AUTHOR
J. P. Mocquard and D. Bouchon (BOUCHON(AT)zeus.univ-poitiers.fr)
STATUS
approved