A029837 -id:A029837 - cp's OEIS frontend

A252464 a(1) = 0, a(2n) = 1 + a(n), a(2n+1) = 1 + a(A064989(2n+1)); also binary width of terms of A156552 and A243071.

0, 1, 2, 2, 3, 3, 4, 3, 3, 4, 5, 4, 6, 5, 4, 4, 7, 4, 8, 5, 5, 6, 9, 5, 4, 7, 4, 6, 10, 5, 11, 5, 6, 8, 5, 5, 12, 9, 7, 6, 13, 6, 14, 7, 5, 10, 15, 6, 5, 5, 8, 8, 16, 5, 6, 7, 9, 11, 17, 6, 18, 12, 6, 6, 7, 7, 19, 9, 10, 6, 20, 6, 21, 13, 5, 10, 6, 8, 22, 7, 5, 14, 23, 7, 8, 15, 11, 8, 24, 6, 7, 11, 12, 16, 9, 7, 25, 6, 7, 6, 26, 9, 27

Offset: 1

Antti Karttunen, Dec 20 2014

nonn

a(n) tells how many iterations of A252463 are needed before 1 is reached, i.e., the distance of n from 1 in binary trees like A005940 and A163511.

Similarly for A253553 in trees A253563 and A253565. - Antti Karttunen, Apr 14 2019

			From _Gus Wiseman_, Apr 02 2019: (Start)
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so a(n) is the size of the inner lining of the integer partition with Heinz number n, which is also the size of the largest hook of the same partition. For example, the partition with Heinz number 715 is (6,5,3), with diagram
  o o o o o o
  o o o o o
  o o o
which has inner lining
          o o
      o o o
  o o o
and largest hook
  o o o o o o
  o
  o
both of which have size 8, so a(715) = 8.
(End)

Antti Karttunen, Table of n, a(n) for n = 1..8192

Cf. A000040, A000079, A001221, A001222, A005940, A029837, A061395, A156552 (A005941), A163511, A243071, A252461, A252463, A252735, A252736, A252759, A253553, A253563, A253565, A297113, A297155, A297167, A324870, A324872.
Right edge of irregular triangle A265146.
Cf. also A246348.
Cf. A052126, A056239, A093641, A112798, A257990, A325133, A325134, A325135.

Table[If[n==1,1,PrimeOmega[n]+PrimePi[FactorInteger[n][[-1,1]]]]-1,{n,100}] (* Gus Wiseman, Apr 02 2019 *)

A061395(n) = if(n>1, primepi(vecmax(factor(n)[, 1])), 0);
A252464(n) = (bigomega(n) + A061395(n) - 1); \\ Antti Karttunen, Apr 14 2019

from sympy import primepi, primeomega, primefactors
def A252464(n): return primeomega(n)+primepi(max(primefactors(n)))-1 if n>1 else 0 # Chai Wah Wu, Jul 17 2023

a(1) = 0; for n > 1: a(n) = 1 + a(A252463(n)).
a(n) = A029837(1+A243071(n)). [a(n) = binary width of terms of A243071.]
a(n) = A029837(A005941(n)) = A029837(1+A156552(n)). [Also binary width of terms of A156552.]
Other identities. For all n >= 1:
a(A000040(n)) = n.
a(A001248(n)) = n+1.
a(A030078(n)) = n+2.
And in general, a(prime(n)^k) = n+k-1.
a(A000079(n)) = n. [I.e., a(2^n) = n.]
For all n >= 2:
a(n) = A001222(n) + A061395(n) - 1 = A001222(n) + A252735(n) = A061395(n) + A252736(n) = 1 + A252735(n) + A252736(n).
a(n) = A325134(n) - 1. - Gus Wiseman, Apr 02 2019
From Antti Karttunen, Apr 14 2019: (Start)
a(1) = 0; for n > 1: a(n) = 1 + a(A253553(n)).
a(n) = A001221(n) + A297167(n) = A297113(n) + A297155(n).
(End).

