A008469 -id:A008469 - cp's OEIS frontend

A008467 a(n) = n OR n^2 (applied to ternary expansions).

1, 5, 12, 16, 26, 42, 52, 71, 90, 100, 122, 147, 178, 206, 231, 259, 296, 342, 370, 401, 453, 484, 539, 591, 646, 701, 756, 784, 842, 930, 961, 1025, 1095, 1186, 1259, 1332, 1369, 1445, 1524, 1609, 1691, 1770, 1852, 1943, 2070, 2152, 2237, 2316, 2401, 2510

Offset: 1

Olivier Gérard

OR(a, b) = max(a, b) for ternary digits a and b. - Sean A. Irvine, Mar 29 2018

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A007745, A008468, A008469.

a:= proc(n) option remember; local g, h, i, r;
      g, h, r:= n, n^2, 0; for i from 0 while h>0 do
      r:= r + 3^i * max(irem(g, 3, 'g'), irem(h, 3, 'h')) od; r
    end:
seq(a(n), n=1..140);  # Alois P. Heinz, Mar 29 2018

(* program should not be used to extend this sequence *)
b=3; d=9; Mb=Array[ b^(d-#)&, d ]; Map[ Plus@@(Mb*#)&, Array[ Max /@ Transpose[ {IntegerDigits[ #, b, d ], IntegerDigits[ #^2, b, d ]} ]&, 140 ] ]

A008468 a(n) = n OR n^3 (applied to binary expansions).

Original entry on oeis.org

1, 10, 27, 68, 125, 222, 343, 520, 729, 1002, 1339, 1740, 2205, 2750, 3375, 4112, 4913, 5850, 6875, 8020, 9277, 10654, 12183, 13848, 15641, 17594, 19707, 21980, 24413, 27006, 29791, 32800, 35937, 39338, 42875, 46692, 50685, 54910, 59319, 64040

Offset: 1

Olivier Gérard

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A007745, A008467, A008469.

[BitwiseOr(n,n^3): n in [1..40]]; // Bruno Berselli, Mar 30 2018

a:= n-> Bits[Or](n, n^3):
seq(a(n), n=1..100);  # Alois P. Heinz, Mar 29 2018

(* program should not be used to extend this sequence *)
d=16; Mb=Array[ 2^(d-#)&, d ]; Map[ Plus@@(Mb*#)&, Array[ IntegerDigits[ #, 2, d ]+IntegerDigits[ #^3, 2, d ]&, 40 ]/. (2->1) ]
f[n_]:=BitOr[n, n^3]; Array[f, 10000, 1] (* Vincenzo Librandi, Mar 30 2018 *)

a(n) = bitor(n, n^3); \\ Michel Marcus, Mar 30 2018

cp's OEIS Frontend

A008467 a(n) = n OR n^2 (applied to ternary expansions).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica