%I #16 Jul 02 2022 14:43:46

%S 1,5,11,25,101,55,111,125,1001,101,1011,275,1101,555,1111,625,10001,

%T 5005,10011,505,481,5055,10111,1375,11001,5505,11011,2775,11101,1111,

%U 11111,3125,9091,50005,100011,25025,100101,50055,100111,2525,101001,2405,101011,25275

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

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

%t a[n_] := (b = FromDigits[IntegerDigits[n, 2]])/GCD[n, b]; Array[a, 50] (* _Amiram Eldar_, Jun 08 2022 *)

%o (PARI) a(n) = numerator(fromdigits(binary(n))/n); \\ _Kevin Ryde_, Jun 09 2022

%Y Cf. A007088, A076775, A354837.

%K nonn,base,easy

%O 1,2

%A _Ctibor O. Zizka_, Jun 08 2022