login
Prime factorization of n+1 encoded with the run lengths of binary expansion.
17

%I #7 Aug 04 2014 15:19:58

%S 0,1,2,3,5,4,10,7,6,11,21,8,42,20,9,15,85,12,170,23,22,43,341,16,13,

%T 84,14,40,682,19,1365,31,41,171,18,24,2730,340,86,47,5461,44,10922,87,

%U 17,683,21845,32,26,27,169,168,43690,28,45,80,342,1364,87381,39,174762

%N Prime factorization of n+1 encoded with the run lengths of binary expansion.

%C a(2n) = 1 or 2 mod 4 and a(2n+1) = 0 or 3 mod 4 for all n > 1

%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>

%e a(1) = 1 as 2 = 2^1, a(2) = 2 (10 in binary) as 3 = 3^1 * 2^0, a(3) = 3 (11) as 4 = 2^2, a(4) = 5 (101) as 5 = 5^1 * 3^0 * 2^0, a(5) = 4 (100) as 6 = 3^1 * 2^1, a(8) = 6 (110) as 9 = 3^2 * 2^0, a(11) = 8 (1000) as 12 = 3^1 * 2^2, a(89) = 35 (100011) as 90 = 5^1 * 3^2 * 2^1, a(90) = 90 (1011010) as 91 = 13^1 * 11^0 * 7^1 * 5^0 * 3^0 * 2^0.

%e The binary expansion of a(n) begins from the left with as many 1's as is the exponent of the largest prime present in the factorization of n+1 and from then on follows runs of ej+1 zeros and ones alternatively, where ej are the corresponding exponents of the successively lesser primes (0 if that prime does not divide n+1).

%o (Haskell)

%o import Data.List (elemIndex); import Data.Maybe (fromJust)

%o a075158 = fromJust . (`elemIndex` a075157_list)

%o -- _Reinhard Zumkeller_, Aug 04 2014

%Y Inverse of A075157. a(n) = A075160(n+1)-1. a(A006093(n)) = A000975(n). Cf. A059884.

%K nonn

%O 0,3

%A _Antti Karttunen_, Sep 13 2002