A326669 -id:A326669 - cp's OEIS frontend

A048793 List giving all subsets of natural numbers arranged in standard statistical (or Yates) order.

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

Offset: 0

N. J. A. Sloane

For n>0: first occurrence of n in row 2^(n-1), and when the table is seen as a flattened list at position n*2^(n-1)+1, cf. A005183. - Reinhard Zumkeller, Nov 16 2013

Row n lists the positions of 1's in the reversed binary expansion of n. Compare to triangles A112798 and A213925. - Gus Wiseman, Jul 22 2019

			From _Gus Wiseman_, Jul 22 2019: (Start)
Triangle begins:
  {}
  1
  2
  1  2
  3
  1  3
  2  3
  1  2  3
  4
  1  4
  2  4
  1  2  4
  3  4
  1  3  4
  2  3  4
  1  2  3  4
  5
  1  5
  2  5
  1  2  5
  3  5
(End)

S. Hedayat, N. J. A. Sloane and J. Stufken, Orthogonal Arrays, Springer-Verlag, NY, 1999, p. 249.

Reinhard Zumkeller, Rows n = 0..1000 of triangle, flattened

Cf. A048794.
Row lengths are A000120.
First column is A001511.
Heinz numbers of rows are A019565.
Row sums are A029931.
Reversing rows gives A272020.
Subtracting 1 from each term gives A133457; subtracting 1 and reversing rows gives A272011.
Indices of relatively prime rows are A291166 (see also A326674); arithmetic progressions are A295235; rows with integer average are A326669 (see also A326699/A326700); pairwise coprime rows are A326675.
Cf. A035327, A070939.

#include 
#include 
#define USAGE "Usage: 'A048793 num' where num is the largest number to use creating sets.\n"
#define MAX_NUM 10
#define MAX_ROW 1024
int main(int argc, char *argv[]) {
  unsigned short a[MAX_ROW][MAX_NUM]; signed short old_row, new_row, i, j, end;
  if (argc < 2) { fprintf(stderr, USAGE); return EXIT_FAILURE; }
  end = atoi(argv[1]); end = (end > MAX_NUM) ? MAX_NUM: end;
  for (i = 0; i < MAX_ROW; i++) for ( j = 0; j < MAX_NUM; j++) a[i][j] = 0;
  a[1][0] = 1; new_row = 2;
  for (i = 2; i <= end; i++) {
    a[new_row++ ][0] = i;
    for (old_row = 1; a[old_row][0] != i; old_row++) {
      for (j = 0; a[old_row][j] != 0; j++) { a[new_row][j] = a[old_row][j]; }
      a[new_row++ ][j] = i;
    }
  }
  fprintf(stdout, "Values: 0");
  for (i = 1; a[i][0] != 0; i++) for (j = 0; a[i][j] != 0; j++) fprintf(stdout, ",%d", a[i][j]);
  fprintf(stdout, "\n"); return EXIT_SUCCESS
}

a048793 n k = a048793_tabf !! n !! k
a048793_row n = a048793_tabf !! n
a048793_tabf = [0] : [1] : f [[1]] where
   f xss = yss ++ f (xss ++ yss) where
     yss = [y] : map (++ [y]) xss
     y = last (last xss) + 1
-- Reinhard Zumkeller, Nov 16 2013

T:= proc(n) local i, l, m; l:= NULL; m:= n;
      if n=0 then return 0 fi; for i while m>0 do
      if irem(m, 2, 'm')=1 then l:=l, i fi od; l
    end:
seq(T(n), n=0..50);  # Alois P. Heinz, Sep 06 2014

