A354838 For n >= 1, a(n) = A007088(n)/gcd(n, A007088(n)).

1, 5, 11, 25, 101, 55, 111, 125, 1001, 101, 1011, 275, 1101, 555, 1111, 625, 10001, 5005, 10011, 505, 481, 5055, 10111, 1375, 11001, 5505, 11011, 2775, 11101, 1111, 11111, 3125, 9091, 50005, 100011, 25025, 100101, 50055, 100111, 2525, 101001, 2405, 101011, 25275

Offset

1,2

Links

Table of n, a(n) for n=1..44.

Example

a(6) = A007088(6)/gcd(6, A007088(6)) = 110/gcd(6,110) = 110/2 = 55.

Mathematica

a[n_] := (b = FromDigits[IntegerDigits[n, 2]])/GCD[n, b]; Array[a, 50] (* Amiram Eldar, Jun 08 2022 *)

Prog

(PARI) a(n) = numerator(fromdigits(binary(n))/n); \\ Kevin Ryde, Jun 09 2022

Crossrefs

Cf. A007088, A076775, A354837.

Sequence in context: A050344 A318033 A326161 * A038253 A372918 A274681

Adjacent sequences: A354835 A354836 A354837 * A354839 A354840 A354841

Keyword

nonn,base,easy

Author

Ctibor O. Zizka, Jun 08 2022

Status

approved