A035100 Number of bits in binary expansion of n-th prime.

Original entry on oeis.org

2, 2, 3, 3, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9

Offset: 1

N. J. A. Sloane

Jamie Morken, Table of n, a(n) for n = 1..10000 (first 1000 terms from G. C. Greubel)

Cf. A004676, A000523, A029837, A070939, A162145.

Digits := 40: [ seq(1+floor(evalf(log(ithprime(n))/log(2))), n=1..100) ];

f[n_]:=Length[IntegerDigits[Prime[n],2]]; Table[f[n],{n,5!}] (* Vladimir Joseph Stephan Orlovsky, Feb 07 2010 *)
IntegerLength[Prime[Range[100]],2] (* Harvey P. Dale, Nov 10 2022 *)

a(n)=#binary(prime(n)) \\ Charles R Greathouse IV, Sep 23 2012

from sympy import prime
def A035100(n): return prime(n).bit_length() # Chai Wah Wu, Nov 19 2024

A053738 If k is in sequence then 2k and 2k+1 are not (and 1 is in the sequence); numbers with an odd number of digits in binary.

Original entry on oeis.org

1, 4, 5, 6, 7, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109

Offset: 1

Henry Bottomley, Apr 06 2000

Runs of successive numbers have lengths which are powers of 4.
Apparently, for any m>=1, 2^m is the largest power of 2 dividing sum(k=1,n,binomial(2k,k)^m) if and only if n is in the sequence. - Benoit Cloitre, Apr 27 2003
Numbers that begin with a 1 in base 4. - Michel Marcus, Dec 05 2013
The lower and upper asymptotic densities of this sequence are 1/3 and 2/3, respectively. - Amiram Eldar, Feb 01 2021

Robert Israel, Table of n, a(n) for n = 1..10000
Manfred Madritsch and Stephan Wagner, A central limit theorem for integer partitions, Monatsh. Math., Vol. 161, No. 1 (2010), pp. 85-114; alternative link. Section 4.3.

Complement of A053754.
Cf. A029837, A079112.
Cf. A000302, A083420.

seq(seq(i,i=4^k..2*4^k-1),k=0..5); # Robert Israel, Dec 28 2016