s[0] = {{}}; s[n_] := s[n] = Join[s[n - 1], Append[#, n]& /@ s[n - 1]]; Join[{0}, Flatten[s[6]]] (* Jean-François Alcover, May 24 2012 *)
Table[Join@@Position[Reverse[IntegerDigits[n,2]],1],{n,30}] (* Gus Wiseman, Jul 22 2019 *)

Constructed recursively: subsets that include n are obtained by appending n to all earlier subsets.

More terms from Larry Reeves (larryr(AT)acm.org), Apr 11 2000

A029931 If 2n = Sum 2^e_i, a(n) = Sum e_i.

Original entry on oeis.org

0, 1, 2, 3, 3, 4, 5, 6, 4, 5, 6, 7, 7, 8, 9, 10, 5, 6, 7, 8, 8, 9, 10, 11, 9, 10, 11, 12, 12, 13, 14, 15, 6, 7, 8, 9, 9, 10, 11, 12, 10, 11, 12, 13, 13, 14, 15, 16, 11, 12, 13, 14, 14, 15, 16, 17, 15, 16, 17, 18, 18, 19, 20, 21, 7, 8, 9, 10, 10, 11, 12, 13, 11, 12, 13, 14, 14, 15, 16

Offset: 0

N. J. A. Sloane

Write n in base 2, n = sum b(i)*2^(i-1), then a(n) = sum b(i)*i. - Benoit Cloitre, Jun 09 2002
May be regarded as a triangular array read by rows, giving weighted sum of compositions in standard order. The standard order of compositions is given by A066099. - Franklin T. Adams-Watters, Nov 06 2006
Sum of all positive integer roots m_i of polynomial {m,k} - see link [Shevelev]; see also A264613. - Vladimir Shevelev, Dec 13 2015
Also the sum of binary indices of n, where a binary index of n (A048793) is any position of a 1 in its reversed binary expansion. For example, the binary indices of 11 are {1,2,4}, so a(11) = 7. - Gus Wiseman, May 22 2024

			14 = 8+4+2 so a(7) = 3+2+1 = 6.
Composition number 11 is 2,1,1; 1*2+2*1+3*1 = 7, so a(11) = 7.
The triangle starts:
  0
  1
  2 3
  3 4 5 6
The reversed binary expansion of 18 is (0,1,0,0,1) with 1's at positions {2,5}, so a(18) = 2 + 5 = 7. - _Gus Wiseman_, Jul 22 2019

T. D. Noe, Table of n, a(n) for n = 0..1023
J.-P. Allouche and J. Shallit, The ring of k-regular sequences, Theoretical Computer Sci., 98 (1992), 163-197, ex. 10. See also DOI.
Vladimir Shevelev, The number of permutations with prescribed up-down structure as a function of two variables, INTEGERS, 12 (2012), #A1. (See Section 3, Theorem 21 and Section 8, Theorem 50)

Other sequences that are built by replacing 2^k in the binary representation with other numbers: A022290 (Fibonacci), A059590 (factorials), A073642, A089625 (primes), A116549, A326031.
Cf. A001793 (row sums), A011782 (row lengths), A059867, A066099, A124757.
Row sums of A048793 and A272020.
Cf. A035327, A056239, A291166, A295235, A326669, A326674, A326675, A326699/A326700.
Contains exactly A000009(n) copies of n.
For length instead of sum we have A000120, complement A023416.
For minimum instead of sum we have A001511, opposite A000012.
For maximum instead of sum we have A029837 or A070939, opposite A070940.
For product instead of sum we have A096111.
The reverse version is A230877, row sums of A371572.
The reverse complement is A359359, row sums of A371571.
The complement is A359400, row sums of A368494.
Numbers k such that a(k) is prime are A372689.
A014499 lists binary indices of prime numbers.
A019565 gives Heinz number of binary indices, inverse A048675.
A372471 lists binary indices of primes, row-sums A372429.
Cf. A000009, A030101, A359401, A359402.

a029931 = sum . zipWith (*) [1..] . a030308_row
-- Reinhard Zumkeller, Feb 28 2014

HammingWeight := n -> add(i, i = convert(n, base, 2)):
a := proc(n) option remember; `if`(n = 0, 0,
ifelse(n::even, a(n/2) + HammingWeight(n/2), a(n-1) + 1)) end:
seq(a(n), n = 0..78); # Peter Luschny, Oct 30 2021

a[n_] := (b = IntegerDigits[n, 2]).Reverse @ Range[Length @ b]; Array[a,78,0] (* Jean-François Alcover, Apr 28 2011, after B. Cloitre *)

for(n=0,100,l=length(binary(n)); print1(sum(i=1,l, component(binary(n),i)*(l-i+1)),","))

a(n) = my(b=binary(n)); b*-[-#b..-1]~; \\ Ruud H.G. van Tol, Oct 17 2023

def A029931(n): return sum(i if j == '1' else 0 for i, j in enumerate(bin(n)[:1:-1],1)) # Chai Wah Wu, Dec 20 2022
(C#)
ulong A029931(ulong n) {
    ulong result = 0, counter = 1;
    while(n > 0) {
        if (n % 2 == 1)
          result += counter;
        counter++;
        n /= 2;
    }
    return result;
} // Frank Hollstein, Jan 07 2023

a(n) = a(n - 2^L(n)) + L(n) + 1 [where L(n) = floor(log_2(n)) = A000523(n)] = sum of digits of A048794 [at least for n < 512]. - Henry Bottomley, Mar 09 2001
a(0) = 0, a(2n) = a(n) + e1(n), a(2n+1) = a(2n) + 1, where e1(n) = A000120(n). a(n) = log_2(A029930(n)). - Ralf Stephan, Jun 19 2003
G.f.: (1/(1-x)) * Sum_{k>=0} (k+1)*x^2^k/(1+x^2^k). - Ralf Stephan, Jun 23 2003
a(n) = Sum_{k>=0} A030308(n,k)*A000027(k+1). - Philippe Deléham, Oct 15 2011
a(n) = sum of n-th row of the triangle in A213629. - Reinhard Zumkeller, Jun 17 2012
From Reinhard Zumkeller, Feb 28 2014: (Start)
a(A089633(n)) = n and a(m) != n for m < A089633(n).
a(n) = Sum_{k=1..A070939(n)} k*A030308(n,k-1). (End)
a(n) = A073642(n) + A000120(n). - Peter Kagey, Apr 04 2016

More terms from Erich Friedman

A102627 Number of partitions of n into distinct parts in which the number of parts divides n.

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 4, 4, 5, 1, 15, 1, 7, 14, 17, 1, 28, 1, 40, 28, 11, 1, 99, 31, 13, 49, 99, 1, 186, 1, 152, 76, 17, 208, 425, 1, 19, 109, 699, 1, 584, 1, 433, 823, 23, 1, 1625, 437, 1140, 193, 746, 1, 2003, 1748, 2749, 244, 29, 1, 7404, 1, 31, 4158, 3258, 3766, 6307, 1

Offset: 1

Vladeta Jovovic, Feb 01 2005

			From _Gus Wiseman_, Sep 24 2019: (Start)
The a(1) = 1 through a(12) = 15 strict integer partitions whose average is an integer (A = 10, B = 11, C = 12):
  (1)  (2)  (3)  (4)   (5)  (6)    (7)  (8)   (9)    (A)   (B)  (C)
                 (31)       (42)        (53)  (432)  (64)       (75)
                            (51)        (62)  (531)  (73)       (84)
                            (321)       (71)  (621)  (82)       (93)
                                                     (91)       (A2)
                                                                (B1)
                                                                (543)
                                                                (642)
                                                                (651)
                                                                (732)
                                                                (741)
                                                                (831)
                                                                (921)
                                                                (5421)
                                                                (6321)
(End)

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

The BI-numbers of these partitions are given by A326669 (numbers whose binary indices have integer mean).
The non-strict case is A067538.
Strict partitions with integer geometric mean are A326625.
Strict partitions whose maximum divides their sum are A326850.
Cf. A018818, A033630, A316413, A326622, A326843, A326851.

a:= proc(m) option remember; local b; b:=
      proc(n, i, t) option remember; `if`(i*(i+1)/2Alois P. Heinz, Sep 25 2019

npdp[n_]:=Count[Select[IntegerPartitions[n],Length[#]==Length[ Union[ #]]&], ?(Divisible[n,Length[#]]&)]; Array[npdp,70] (* _Harvey P. Dale, Feb 12 2016 *)
a[m_] := a[m] = Module[{b}, b[n_, i_, t_] := b[n, i, t] = If[i(i+1)/2 < n, 0, If[n == 0, If[Mod[m, t] == 0, 1, 0], b[n, i - 1, t] + b[n - i, Min[n - i, i - 1], t + 1]]]; If[PrimeQ[m], 1, b[m, m, 0]]];
Array[a, 100] (* Jean-François Alcover, May 21 2021, after Alois P. Heinz *)

A359907 Number of strict integer partitions of n with integer median.

Original entry on oeis.org

0, 1, 1, 1, 2, 1, 4, 2, 6, 4, 9, 6, 14, 10, 18, 16, 27, 23, 36, 34, 51, 49, 67, 68, 94, 95, 122, 129, 166, 174, 217, 233, 287, 308, 371, 405, 487, 528, 622, 683, 805, 880, 1024, 1127, 1305, 1435, 1648, 1818, 2086, 2295, 2611, 2882, 3273, 3606, 4076, 4496, 5069

Offset: 0

Gus Wiseman, Jan 21 2023

nonn

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

			The a(1) = 1 through a(14) = 18 partitions (A..E = 10..14):
  1  2  3  4   5  6    7    8    9    A    B    C     D     E
           31     42   421  53   432  64   542  75    643   86
                  51        62   531  73   632  84    652   95
                  321       71   621  82   641  93    742   A4
                            431       91   731  A2    751   B3
                            521       532  821  B1    832   C2
                                      541       543   841   D1
                                      631       642   931   653
                                      721       651   A21   743
                                                732   6421  752
                                                741         761
                                                831         842
                                                921         851
                                                5421        932
                                                            941
                                                            A31
                                                            B21
                                                            7421

For mean instead of median: A102627, non-strict A067538 (ranked by A316413).
This is the strict case of A325347, ranked by A359908.
The median statistic is ranked by A360005(n)/2.
A000041 counts partitions, strict A000009.
A051293 counts subsets with integer mean, median A000975, cf. A005578.
A058398 counts partitions by mean, see also A008284, A327482.
A326567/A326568 gives the mean of prime indices.
A359893, A359901, A359902 count partitions by median.
Cf. A000016, A082550, A240219, A240850, A326669, A327475, A328966, A359897.

Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&IntegerQ[Median[#]]&]],{n,0,30}]

A291166 Connected Haar graph numbers.

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 33, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81

Offset: 1

Eric W. Weisstein, Aug 19 2017

nonn

Complement of A291165.

These appear to be numbers whose positions of 1's in their reversed binary expansion are relatively prime. If so, this sequence lists all positions of 1's in A326674. Numbers whose positions of 1's in their reversed binary expansion are pairwise coprime (as opposed to relatively prime) are A326675. - Gus Wiseman, Jul 19 2019

Eric Weisstein's World of Mathematics, Connected Graph
Eric Weisstein's World of Mathematics, Haar Graph

Cf. A000120, A029931, A048793, A059519, A070939, A289508, A291165, A295235, A319826, A326031, A326669, A326672, A326673.

A326675 The positions of 1's in the reversed binary expansion of n are pairwise coprime, where a singleton is not coprime unless it is {1}.

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 12, 13, 17, 18, 19, 20, 21, 22, 23, 24, 25, 28, 29, 33, 48, 49, 65, 66, 67, 68, 69, 70, 71, 72, 73, 76, 77, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 92, 93, 96, 97, 112, 113, 129, 132, 133, 144, 145, 148, 149, 192, 193, 196, 197, 208, 209, 212

Offset: 1

Gus Wiseman, Jul 17 2019

			41 has reversed binary expansion (1,0,0,1,0,1) with positions of 1's being {1,4,6}, which are not pairwise coprime, so 41 is not in the sequence.

Robert Israel, Table of n, a(n) for n = 1..10000

Equals the complement of A131577 in A087087.
Numbers whose prime indices are pairwise coprime are A302696.
Taking relatively prime instead of pairwise coprime gives A291166.
Cf. A000120, A051293, A070939, A289509, A291165, A319826, A326669, A326673, A326674.

extend:= proc(L) local C,c;
  C:= select(t -> andmap(s -> igcd(s,t)=1, L), [$1..L[-1]-1]);
  L, seq(procname([op(L),c]),c=C)
end proc:
g:= proc(L) local i;
  add(2^(i-1),i=L)
end proc:
map(g, [[1],seq(extend([k])[2..-1], k=2..10)]); # Robert Israel, Jul 19 2019

Select[Range[100],CoprimeQ@@Join@@Position[Reverse[IntegerDigits[#,2]],1]&]

is(n) = my (p=1); while (n, my (o=1+valuation(n,2)); if (gcd(p,o)>1, return (0), n-=2^(o-1); p*=o)); return (1) \\ Rémy Sigrist, Jul 19 2019

A326674 GCD of the set of positions of 1's in the reversed binary expansion of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 17 2019

a(n) is even if and only if n is in A062880. - Robert Israel, Oct 13 2020

			The reversed binary expansion of 40 is (0,0,0,1,0,1), with positions of 1's being {4,6}, so a(40) = GCD(4,6) = 2.

Robert Israel, Table of n, a(n) for n = 1..10000

Positions of 1's are A291166, and non-1's are A291165.
GCDs of prime indices are A289508.
GCDs of strict partitions encoded by FDH numbers are A319826.
Numbers whose binary positions are pairwise coprime are A326675.
Cf. A000120, A051293, A062880, A070939, A326667, A326668, A326669, A326670, A326672, A326673.

f:= proc(n) local B;
  B:= convert(n,base,2);
  igcd(op(select(t -> B[t]=1, [$1..ilog2(n)+1])))
end proc:
map(f, [$1..100]); # Robert Israel, Oct 13 2020

Table[GCD@@Join@@Position[Reverse[IntegerDigits[n,2]],1],{n,100}]

Trivially, a(n) <= log_2(n). - Charles R Greathouse IV, Nov 15 2022

A360241 Number of integer partitions of n whose distinct parts have integer mean.

Original entry on oeis.org

0, 1, 2, 2, 4, 3, 8, 6, 13, 13, 22, 19, 43, 34, 56, 66, 97, 92, 156, 143, 233, 256, 322, 341, 555, 542, 710, 831, 1098, 1131, 1644, 1660, 2275, 2484, 3035, 3492, 4731, 4848, 6063, 6893, 8943, 9378, 12222, 13025, 16520, 18748, 22048, 24405, 31446, 33698, 41558

Offset: 0

Gus Wiseman, Feb 02 2023

nonn

			The a(1) = 1 through a(8) = 13 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (311)    (33)      (331)      (44)
                    (31)    (11111)  (42)      (511)      (53)
                    (1111)           (51)      (3211)     (62)
                                     (222)     (31111)    (71)
                                     (321)     (1111111)  (422)
                                     (3111)               (2222)
                                     (111111)             (3221)
                                                          (3311)
                                                          (5111)
                                                          (32111)
                                                          (311111)
                                                          (11111111)
For example, the partition (32111) has distinct parts {1,2,3} with mean 2, so is counted under a(8).

For parts instead of distinct parts we have A067538, ranked by A316413.
The strict case is A102627.
These partitions are ranked by A326621.
For multiplicities instead of distinct parts: A360069, ranked by A067340.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by number of parts.
A051293 counts subsets with integer mean, median A000975.
A058398 counts partitions by mean, also A327482.
A116608 counts partitions by number of distinct parts.
A326619/A326620 gives mean of distinct prime indices.
A326622 counts factorizations with integer mean, strict A328966.
A360071 counts partitions by number of parts and number of distinct parts.
The following count partitions:
- A360242 mean(parts) != mean(distinct parts), ranked by A360246.
- A360243 mean(parts) = mean(distinct parts), ranked by A360247.
- A360250 mean(parts) > mean(distinct parts), ranked by A360252.
- A360251 mean(parts) < mean(distinct parts), ranked by A360253.
Cf. A067539, A078174, A082550, A240219, A316313, A325347, A326567/A326568, A326669, A327475, A349156.

Table[Length[Select[IntegerPartitions[n],IntegerQ[Mean[Union[#]]]&]],{n,0,30}]

A295235 Numbers k such that the positions of the ones in the binary representation of k are in arithmetic progression.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 17, 18, 20, 21, 24, 28, 30, 31, 32, 33, 34, 36, 40, 42, 48, 56, 60, 62, 63, 64, 65, 66, 68, 72, 73, 80, 84, 85, 96, 112, 120, 124, 126, 127, 128, 129, 130, 132, 136, 144, 146, 160, 168, 170, 192, 224, 240, 248

Offset: 1

Rémy Sigrist, Nov 18 2017

Also numbers k of the form Sum_{b=0..h-1} 2^(i+j*b) for some h >= 0, i >= 0, j > 0 (in fact, h = A000120(k), and if k > 0, i = A007814(k)).
There is a simple bijection between the finite sets of nonnegative integers in arithmetic progression and the terms of this sequence: s -> Sum_{i in s} 2^i; the term 0 corresponds to the empty set.
For any n > 0, A054519(n) gives the numbers of terms with n+1 digits in binary representation.
For any n >= 0, n is in the sequence iff 2*n is in the sequence.
For any n > 0, A000695(a(n)) is in the sequence.
The first prime numbers in the sequence are: 2, 3, 5, 7, 17, 31, 73, 127, 257, 8191, 65537, 131071, 262657, 524287, ...
This sequence contains the following sequences: A000051, A000079, A000225, A000668, A002450, A019434, A023001, A048645.
For any k > 0, 2^k - 2, 2^k - 1, 2^k, 2^k + 1 and 2^k + 2 are in the sequence (e.g., 14, 15, 16, 17, and 18).
Every odd term is a binary palindrome (and thus belongs to A006995).
Odd terms are A064896. - Robert Israel, Nov 20 2017

			The binary representation of the number 42 is "101010" and has ones evenly spaced, hence 42 appears in the sequence.
The first terms, alongside their binary representations, are:
   n  a(n)  a(n) in binary
  --  ----  --------------
   1    0           0
   2    1           1
   3    2          10
   4    3          11
   5    4         100
   6    5         101
   7    6         110
   8    7         111
   9    8        1000
  10    9        1001
  11   10        1010
  12   12        1100
  13   14        1110
  14   15        1111
  15   16       10000
  16   17       10001
  17   18       10010
  18   20       10100
  19   21       10101
  20   24       11000

Rémy Sigrist, Table of n, a(n) for n = 1..1000

Cf. A000051, A000079, A000120, A000225, A000668, A000695, A002450, A006995, A007814, A019434, A023001, A048645, A054519, A064896.

Cf. A029931, A048793 (binary indices triangle), A070939, A291166, A325328 (prime indices rather than binary indices), A326669, A326675.

f:= proc(d) local i,j,k;
  op(sort([seq(seq(add(2^(d-j*k),k=0..m),m=1..d/j),j=1..d),2^(d+1)]))
end proc:
0,1,seq(f(d),d=0..10); # Robert Israel, Nov 20 2017

bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
Select[Range[100],SameQ@@Differences[bpe[#]]&] (* Gus Wiseman, Jul 22 2019 *)

is(n) = my(h=hammingweight(n)); if(h<3, return(1), my(i=valuation(n,2),w=#binary(n)); if((w-i-1)%(h-1)==0, my(j=(w-i-1)/(h-1)); return(sum(k=0,h-1,2^(i+j*k))==n), return(0)))

A359402 Numbers whose binary expansion and reversed binary expansion have the same sum of positions of 1's, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

0, 1, 3, 5, 7, 9, 15, 17, 21, 27, 31, 33, 45, 51, 63, 65, 70, 73, 78, 85, 93, 99, 107, 119, 127, 129, 150, 153, 165, 189, 195, 219, 231, 255, 257, 266, 273, 282, 294, 297, 310, 313, 325, 334, 341, 350, 355, 365, 371, 381, 387, 397, 403, 413, 427, 443, 455, 471

Offset: 1

Gus Wiseman, Jan 05 2023

nonn

Also numbers whose binary expansion and reversed binary expansion have the same sum of partial sums.
Also numbers whose average position of a 1 in their binary expansion is (c+1)/2, where c is the number of digits.
Conjecture: Also numbers whose binary expansion has as least squares fit a line of zero slope, counted by A222955.

			The binary expansion of 70 is (1,0,0,0,1,1,0), with positions of 1's {1,5,6}, while the reverse positions are {2,3,7}. Both sum to 12, so 70 is in the sequence.

Binary words of this type appear to be counted by A222955.
For greater instead of equal sums we have A359401.
These are the indices of 0's in A359495.
A030190 gives binary expansion, reverse A030308.
A048793 lists partial sums of reversed standard compositions, sums A029931.
A070939 counts binary digits, 1's A000120.
A326669 lists numbers with integer mean position of a 1 in binary expansion.
Cf. A051293, A053632, A231204, A291166, A304818, A318283, A326672, A326673, A358134, A359042.

Select[Range[0,100],#==0||Mean[Join@@Position[IntegerDigits[#,2],1]]==(IntegerLength[#,2]+1)/2&]

from functools import reduce
from itertools import count, islice
def A359402_gen(startvalue=0): # generator of terms
    return filter(lambda n:(r:=reduce(lambda c, d:(c[0]+d[0]*(e:=int(d[1])),c[1]+e),enumerate(bin(n)[2:],start=1),(0,0)))[0]<<1==(n.bit_length()+1)*r[1],count(max(startvalue,0)))
A359402_list = list(islice(A359402_gen(),30)) # Chai Wah Wu, Jan 08 2023

A230877(a(n)) = A029931(a(n)).

cp's OEIS Frontend

A048793 List giving all subsets of natural numbers arranged in standard statistical (or Yates) order.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

C

Haskell

Maple

Mathematica

Formula

Extensions

A029931 If 2n = Sum 2^e_i, a(n) = Sum e_i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A102627 Number of partitions of n into distinct parts in which the number of parts divides n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A359907 Number of strict integer partitions of n with integer median.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A291166 Connected Haar graph numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

A326675 The positions of 1's in the reversed binary expansion of n are pairwise coprime, where a singleton is not coprime unless it is {1}.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A326674 GCD of the set of positions of 1's in the reversed binary expansion of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs