cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A373360 First differences of A373347.

Original entry on oeis.org

2, 2, 2, 2, 2, 1, 1, 2, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1
Offset: 1

Views

Author

Paolo Xausa, Jun 07 2024

Keywords

Comments

This sequence contains only 1's and 2's.

Crossrefs

Programs

  • Mathematica
    Differences[Select[Range[200], DigitSum[#, 2] > IntegerExponent[#, 2] + 1 &]]

Formula

a(n) = A373347(n+1) - A373347(n).

A371176 Numbers k such that A000120(k) <= A001511(k).

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 16, 18, 20, 24, 28, 32, 34, 36, 40, 44, 48, 52, 56, 64, 66, 68, 72, 76, 80, 84, 88, 96, 100, 104, 112, 120, 128, 130, 132, 136, 140, 144, 148, 152, 160, 164, 168, 176, 184, 192, 196, 200, 208, 216, 224, 232, 240, 256, 258, 260, 264, 268
Offset: 1

Views

Author

Mikhail Kurkov, Mar 14 2024

Keywords

Comments

It appears that this sequence is obtained when ordering Schreier sets as explained in the Bird link. See decM(n) PARI code. - Michel Marcus, May 31 2024
That is correct since the binary representation of these numbers can be put into 1-to-1 correspondence with Schreier sets, which satisfy |X| <= min X, using the indicator function of X as the bits (starting from the right, LSB). The reason is that A000120 then computes |X| and A001511 computes min X. For example, the Schreier set X = {2, 5} can be mapped to 10010_2 = 18. - Michael S. Branicky, May 31 2024
From David A. Corneth, May 31 2024: (Start)
If k is in the sequence then so is 2*k.
a(A000045(k)) = 2^(k-2) for k >= 2. (End)
Apart from a(1), all terms are even. - Paolo Xausa, May 31 2024
Zeckendorf representation of n with rewrite 0 -> 0, {0, 1} -> 1 and k-1 zeros appended to the right side (where k is the number of ones in the given representation) and then interpreted as binary expansion is the same as a(n) (see the first formula). - Mikhail Kurkov, Oct 21 2024

Crossrefs

Programs

  • Maple
    filter:= proc(n) convert(convert(n,base,2),`+`) <= 1+padic:-ordp(n,2) end proc:
    select(filter, [1,seq(i,i=2..1000,2)]); # Robert Israel, Oct 20 2024
  • Mathematica
    Join[{1}, Select[Range[2, 1000, 2], DigitSum[#, 2] <= IntegerExponent[#, 2] + 1 &]] (* Paolo Xausa, Aug 12 2025 *)
  • PARI
    isok(n) = hammingweight(n) <= (valuation(n, 2) + 1)
    
  • PARI
    M(n) = my(list=List()); for (i=1, n, forsubset(i, s, my(bOk = if (#s && (vecmax(s) == n), #s <= vecmin(s), 0)); if (bOk, listput(list, vecsort(Vec(s),,4))););); Vec(list);
    decM(nn) = my(v = vector(nn, k, M(k)), list=List()); for (i=1, #v, my(vi = v[i]); for (j=1, #vi, my(s = vecsort(vi[j]), slist=List(), m = vecmax(s)); forstep(k=m, 1, -1, listput(slist, sign(vecsearch(s, k)))); listput(list, fromdigits(Vec(slist), 2)););); vecsort(Vec(list)); \\ Michel Marcus, May 31 2024
    
  • Python
    def ok(n): return n.bit_count() <= (-n&n).bit_length()
    print([k for k in range(1, 300) if ok(k)]) # Michael S. Branicky, May 31 2024
    
  • Python
    # Assuming the list starts with 0.
    def a():
        n = na = nb = 1
        while True:
            yield not(nb < (na - 1) << 1)
            nb, na = na, n.bit_count()
            n += 1
    aList = a(); print([n for n in range(77) if next(aList)]) # Peter Luschny, Jun 07 2024

Formula

a(n) = b(n)*A001316(b(n))/2 where b(n) = A048679(n).
a(n) = Sum_{i=0..n-1} 2^A213911(i).
a(n) = 2^(A072649(n) - 1) + [c(n) > 0]*2*a(c(n)) where c(n) = A066628(n).
a(n) = 2*a(A005206(n)) - A003849(n)*2^A007895(n-1) for n > 1 with a(1) = 1.
Conjecture: lim sup_{n -> oo} log(a(n))/log(n) = log(2) / log((1 + sqrt(5))/2) = 1.440420090412556479... = A104287. - Vaclav Kotesovec, Aug 12 2025

A373345 Irregular triangle read by rows where row n lists (in decreasing order) the elements of the Schreier set encoded by A371176(n).

Original entry on oeis.org

1, 2, 3, 3, 2, 4, 4, 2, 4, 3, 5, 5, 2, 5, 3, 5, 4, 5, 4, 3, 6, 6, 2, 6, 3, 6, 4, 6, 4, 3, 6, 5, 6, 5, 3, 6, 5, 4, 7, 7, 2, 7, 3, 7, 4, 7, 4, 3, 7, 5, 7, 5, 3, 7, 5, 4, 7, 6, 7, 6, 3, 7, 6, 4, 7, 6, 5, 7, 6, 5, 4, 8, 8, 2, 8, 3, 8, 4, 8, 4, 3, 8, 5, 8, 5, 3, 8, 5, 4, 8, 6
Offset: 1

Views

Author

Paolo Xausa, Jun 01 2024

Keywords

Comments

A Schreier set is a subset of the positive integers with cardinality less than or equal to the minimum element in the set.
Each term k of A371176 can be put into a one-to-one correspondence with a Schreier set by interpreting the 1-based position of the ones in the binary expansion of k (where position 1 corresponds to the least significant bit) as the elements of the corresponding Schreier set (see A371176 and Bird link).
See A373359 for the elements in each set arranged in increasing order.
The number of sets having maximum element m is A000045(m).

Examples

			Triangle begins:
                                   Corresponding Schreier
   n  A371176(n)  bin(A371176(n))  set (this sequence)
  -------------------------------------------------------
   1      1              1         {1}
   2      2             10         {2}
   3      4            100         {3}
   4      6            110         {3, 2}
   5      8           1000         {4}
   6     10           1010         {4, 2}
   7     12           1100         {4, 3}
   8     16          10000         {5}
   9     18          10010         {5, 2}
  10     20          10100         {5, 3}
  11     24          11000         {5, 4}
  12     28          11100         {5, 4, 3}
  ...
		

Crossrefs

Cf. A007895 (conjectured row lengths), A072649 (first column), A373346 (row sums), A373347.

Programs

  • Mathematica
    Join[{{1}}, Map[Reverse[PositionIndex[Reverse[IntegerDigits[#, 2]]][1]] &, Select[Range[2, 200, 2], DigitCount[#, 2, 1] <= IntegerExponent[#, 2] + 1 &]]]

Formula

T(n,k) = A373557(n,k) - 1.

A373359 Irregular triangle read by rows where row n lists (in increasing order) the elements of the Schreier set encoded by A371176(n).

Original entry on oeis.org

1, 2, 3, 2, 3, 4, 2, 4, 3, 4, 5, 2, 5, 3, 5, 4, 5, 3, 4, 5, 6, 2, 6, 3, 6, 4, 6, 3, 4, 6, 5, 6, 3, 5, 6, 4, 5, 6, 7, 2, 7, 3, 7, 4, 7, 3, 4, 7, 5, 7, 3, 5, 7, 4, 5, 7, 6, 7, 3, 6, 7, 4, 6, 7, 5, 6, 7, 4, 5, 6, 7, 8, 2, 8, 3, 8, 4, 8, 3, 4, 8, 5, 8, 3, 5, 8, 4, 5, 8, 6, 8
Offset: 1

Views

Author

Paolo Xausa, Jun 04 2024

Keywords

Comments

See A373345 (where elements in each set are listed in decreasing order) for more information.

Examples

			Triangle begins:
                                   Corresponding Schreier
   n  A371176(n)  bin(A371176(n))  set (this sequence)
  -------------------------------------------------------
   1      1              1         {1}
   2      2             10         {2}
   3      4            100         {3}
   4      6            110         {2, 3}
   5      8           1000         {4}
   6     10           1010         {2, 4}
   7     12           1100         {3, 4}
   8     16          10000         {5}
   9     18          10010         {2, 5}
  10     20          10100         {3, 5}
  11     24          11000         {4, 5}
  12     28          11100         {3, 4, 5}
  ...
		

Crossrefs

Cf. A007895 (conjectured row lengths), A373346 (row sums), A373347.

Programs

  • Mathematica
    Join[{{1}}, Map[PositionIndex[Reverse[IntegerDigits[#, 2]]][1] &, Select[Range[2, 200, 2], DigitCount[#, 2, 1] <= IntegerExponent[#, 2] + 1 &]]]

Formula

T(n,k) = A373579(n,k) - 1.
Showing 1-4 of 4 results.