Select[Range[110],OddQ[IntegerLength[#,2]]&] (* Harvey P. Dale, Sep 29 2012 *)

isok(n, b=4) = digits(n, b)[1] == 1; \\ Michel Marcus, Dec 05 2013

PARI
```
a(n) = n + 1<Kevin Ryde, Mar 27 2021
```

G.f.: x/(1-x)^2 + Sum_{k>=1} 2^(2k-1)*x^((4^k+2)/3)/(1-x). - Robert Israel, Dec 28 2016

A373948 Run-compression encoded as a transformation of compositions in standard order.

Original entry on oeis.org

0, 1, 2, 1, 4, 5, 6, 1, 8, 9, 2, 5, 12, 13, 6, 1, 16, 17, 18, 9, 20, 5, 22, 5, 24, 25, 6, 13, 12, 13, 6, 1, 32, 33, 34, 17, 4, 37, 38, 9, 40, 41, 2, 5, 44, 45, 22, 5, 48, 49, 50, 25, 52, 13, 54, 13, 24, 25, 6, 13, 12, 13, 6, 1, 64, 65, 66, 33, 68, 69, 70, 17, 72

Offset: 0

Gus Wiseman, Jun 24 2024

nonn

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.
We define the (run-) compression of a sequence to be the anti-run obtained by reducing each run of repeated parts to a single part. Alternatively, compression removes all parts equal to the part immediately to their left. For example, (1,1,2,2,1) has compression (1,2,1).
For the present sequence, the a(n)-th composition in standard order is obtained by compressing the n-th composition in standard order.

			The standard compositions and their compressions begin:
   0: ()        -->  0: ()
   1: (1)       -->  1: (1)
   2: (2)       -->  2: (2)
   3: (1,1)     -->  1: (1)
   4: (3)       -->  4: (3)
   5: (2,1)     -->  5: (2,1)
   6: (1,2)     -->  6: (1,2)
   7: (1,1,1)   -->  1: (1)
   8: (4)       -->  8: (4)
   9: (3,1)     -->  9: (3,1)
  10: (2,2)     -->  2: (2)
  11: (2,1,1)   -->  5: (2,1)
  12: (1,3)     --> 12: (1,3)
  13: (1,2,1)   --> 13: (1,2,1)
  14: (1,1,2)   -->  6: (1,2)
  15: (1,1,1,1) -->  1: (1)

Gus Wiseman, Statistics, classes, and transformations of standard compositions

Positions of 1's are A000225.
The image is A333489, counted by A003242.
Sum of standard composition for a(n) is given by A373953, length A124767.
A037201 gives compression of first differences of primes, halved A373947.
A066099 lists the parts of all compositions in standard order.
A114901 counts compositions with no isolated parts.
A116861 counts partitions by compressed sum, by length A116608.
A240085 counts compositions with no unique parts.
A333755 counts compositions by compressed length.
A373949 counts compositions by compressed sum, opposite A373951.
Cf. A106356, A124762, A238130, A238279, A272919, A333381, A333382, A333627, A373952, A373954.

stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
stcinv[q_]:=Total[2^(Accumulate[Reverse[q]])]/2;
Table[stcinv[First/@Split[stc[n]]],{n,0,30}]

A029837(a(n)) = A373953(n).

A000120(a(n)) = A124767(n).

A087207 A binary representation of the primes that divide a number, shown in decimal.

Original entry on oeis.org

0, 1, 2, 1, 4, 3, 8, 1, 2, 5, 16, 3, 32, 9, 6, 1, 64, 3, 128, 5, 10, 17, 256, 3, 4, 33, 2, 9, 512, 7, 1024, 1, 18, 65, 12, 3, 2048, 129, 34, 5, 4096, 11, 8192, 17, 6, 257, 16384, 3, 8, 5, 66, 33, 32768, 3, 20, 9, 130, 513, 65536, 7, 131072, 1025, 10, 1, 36, 19, 262144, 65, 258

Offset: 1

Mitch Cervinka (puritan(AT)planetkc.com), Oct 26 2003

The binary representation of a(n) shows which prime numbers divide n, but not the multiplicities. a(2)=1, a(3)=10, a(4)=1, a(5)=100, a(6)=11, a(10)=101, a(30)=111, etc.
For n > 1, a(n) gives the (one-based) index of the column where n is located in array A285321. A008479 gives the other index. - Antti Karttunen, Apr 17 2017
From Antti Karttunen, Jun 18 & 20 2017: (Start)
A268335 gives all n such that a(n) = A248663(n); the squarefree numbers (A005117) are all the n such that a(n) = A285330(n) = A048675(n).
For all n > 1 for which the value of A285331(n) is well-defined, we have A285331(a(n)) <= floor(A285331(n)/2), because then n is included in the binary tree A285332 and a(n) is one of its ancestors (in that tree), and thus must be at least one step nearer to its root than n itself.
Conjecture: Starting at any n and iterating the map n -> a(n), we will always reach 0 (see A288569). This conjecture is equivalent to the conjecture that at any n that is neither a prime nor a power of two, we will eventually hit a prime number (which then becomes a power of two in the next iteration). If this conjecture is false then sequence A285332 cannot be a permutation of natural numbers. On the other hand, if the conjecture is true, then A285332 must be a permutation of natural numbers, because all primes and powers of 2 occur in definite positions in that tree. This conjecture also implies the conjectures made in A019565 and A285320 that essentially claim that there are neither finite nor infinite cycles in A019565.
If there are any 2-cycles in this sequence, then both terms of the cycle should be present in A286611 and the larger one should be present in A286612.
(End)
Binary rank of the distinct prime indices of n, where the binary rank of an integer partition y is given by Sum_i 2^(y_i-1). For all prime indices (with multiplicity) we have A048675. - Gus Wiseman, May 25 2024

			a(38) = 129 because 38 = 2*19 = prime(1)*prime(8) and 129 = 2^0 + 2^7 (in binary 10000001).
a(140) = 13, binary 1101 because 140 is divisible by the first, third and fourth primes and 2^(1-1) + 2^(3-1) + 2^(4-1) = 13.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000 [First 1000 terms from _T. D. Noe_]
Index entries for sequences related to binary expansion of n
Index entries for sequences computed from indices in prime factorization

For partial sums see A288566.
Cf. A000040, A000120, A001221, A005117, A008479, A019565, A055396, A285320, A285321, A285329, A285330, A285332.
Sequences with related definitions: A007947, A008472, A027748, A048675, A248663, A276379 (same sequence shown in base 2), A288569, A289271, A297404.
Cf. A286608 (numbers n for which a(n) < n), A286609 (n for which a(n) > n), and also A286611, A286612.
A003986, A003961, A059896 are used to express relationship between terms of this sequence.
Related to A267116 via A225546.
Positions of particular values are: A000079\{1} (1), A000244\{1} (2), A033845 (3), A000351\{1} (4), A033846 (5), A033849 (6), A143207 (7), A000420\{1} (8), A033847 (9), A033850 (10), A033851 (12), A147576 (14), A147571 (15), A001020\{1} (16), A033848 (17).
A048675 gives binary rank of prime indices.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Binary indices (listed A048793):
- length A000120, complement A023416
- min A001511, opposite A000012
- sum A029931, product A096111
- max A029837 or A070939, opposite A070940
- complement A368494, sum A359400
- opposite complement A371571, sum A359359
- opposite A371572, sum A230877
Cf. A000720, A005940, A018819, A023506, A071814, A225620, A277319, A277905, A304818, A372689, A372890.

a087207 = sum . map ((2 ^) . (subtract 1) . a049084) . a027748_row
-- Reinhard Zumkeller, Jul 16 2013

a[n_] := Total[ 2^(PrimePi /@ FactorInteger[n][[All, 1]] - 1)]; a[1] = 0; Table[a[n], {n, 1, 69}] (* Jean-François Alcover, Dec 12 2011 *)

a(n) = {if (n==1, 0, my(f=factor(n), v = []); forprime(p=2, vecmax(f[,1]), v = concat(v, vecsearch(f[,1], p)!=0);); fromdigits(Vecrev(v), 2));} \\ Michel Marcus, Jun 05 2017

A087207(n)=vecsum(apply(p->1<M. F. Hasler, Jun 23 2017

from sympy import factorint, primepi
def a(n):
    return sum(2**primepi(i - 1) for i in factorint(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 06 2017

(definec (A087207 n) (if (= 1 n) 0 (+ (A000079 (+ -1 (A055396 n))) (A087207 (A028234 n))))) ;; This uses memoization-macro definec
(define (A087207 n) (A048675 (A007947 n))) ;; Needs code from A007947 and A048675. - Antti Karttunen, Jun 19 2017

Additive with a(p^e) = 2^(i-1) where p is the i-th prime. - Vladeta Jovovic, Oct 29 2003
a(n) gives the m such that A019565(m) = A007947(n). - Naohiro Nomoto, Oct 30 2003
A000120(a(n)) = A001221(n); a(n) = Sum(2^(A049084(p)-1): p prime-factor of n). - Reinhard Zumkeller, Nov 30 2003
G.f.: Sum_{k>=1} 2^(k-1)*x^prime(k)/(1-x^prime(k)). - Franklin T. Adams-Watters, Sep 01 2009
From Antti Karttunen, Apr 17 2017, Jun 19 2017 & Dec 06 2018: (Start)
a(n) = A048675(A007947(n)).
a(1) = 0; for n > 1, a(n) = 2^(A055396(n)-1) + a(A028234(n)).
A000035(a(n)) = 1 - A000035(n). [a(n) and n are of opposite parity.]
A248663(n) <= a(n) <= A048675(n). [XOR-, OR- and +-variants.]
a(A293214(n)) = A218403(n).
a(A293442(n)) = A267116(n).
A069010(a(n)) = A287170(n).
A007088(a(n)) = A276379(n).
A038374(a(n)) = A300820(n) for n >= 1.
(End)
From Peter Munn, Jan 08 2020: (Start)
a(A059896(n,k)) = a(n) OR a(k) = A003986(a(n), a(k)).
a(A003961(n)) = 2*a(n).
a(n^2) = a(n).
a(n) = A267116(A225546(n)).
a(A225546(n)) = A267116(n).
(End)

More terms from Don Reble, Ray Chandler and Naohiro Nomoto, Oct 28 2003

Name clarified by Antti Karttunen, Jun 18 2017

A053754 If k is in the sequence then 2k and 2k+1 are not (and 0 is in the sequence); when written in binary k has an even number of bits (0 has 0 digits).

Original entry on oeis.org

0, 2, 3, 8, 9, 10, 11, 12, 13, 14, 15, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148

Offset: 1

Henry Bottomley, Apr 06 2000

Runs of successive terms with same number of bits have length twice powers of 4 (A081294). [Clarified by Michel Marcus, Oct 21 2020]
The sequence A081294 counts compositions of even numbers - Gus Wiseman, Aug 12 2021
A031443 is a subsequence; A179888 is the intersection of this sequence and A032925. - Reinhard Zumkeller, Jul 31 2010
The lower and upper asymptotic densities of this sequence are 1/3 and 2/3, respectively. - Amiram Eldar, Feb 01 2021
From Gus Wiseman, Aug 10 2021: (Start)
Also numbers k such that the k-th composition in standard order (row k of A066099) has even sum. The terms and corresponding compositions begin:
  0: ()   2: (2)      8: (4)
          3: (1,1)    9: (3,1)
                     10: (2,2)
                     11: (2,1,1)
                     12: (1,3)
                     13: (1,2,1)
                     14: (1,1,2)
                     15: (1,1,1,1)
The following pertain to compositions in standard order: A000120, A029837, A070939, A066099, A124767.
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10001

Cf. A031443, A032925, A179888.
Positions of even terms in A029837 with offset 0.
The complement (the odd version) is A053738, counted by A000302.
The version for Heinz numbers of partitions is A300061, counted by A058696.
Cf. A000120, A001969, A008549, A070939, A080791, A114121.

a053754 n = a053754_list !! (n-1)
a053754_list = 0 : filter (even . a070939) [1..]
-- Reinhard Zumkeller, Apr 18 2015

Select[Range[0, 150], EvenQ @ IntegerLength[#, 2] &] (* Amiram Eldar, Feb 01 2021 *)

lista(nn) = {my(va = vector(nn)); for (n=2, nn, my(k=va[n-1]+1); while (#select(x->(x==k\2), va), k++); va[n] = k;); va;} \\ Michel Marcus, Oct 20 2020

PARI
```
a(n) = n-1 + (1<Kevin Ryde, Apr 30 2021
```

Offset corrected by Reinhard Zumkeller, Jul 30 2010

A373953 Sum of run-compression of the n-th integer composition in standard order.

Original entry on oeis.org

0, 1, 2, 1, 3, 3, 3, 1, 4, 4, 2, 3, 4, 4, 3, 1, 5, 5, 5, 4, 5, 3, 5, 3, 5, 5, 3, 4, 4, 4, 3, 1, 6, 6, 6, 5, 3, 6, 6, 4, 6, 6, 2, 3, 6, 6, 5, 3, 6, 6, 6, 5, 6, 4, 6, 4, 5, 5, 3, 4, 4, 4, 3, 1, 7, 7, 7, 6, 7, 7, 7, 5, 7, 4, 5, 6, 7, 7, 6, 4, 7, 7, 7, 6, 5, 3, 5

Offset: 0

Gus Wiseman, Jun 25 2024

nonn

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

We define the (run-) compression of a sequence to be the anti-run obtained by reducing each run of repeated parts to a single part. Alternatively, compression removes all parts equal to the part immediately to their left. For example, (1,1,2,2,1) has compression (1,2,1).

			The standard compositions and their compressions and compression sums begin:
   0: ()        --> ()      --> 0
   1: (1)       --> (1)     --> 1
   2: (2)       --> (2)     --> 2
   3: (1,1)     --> (1)     --> 1
   4: (3)       --> (3)     --> 3
   5: (2,1)     --> (2,1)   --> 3
   6: (1,2)     --> (1,2)   --> 3
   7: (1,1,1)   --> (1)     --> 1
   8: (4)       --> (4)     --> 4
   9: (3,1)     --> (3,1)   --> 4
  10: (2,2)     --> (2)     --> 2
  11: (2,1,1)   --> (2,1)   --> 3
  12: (1,3)     --> (1,3)   --> 4
  13: (1,2,1)   --> (1,2,1) --> 4
  14: (1,1,2)   --> (1,2)   --> 3
  15: (1,1,1,1) --> (1)     --> 1

Gus Wiseman, Statistics, classes, and transformations of standard compositions

Positions of 1's are A000225.
Counting partitions by this statistic gives A116861, by length A116608.
For length instead of sum we have A124767, counted by A238279 and A333755.
Compositions counted by this statistic are A373949, opposite A373951.
A037201 gives compression of first differences of primes, halved A373947.
A066099 lists the parts of all compositions in standard order.
A114901 counts compositions with no isolated parts.
A240085 counts compositions with no unique parts.
A333489 ranks anti-runs, counted by A003242.
Cf. A106356, A124762, A238130, A238343, A272919, A285981, A333381, A333382, A333627, A373952, A373954.

stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Table[Total[First/@Split[stc[n]]],{n,0,100}]

a(n) = A029837(A373948(n)).

A037861 (Number of 0's) - (number of 1's) in the base-2 representation of n.

Original entry on oeis.org

1, -1, 0, -2, 1, -1, -1, -3, 2, 0, 0, -2, 0, -2, -2, -4, 3, 1, 1, -1, 1, -1, -1, -3, 1, -1, -1, -3, -1, -3, -3, -5, 4, 2, 2, 0, 2, 0, 0, -2, 2, 0, 0, -2, 0, -2, -2, -4, 2, 0, 0, -2, 0, -2, -2, -4, 0, -2, -2, -4, -2, -4, -4, -6, 5, 3, 3, 1, 3, 1, 1, -1, 3

Offset: 0

Clark Kimberling

-Sum_{n>=1} a(n)/((2*n)*(2*n+1)) = the "alternating Euler constant" log(4/Pi) = 0.24156... - (see A094640 and Sondow 2005, 2010).

a(A072600(n)) < 0; a(A072601(n)) <= 0; a(A031443(n)) = 0; a(A072602(n)) >= 0; a(A072603(n)) > 0; a(A031444(n)) = 1; a(A031448(n)) = -1; abs(a(A089648(n))) <= 1. - Reinhard Zumkeller, Feb 07 2015

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000
Jonathan Sondow, Double integrals for Euler's constant and ln(4/Pi) and an analog of Hadjicostas's formula, arXiv:math/0211148 [math.CA], 2002-2004; Amer. Math. Monthly 112 (2005), 61-65.
Jonathan Sondow, New Vacca-Type Rational Series for Euler's Constant and Its "Alternating" Analog ln(4/Pi), arXiv:math/0508042 [math.NT], 2005; Additive Number Theory, Festschrift In Honor of the Sixtieth Birthday of Melvyn B. Nathanson (D. Chudnovsky and G. Chudnovsky, eds.), Springer, 2010, pp. 331-340.
Ralf Stephan, Divide-and-conquer generating functions. I. Elementary sequences, arXiv:math/0307027 [math.CO], 2003.
Ralf Stephan, Some divide-and-conquer sequences with (relatively) simple ordinary generating functions, 2004.
Ralf Stephan, Table of generating functions (ps file).
Ralf Stephan, Table of generating functions (pdf file).
Index entries for sequences related to binary expansion of n

Cf. A031443 for n when a(n)=0, A053738 for n when a(n) odd, A053754 for n when a(n) even, A030300 for a(n+1) mod 2.
Cf. A066879, A094640, A110625, A145037.
Cf. A072600, A072601, A072602, A072603, A031444, A031448, A089648.
See A268289 for a recurrence based on this sequence.

a037861 n = a023416 n - a000120 n  -- Reinhard Zumkeller, Aug 01 2013

A037861:= proc(n) local L;
     L:= convert(n,base,2);
     numboccur(0,L) - numboccur(1,L)
end proc:
map(A037861, [$0..100]); # Robert Israel, Mar 08 2016

Table[Count[ IntegerDigits[n, 2], 0] - Count[IntegerDigits[n, 2], 1], {n, 0, 75}]

a(n) = if (n==0, 1, 1 + logint(n, 2) - 2*hammingweight(n)); \\ Michel Marcus, May 15 2020 and Jun 16 2020

def A037861(n):
    return 2*format(n,'b').count('0')-len(format(n,'b')) # Chai Wah Wu, Mar 07 2016

From Henry Bottomley, Oct 27 2000: (Start)
a(n) = A023416(n) - A000120(n) = A029837(n) - 2*A000120(n) = 2*A023416(n) - A029837(n).
a(2*n) = a(n) + 1; a(2*n + 1) = a(2*n) - 2 = a(n) - 1. (End)
G.f. satisfies A(x) = (1 + x)*A(x^2) - x*(2 + x)/(1 + x). - Franklin T. Adams-Watters, Dec 26 2006
a(n) = b(n) for n > 0 with b(0) = 0 and b(n) = b(floor(n/2)) + (-1)^(n mod 2). - Reinhard Zumkeller, Dec 31 2007
G.f.: 1 + (1/(1 - x))*Sum_{k>=0} x^(2^k)*(x^(2^k) - 1)/(1 + x^(2^k)). - Ilya Gutkovskiy, Apr 07 2018

A353932 Irregular triangle read by rows where row k lists the run-sums of the k-th composition in standard order.

Original entry on oeis.org

1, 2, 2, 3, 2, 1, 1, 2, 3, 4, 3, 1, 4, 2, 2, 1, 3, 1, 2, 1, 2, 2, 4, 5, 4, 1, 3, 2, 3, 2, 2, 3, 4, 1, 2, 1, 2, 2, 3, 1, 4, 1, 3, 1, 1, 4, 1, 2, 2, 2, 3, 2, 2, 1, 3, 2, 5, 6, 5, 1, 4, 2, 4, 2, 6, 3, 2, 1, 3, 1, 2, 3, 3, 2, 4, 2, 3, 1, 6, 4, 2, 2, 1, 3

Offset: 1

Gus Wiseman, Jun 10 2022

Every sequence can be uniquely split into a sequence of non-overlapping runs. For example, the runs of (2,2,1,1,1,3,2,2) are ((2,2),(1,1,1),(3),(2,2)), with sums (4,3,3,4).

			Triangle begins:
  1
  2
  2
  3
  2 1
  1 2
  3
  4
  3 1
  4
  2 2
  1 3
  1 2 1
For example, composition 350 in standard order is (2,2,1,1,1,2), so row 350 is (4,3,2).

Row-sums are A029837.
Standard compositions are listed by A066099.
Row-lengths are A124767.
These compositions are ranked by A353847.
Row k has A353849(k) distinct parts.
The version for partitions is A354584, ranked by A353832.
A005811 counts runs in binary expansion.
A300273 ranks collapsible partitions, counted by A275870.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353851 counts compositions with all equal run-sums, ranked by A353848.
A353840-A353846 pertain to partition run-sum trajectory.
A353852 ranks compositions with all distinct run-sums, counted by A353850.
A353853-A353859 pertain to composition run-sum trajectory.
A353860 counts collapsible compositions.
Cf. A003242, A175413, A181819, A238279, A274174, A333381, A333489, A333755, A353835, A353839, A353863, A353864.

stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Table[Total/@Split[stc[n]],{n,0,30}]

A373954 Excess run-compression of standard compositions. Sum of all parts minus sum of compressed parts of the n-th integer composition in standard order.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 2, 1, 0, 0, 1, 3, 0, 0, 0, 1, 0, 2, 0, 2, 0, 0, 2, 1, 1, 1, 2, 4, 0, 0, 0, 1, 3, 0, 0, 2, 0, 0, 4, 3, 0, 0, 1, 3, 0, 0, 0, 1, 0, 2, 0, 2, 1, 1, 3, 2, 2, 2, 3, 5, 0, 0, 0, 1, 0, 0, 0, 2, 0, 3, 2, 1, 0, 0, 1, 3, 0, 0, 0, 1, 2, 4, 2

Offset: 0

Gus Wiseman, Jun 27 2024

nonn

			The excess compression of (2,1,1,3) is 1, so a(92) = 1.

Gus Wiseman, Statistics, classes, and transformations of standard compositions

For length instead of sum we have A124762, counted by A106356.
The opposite for length is A124767, counted by A238279 and A333755.
Positions of zeros are A333489, counted by A003242.
Positions of nonzeros are A348612, counted by A131044.
Compositions counted by this statistic are A373951, opposite A373949.
Compression of standard compositions is A373953.
Positions of ones are A373955.
A037201 gives compression of first differences of primes, halved A373947.
A066099 lists the parts of all compositions in standard order.
A114901 counts compositions with no isolated parts.
A116861 counts partitions by this statistic, by length A116608.
A240085 counts compositions with no unique parts.
A333627 takes the rank of a composition to the rank of its run-lengths.
Cf. A070939, A238130, A238343, A272919, A285981, A333213, A333381, A333382.

stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Table[Total[stc[n]]-Total[First/@Split[stc[n]]],{n,0,100}]

a(n) = A029837(n) - A373953(n).

cp's OEIS Frontend

A252464 a(1) = 0, a(2n) = 1 + a(n), a(2n+1) = 1 + a(A064989(2n+1)); also binary width of terms of A156552 and A243071.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A035100 Number of bits in binary expansion of n-th prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

A053738 If k is in sequence then 2*k and 2*k+1 are not (and 1 is in the sequence); numbers with an odd number of digits in binary.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

A373948 Run-compression encoded as a transformation of compositions in standard order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A087207 A binary representation of the primes that divide a number, shown in decimal.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Python

Scheme

Formula

Extensions

A053754 If k is in the sequence then 2*k and 2*k+1 are not (and 0 is in the sequence); when written in binary k has an even number of bits (0 has 0 digits).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Extensions

A373953 Sum of run-compression of the n-th integer composition in standard order.

A053738 If k is in sequence then 2k and 2k+1 are not (and 1 is in the sequence); numbers with an odd number of digits in binary.

A053754 If k is in the sequence then 2k and 2k+1 are not (and 0 is in the sequence); when written in binary k has an even number of bits (0 has 0 digits).