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-10 of 31 results. Next

A331579 Position of first appearance of n in A124758 (products of compositions in standard order).

Original entry on oeis.org

1, 2, 4, 8, 16, 18, 64, 34, 36, 66, 1024, 68, 4096, 258, 132, 136, 65536, 146, 262144, 264, 516, 4098
Offset: 1

Views

Author

Gus Wiseman, Mar 20 2020

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n. The k-th composition in standard order (row k of 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.

Examples

			The list of terms together with the corresponding compositions begins:
       1: (1)
       2: (2)
       4: (3)
       8: (4)
      16: (5)
      18: (3,2)
      64: (7)
      34: (4,2)
      36: (3,3)
      66: (5,2)
    1024: (11)
      68: (4,3)
    4096: (13)
     258: (7,2)
     132: (5,3)
     136: (4,4)
   65536: (17)
     146: (3,3,2)
  262144: (19)
     264: (5,4)
		

Crossrefs

The product of prime indices is A003963.
The sum of binary indices is A029931.
The sum of prime indices is A056239.
Sums of compositions in standard order are A070939.
The product of binary indices is A096111.
All terms belong to A114994.
Products of compositions in standard order are A124758.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    q=Table[Times@@stc[n],{n,1000}];
    Table[Position[q,i][[1,1]],{i,First[Split[Union[q],#1+1==#2&]]}]

A333227 Numbers k such that the k-th composition in standard order is pairwise coprime, where a singleton is not coprime unless it is (1).

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 11, 12, 13, 14, 15, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 30, 31, 33, 35, 37, 38, 39, 41, 44, 47, 48, 49, 50, 51, 52, 55, 56, 57, 59, 60, 61, 62, 63, 65, 66, 67, 68, 71, 72, 75, 77, 78, 79, 80, 83, 89, 92, 95, 96, 97, 99, 101, 102, 103, 105
Offset: 1

Views

Author

Gus Wiseman, Mar 27 2020

Keywords

Comments

This is the definition used for CoprimeQ in Mathematica.
The k-th composition in standard order (row k of 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.

Examples

			The sequence together with the corresponding compositions begins:
   1: (1)          27: (1,2,1,1)      55: (1,2,1,1,1)
   3: (1,1)        28: (1,1,3)        56: (1,1,4)
   5: (2,1)        29: (1,1,2,1)      57: (1,1,3,1)
   6: (1,2)        30: (1,1,1,2)      59: (1,1,2,1,1)
   7: (1,1,1)      31: (1,1,1,1,1)    60: (1,1,1,3)
   9: (3,1)        33: (5,1)          61: (1,1,1,2,1)
  11: (2,1,1)      35: (4,1,1)        62: (1,1,1,1,2)
  12: (1,3)        37: (3,2,1)        63: (1,1,1,1,1,1)
  13: (1,2,1)      38: (3,1,2)        65: (6,1)
  14: (1,1,2)      39: (3,1,1,1)      66: (5,2)
  15: (1,1,1,1)    41: (2,3,1)        67: (5,1,1)
  17: (4,1)        44: (2,1,3)        68: (4,3)
  18: (3,2)        47: (2,1,1,1,1)    71: (4,1,1,1)
  19: (3,1,1)      48: (1,5)          72: (3,4)
  20: (2,3)        49: (1,4,1)        75: (3,2,1,1)
  23: (2,1,1,1)    50: (1,3,2)        77: (3,1,2,1)
  24: (1,4)        51: (1,3,1,1)      78: (3,1,1,2)
  25: (1,3,1)      52: (1,2,3)        79: (3,1,1,1,1)
		

Crossrefs

A different ranking of the same compositions is A326675.
Ignoring repeated parts gives A333228.
Let q(k) be the k-th composition in standard order:
- The terms of q(k) are row k of A066099.
- The sum of q(k) is A070939(k).
- The product of q(k) is A124758(k).
- q(k) has A124767(k) runs and A333381(k) anti-runs.
- The GCD of q(k) is A326674(k).
- The Heinz number of q(k) is A333219(k).
- The LCM of q(k) is A333226(k).
Coprime or singleton sets are ranked by A087087.
Strict compositions are ranked by A233564.
Constant compositions are ranked by A272919.
Relatively prime compositions appear to be ranked by A291166.
Normal compositions are ranked by A333217.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,120],CoprimeQ@@stc[#]&]

A333228 Numbers k such that the distinct parts of the k-th composition in standard order (A066099) are pairwise coprime, where a singleton is not considered coprime unless it is (1).

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, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80
Offset: 1

Views

Author

Gus Wiseman, May 28 2020

Keywords

Comments

First differs from A291166 in lacking 69, which corresponds to the composition (4,2,1).
We use the Mathematica definition for CoprimeQ, so a singleton is not considered coprime unless it is (1).
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.

Examples

			The sequence together with the corresponding compositions begins:
   1: (1)          21: (2,2,1)        39: (3,1,1,1)
   3: (1,1)        22: (2,1,2)        41: (2,3,1)
   5: (2,1)        23: (2,1,1,1)      43: (2,2,1,1)
   6: (1,2)        24: (1,4)          44: (2,1,3)
   7: (1,1,1)      25: (1,3,1)        45: (2,1,2,1)
   9: (3,1)        26: (1,2,2)        46: (2,1,1,2)
  11: (2,1,1)      27: (1,2,1,1)      47: (2,1,1,1,1)
  12: (1,3)        28: (1,1,3)        48: (1,5)
  13: (1,2,1)      29: (1,1,2,1)      49: (1,4,1)
  14: (1,1,2)      30: (1,1,1,2)      50: (1,3,2)
  15: (1,1,1,1)    31: (1,1,1,1,1)    51: (1,3,1,1)
  17: (4,1)        33: (5,1)          52: (1,2,3)
  18: (3,2)        35: (4,1,1)        53: (1,2,2,1)
  19: (3,1,1)      37: (3,2,1)        54: (1,2,1,2)
  20: (2,3)        38: (3,1,2)        55: (1,2,1,1,1)
		

Crossrefs

Pairwise coprime or singleton partitions are A051424.
Coprime or singleton sets are ranked by A087087.
The version for relatively prime instead of coprime appears to be A291166.
Numbers whose binary indices are pairwise coprime are A326675.
Coprime partitions are counted by A327516.
Not ignoring repeated parts gives A333227.
The complement is A335238.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Sum is A070939.
- Product is A124758.
- Reverse is A228351
- GCD is A326674.
- Heinz number is A333219.
- LCM is A333226.
- Number of distinct parts is A334028.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,120],CoprimeQ@@Union[stc[#]]&]

A333766 Maximum part of the n-th composition in standard order. a(0) = 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 05 2020

Keywords

Comments

One plus the longest run of 0's in the binary expansion of n.
A composition of n is a finite sequence of positive integers summing to n. The k-th composition in standard order (row k of 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.

Examples

			The 100th composition in standard order is (1,3,3), so a(100) = 3.
		

Crossrefs

Positions of ones are A000225.
Positions of terms <= 2 are A003754.
The version for prime indices is A061395.
Positions of terms > 1 are A062289.
Positions of first appearances are A131577.
The minimum part is given by A333768.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Compositions without 1's are A022340.
- Sum is A070939.
- Product is A124758.
- Runs are counted by A124767.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Runs-resistance is A333628.
- Weakly decreasing compositions are A114994.
- Weakly increasing compositions are A225620.
- Strictly decreasing compositions are A333255.
- Strictly increasing compositions are A333256.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Table[If[n==0,0,Max@@stc[n]],{n,0,100}]

Formula

For n > 0, a(n) = A087117(n) + 1.

A062880 Zero together with the numbers which can be written as a sum of distinct odd powers of 2.

Original entry on oeis.org

0, 2, 8, 10, 32, 34, 40, 42, 128, 130, 136, 138, 160, 162, 168, 170, 512, 514, 520, 522, 544, 546, 552, 554, 640, 642, 648, 650, 672, 674, 680, 682, 2048, 2050, 2056, 2058, 2080, 2082, 2088, 2090, 2176, 2178, 2184, 2186, 2208, 2210, 2216, 2218, 2560, 2562
Offset: 0

Views

Author

Antti Karttunen, Jun 26 2001

Keywords

Comments

Binary expansion of n does not contain 1-bits at even positions.
Integers whose base-4 representation consists of only 0's and 2's.
Every nonnegative even number is a unique sum of the form a(k)+2*a(l); moreover, this sequence is unique with such property. - Vladimir Shevelev, Nov 07 2008
Also numbers such that the digital sum base 2 and the digital sum base 4 are in a ratio of 2:4. - Michel Marcus, Sep 23 2013
From Gus Wiseman, Jun 10 2020: (Start)
Numbers k such that the k-th composition in standard order has all even parts. 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. For example, the sequence of all compositions into even parts begins:
0: () 520: (6,4) 2080: (6,6)
2: (2) 522: (6,2,2) 2082: (6,4,2)
8: (4) 544: (4,6) 2088: (6,2,4)
10: (2,2) 546: (4,4,2) 2090: (6,2,2,2)
32: (6) 552: (4,2,4) 2176: (4,8)
34: (4,2) 554: (4,2,2,2) 2178: (4,6,2)
40: (2,4) 640: (2,8) 2184: (4,4,4)
42: (2,2,2) 642: (2,6,2) 2186: (4,4,2,2)
128: (8) 648: (2,4,4) 2208: (4,2,6)
130: (6,2) 650: (2,4,2,2) 2210: (4,2,4,2)
136: (4,4) 672: (2,2,6) 2216: (4,2,2,4)
138: (4,2,2) 674: (2,2,4,2) 2218: (4,2,2,2,2)
160: (2,6) 680: (2,2,2,4) 2560: (2,10)
162: (2,4,2) 682: (2,2,2,2,2) 2562: (2,8,2)
168: (2,2,4) 2048: (12) 2568: (2,6,4)
170: (2,2,2,2) 2050: (10,2) 2570: (2,6,2,2)
512: (10) 2056: (8,4) 2592: (2,4,6)
514: (8,2) 2058: (8,2,2) 2594: (2,4,4,2)
(End)

Crossrefs

Cf. A000695.
Except for first term, n such that A063694(n) = 0. Binary expansion is given in A062033.
Interpreted as Zeckendorf expansion: A062879.
Central diagonal of arrays A163357 and A163359.
Even partitions are counted by A035363.
Numbers with an even number of 1's in binary expansion are A001969.
Numbers whose binary expansion has even length are A053754.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Compositions without even parts are A060142.
- Sum is A070939.
- Product is A124758.
- Strict compositions are A233564.
- Heinz number is A333219.
- Number of distinct parts is A334028.

Programs

  • C
    uint32_t a_next(uint32_t a_n) { return (a_n + 0x55555556) & 0xaaaaaaaa; } /* Falk Hüffner, Jan 22 2022 */
    
  • Haskell
    a062880 n = a062880_list !! n
    a062880_list = filter f [0..] where
       f 0 = True
       f x = (m == 0 || m == 2) && f x'  where (x', m) = divMod x 4
    -- Reinhard Zumkeller, Nov 20 2012
    
  • Maple
    [seq(a(j),j=0..100)]; a := n -> add((floor(n/(2^i)) mod 2)*(2^((2*i)+1)),i=0..floor_log_2(n+1));
  • Mathematica
    b[n_] := BitAnd[n, Sum[2^k, {k, 0, Log[2, n] // Floor, 2}]]; Select[Range[ 0, 10^4], b[#] == 0&] (* Jean-François Alcover, Feb 28 2016 *)
  • Python
    def A062880(n): return int(bin(n)[2:],4)<<1 # Chai Wah Wu, Aug 21 2023

Formula

a(n) = 2 * A000695(n). - Vladimir Shevelev, Nov 07 2008
From Robert Israel, Apr 10 2018: (Start)
a(2*n) = 4*a(n).
a(2*n+1) = 4*a(n)+2.
G.f. g(x) satisfies: g(x) = 4*(1+x)*g(x^2)+2*x/(1-x^2). (End)

A329369 Number of permutations of {1,2,...,m} with excedance set constructed by taking m-i (0 < i < m) if b(i-1) = 1 where b(k)b(k-1)...b(1)b(0) (0 <= k < m-1) is the binary expansion of n.

Original entry on oeis.org

1, 1, 3, 1, 7, 3, 7, 1, 15, 7, 17, 3, 31, 7, 15, 1, 31, 15, 37, 7, 69, 17, 37, 3, 115, 31, 69, 7, 115, 15, 31, 1, 63, 31, 77, 15, 145, 37, 81, 7, 245, 69, 155, 17, 261, 37, 77, 3, 391, 115, 261, 31, 445, 69, 145, 7, 675, 115, 245, 15, 391, 31, 63, 1, 127, 63
Offset: 0

Views

Author

Mikhail Kurkov, Nov 12 2019

Keywords

Comments

Another version of A152884.
The excedance set of a permutation p of {1,2,...,m} is the set of indices i such that p(i) > i; it is a subset of {1,2,...,m-1}.
Great work on this subject was done by R. Ehrenborg and E. Steingrimsson, so most of the formulas given below are just their results translated into the language of the sequences which are related to the binary expansion of n.
Conjecture 1: equivalently, number of open tours by a biased rook on a specific f(n) X 1 board, which ends on a white cell, where f(n) = A070941(n) = floor(log_2(2n)) + 1 and cells are colored white or black according to the binary representation of 2n. A cell is colored white if the binary digit is 0 and a cell is colored black if the binary digit is 1. A biased rook on a white cell moves only to the left and otherwise moves only to the right. - Mikhail Kurkov, May 18 2021
Conjecture 2: this sequence is an inverse modulo 2 binomial transform of A284005. - Mikhail Kurkov, Dec 15 2021

Examples

			a(1) = 1 because the 1st excedance set is {m-1} and the permutations of {1,2,...,m} with such excedance set are 21, 132, 1243, 12354 and so on, i.e., for a given m we always have 1 permutation.
a(2) = 3 because the 2nd excedance set is {m-2} and the permutations of {1,2,...,m} with such excedance set are 213, 312, 321, 1324, 1423, 1432, 12435, 12534, 12543 and so on, i.e., for a given m we always have 3 permutations.
a(3) = 1 because the 3rd excedance set is {m-2, m-1} and the permutations of {1,2,...,m} with such excedance set are 231, 1342, 12453 and so on, i.e., for a given m we always have 1 permutation.
		

Crossrefs

Programs

  • Maple
    g:= proc(n) option remember;  2^padic[ordp](n, 2) end:
    a:= proc(n) option remember; `if`(n=0, 1, (h-> a(h)+
         `if`(n::odd, 0, (t-> a(h-t)+a(n-t))(g(h))))(iquo(n, 2)))
        end:
    seq(a(n), n=0..100);  # Alois P. Heinz, Jan 30 2023
  • Mathematica
    a[n_] := a[n] = Which[n == 0, 1, OddQ[n], a[(n-1)/2], True, a[n/2] + a[n/2 - 2^IntegerExponent[n/2, 2]] + a[n - 2^IntegerExponent[n/2, 2]]];
    a /@ Range[0, 65] (* Jean-François Alcover, Feb 13 2020 *)
  • PARI
    upto(n) = my(A, v1); v1 = vector(n+1, i, 0); v1[1] = 1; for(i=1, n, v1[i+1] = v1[i\2+1] + if(i%2, 0, A = 1 << valuation(i/2, 2); v1[i/2-A+1] + v1[i-A+1])); v1 \\ Mikhail Kurkov, Jun 06 2024

Formula

a(2n+1) = a(n) for n >= 0.
a(2n) = a(n) + a(n - 2^f(n)) + a(2n - 2^f(n)) for n > 0 with a(0) = 1 where f(n) = A007814(n) (equivalent to proposition 2.1 at the page 286, see R. Ehrenborg and E. Steingrimsson link).
a(2^m*(2n+1)) = Sum_{k=0..m} binomial(m+1,k) a(2^k*n) = a(2^m*n) + a(2^(m-1)*(2n+1)) + a(2^(m-1)*(4n+1)) for m > 0, n >= 0 (equivalent to proposition 2.5 at the page 287, see R. Ehrenborg and E. Steingrimsson link).
a(2n) = a(2*g(n)) + a(2n - 2^h(n)) + a(2*g(n) + 2^h(n)) for n > 0 with a(0) = 1 where g(n) = A053645(n), h(n) = A063250(n) (equivalent to proposition 2.1 at the page 286, see R. Ehrenborg and E. Steingrimsson link).
a(2n) = 2*a(n + g(n)) + a(2*g(n)) for n > 0, floor(n/3) < 2^(floor(log_2(n))-1) (in other words, for 2^m + k where 0 <= k < 2^(m-1), m > 0) with a(0) = 1 (just a special case of the previous formula, because for 2^m + k where 0 <= k < 2^(m-1), m > 0 we have 2^h(n) = n - g(n)).
a(2n) = a(f(n,-1)) + a(f(n,0)) + a(f(n,1)) for n > 0 with a(0) = 1 where f(n,k) = 2*(f(floor(n/2),k) + n mod 2) + k*A036987(n) for n > 1 with f(1,k) = abs(k) (equivalent to a(2n) = a(2*g(n)) + a(2n - 2^h(n)) + a(2*g(n) + 2^h(n))).
a(n) = Sum_{j=0..2^wt(n) - 1} (-1)^(wt(n) - wt(j)) Product_{k=0..wt(n) - 1} (1 + wt(floor(j/2^k)))^T(n,k) for n > 0 with a(0) = 1 where wt(n) = A000120(n), T(n,k) = T(floor(n/2), k - n mod 2) for k > 0 with T(n,0) = A001511(n) (equivalent to theorem 6.3 at page 296, see R. Ehrenborg and E. Steingrimsson link). Here T(n, k) - 1 for k > 0 is the length of the run of zeros between k-th pair of ones from the right side in the binary expansion of n. Conjecture 1: this formula is equivalent to inverse modulo 2 binomial transform of A284005.
Sum_{k=0..2^n-1} a(k) = (n+1)! for n >= 0.
a((4^n-1)/3) = A110501(n+1) for n >= 0.
a(2^2*(2^n-1)) = A091344(n+1),
a(2^3*(2^n-1)) = A091347(n+1),
a(2^4*(2^n-1)) = A091348(n+1).
More generally, a(2^m*(2^n-1)) = a(2^n*(2^m-1)) = S(n+1,m) for n >= 0, m >= 0 where S(n,m) = Sum_{k=1..n} k!*k^m*Stirling2(n,k)*(-1)^(n-k) (equivalent to proposition 6.5 at the page 297, see R. Ehrenborg and E. Steingrimsson link).
Conjecture 2: a(n) = (1 + A023416(n))*a(g(n)) + Sum_{k=0..floor(log_2(n))-1} (1-R(n,k))*a(g(n) + 2^k*(1 - R(n,k))) for n > 1 with a(0) = 1, a(1) = 1, where g(n) = A053645(n) and where R(n,k) = floor(n/2^k) mod 2 (at this moment this is the only formula here, which is not related to R. Ehrenborg's and E. Steingrimsson's work and arises from another definition given above, exactly conjectured definition with a biased rook). Here R(n,k) is the (k+1)-th bit from the right side in the binary expansion of n. - Mikhail Kurkov, Jun 23 2021
From Mikhail Kurkov, Jan 23 2023: (Start)
The formulas below are not related to R. Ehrenborg's and E. Steingrimsson's work.
Conjecture 3: a(n) = A357990(n, 1) for n >= 0.
Conjecture 4: a(2^m*(2k+1)) = Sum_{i=1..wt(k) + 2} i!*i^m*A358612(k, i)*(-1)^(wt(k) - i) for m >= 0, k >= 0 where wt(n) = A000120(n).
Conjecture 5: a(2^m*(2^n - 2^p - 1)) = Sum_{i=1..n} i!*i^m*(-1)^(n - i)*((i - p + 1)*Stirling2(n, i) - Stirling2(n - p, i - p) + Sum_{j=0..p-2} (p - j - 1)*Stirling2(n - p, i - j)/j! Sum_{k=0..j} (i - k)^p*binomial(j, k)*(-1)^k) for n > 2, m >= 0, 0 < p < n - 1. Here we consider that Stirling2(n, k) = 0 for n >= 0, k < 0. (End)
Conjecture 6: a(2^m*n + q) = Sum_{i=A001511(n+1)..A000120(n)+1} A373183(n, i)*a(2^m*(2^(i-1)-1) + q) for n >= 0, m >= 0, q >= 0. Note that this formula is recursive for n != 2^k - 1. Also, it is not related to R. Ehrenborg's and E. Steingrimsson's work. - Mikhail Kurkov, Jun 05 2024
From Mikhail Kurkov, Jul 10 2024: (Start)
a(2^m*(2^n*(2k+1) - 1)) = Sum_{i=1..m+1} a(2^i*k)*(-1)^(m-i+1)*Sum_{j=i..m+1} j^n*Stirling1(j, i)*Stirling2(m+1, j) for m >= 0, n >= 0, k >= 0 with a(0) = 1.
Proof: start with a(2^m*(2n+1)) = Sum_{k=0..m} binomial(m+1,k) a(2^k*n) given above and rewrite it as a(2^m*(2^n*(2k+1) - 1)) = Sum_{i=0..m} binomial(m+1, i) a(2^i*(2^(n-1)*(2k+1) - 1)).
Then conjecture that a(2^m*(2^n*(2k+1) - 1)) = Sum_{i=1..m+1} a(2^i*k)*f(n, m, i). From that it is obvious that f(0, m, i) = [i = (m+1)].
After that use a(2^m*(2^n*(2k+1) - 1)) = Sum_{i=0..m} binomial(m+1, i) Sum_{j=1..i+1} a(2^j*k)*f(n-1, i, j) = Sum_{i=1..m+1} a(2^i*k) Sum_{j=i-1..m} binomial(m+1, j)*f(n-1, j, i). From that it is obvious that f(n, m, i) = Sum_{j=i-1..m} binomial(m+1, j)*f(n-1, j, i).
Finally, all we need is to show that basic conditions and recurrence for f(n, m, i) gives f(n, m, i) = (-1)^(m-i+1)*Sum_{j=i..m+1} j^n*Stirling1(j, i)*Stirling2(m+1, j) (see Max Alekseyev link).
a(2^m*(2k+1)) = a(2^(m-1)*k) + (m+1)*a(2^m*k) + Sum_{i=1..m-1} a(2^m*k + 2^i) for m > 0, k >= 0.
Proof: start with a(2^(m+1)*(2k+1)) = a(2^m*k) + (m+2)*a(2^(m+1)*k) + Sum_{i=1..m} a(2^(m+1)*k + 2^i).
Then use a(2^m*(4k+1)) = a(2^m*k) + (m+1)*a(2^(m+1)*k) + Sum_{i=1..m-1} a(2^(m+1)*k + 2^i).
From that we get a(2^(m+1)*(2k+1)) - a(2^m*k) - (m+2)*a(2^(m+1)*k) - a(2^(m+1)*k + 2^m) = a(2^m*(4k+1)) - a(2^m*k) - (m+1)*a(2^(m+1)*k).
Finally, a(2^(m+1)*(2k+1)) = a(2^(m+1)*k) + a(2^m*(2*k+1)) + a(2^m*(4k+1)) which agrees with the a(2^m*(2n+1)) = a(2^m*n) + a(2^(m-1)*(2n+1)) + a(2^(m-1)*(4n+1)) given above.
This formula can be considered as an alternative to a(2^m*(2n+1)) = Sum_{k=0..m} binomial(m+1,k) a(2^k*n). There are algorithms for both these formulas that allow you to calculate them without recursion. However, even though it is necessary to calculate binomial coefficients in the mentioned formula, the triple-looped algorithm for it still works faster (see Peter J. Taylor link).
It looks like you can also change v2 in the mentioned algorithm to vector with elements a(2^m*(2^(i+A007814(n+1)-1)-1) + q) to get a(2^m*n + q) instead of a(n). This may have common causes with formula that uses A373183 given above. (End)
From Mikhail Kurkov, Jan 27 2025: (Start)
The formulas below are not related to R. Ehrenborg's and E. Steingrimsson's work.
Conjecture 7: A008292(n+1,k+1) = Sum_{i=0..2^n-1} [A000120(i) = k]*a(i) for n >= 0, k >= 0.
Conjecture 8: a(2^m*(2^n*(2k+1)-1)) = Sum_{i=0..m} Sum_{j=0..m-i} Sum_{q=0..i} binomial(m-i,j)*(m-j+1)^n*a(2^(q+1)*k)*L(m,i,q)*(-1)^j for m >= 0, n > 0, k >= 0 where L(n,k,m) = W(n-m,k-m,m+1) for n > 0, 0 <= k < n, 0 <= m <= k and where W(n,k,m) = (k+m)*W(n-1,k,m) + (n-k)*W(n-1,k-1,m) + [m > 1]*W(n,k,m-1) for 0 <= k < n, m > 0 with W(0,0,m) = 1, W(n,k,m) = 0 for n < 0 or k < 0.
In particular, W(n, k, 1) = A173018(n, k), W(n, k, 2) = A062253(n, k), W(n, k, 3) = A062254(n, k) and W(n, k, 4) = A062255(n, k).
Conjecture 9: a(n) = b(n,wt(n)) for n >= 0 where b(2n+1,k) = b(n,k) + (wt(n)-k+2)*b(n,k-1), b(2n,k) = (wt(n)-k+1)*b(2n+1,k) for n > 0, k > 0 with b(n,0) = A341392(n) for n >= 0, b(0,k) = 0 for k > 0 and where wt(n) = A000120(n) (see A379817).
More generally, a(2^m*(2k+1)) = ((m+1)!)^2*b(k,wt(k)-m) - Sum_{j=1..m} Stirling1(m+2,j+1)*a(2^(j-1)*(2k+1)) for m >= 0, k >= 0. Here we also consider that b(n,k) = 0 for k < 0. (End)
Conjecture 10: if we change b(n,0) = A341392(n) given above to b(n,0) = A341392(n)*x^n, then nonzero terms of the resulting polynomials for b(n,wt(n)) form c(n,k) such that a(n) = Sum_{k=0..A080791(n)} c(n,k) for n >= 0 where c(n,k) = (Product_{i=0..k-1} (1 + 1/A000120(floor(n/2^(A000523(n)-i))))) * Sum_{j=max{0,k-A080791(n)+A080791(A053645(n))}..A080791(A053645(n))} c(A053645(n),j) for n > 0, k >= 0 with c(0,0) = 1, c(0,k) = 0 for k > 0. - Mikhail Kurkov, Jun 19 2025

A335235 Numbers k such that the k-th composition in standard order (A066099) is pairwise coprime, where a singleton is always considered coprime.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 30, 31, 32, 33, 35, 37, 38, 39, 41, 44, 47, 48, 49, 50, 51, 52, 55, 56, 57, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 71, 72, 75, 77, 78, 79, 80, 83, 89, 92, 95, 96, 97
Offset: 1

Views

Author

Gus Wiseman, May 28 2020

Keywords

Comments

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.

Examples

			The sequence together with the corresponding compositions begins:
   1: (1)          20: (2,3)          48: (1,5)
   2: (2)          23: (2,1,1,1)      49: (1,4,1)
   3: (1,1)        24: (1,4)          50: (1,3,2)
   4: (3)          25: (1,3,1)        51: (1,3,1,1)
   5: (2,1)        27: (1,2,1,1)      52: (1,2,3)
   6: (1,2)        28: (1,1,3)        55: (1,2,1,1,1)
   7: (1,1,1)      29: (1,1,2,1)      56: (1,1,4)
   8: (4)          30: (1,1,1,2)      57: (1,1,3,1)
   9: (3,1)        31: (1,1,1,1,1)    59: (1,1,2,1,1)
  11: (2,1,1)      32: (6)            60: (1,1,1,3)
  12: (1,3)        33: (5,1)          61: (1,1,1,2,1)
  13: (1,2,1)      35: (4,1,1)        62: (1,1,1,1,2)
  14: (1,1,2)      37: (3,2,1)        63: (1,1,1,1,1,1)
  15: (1,1,1,1)    38: (3,1,2)        64: (7)
  16: (5)          39: (3,1,1,1)      65: (6,1)
  17: (4,1)        41: (2,3,1)        66: (5,2)
  18: (3,2)        44: (2,1,3)        67: (5,1,1)
  19: (3,1,1)      47: (2,1,1,1,1)    68: (4,3)
		

Crossrefs

The version counting partitions is A051424, with strict case A007360.
The version for binary indices is A087087.
The version counting compositions is A101268.
The version for prime indices is A302569.
The case without singletons is A333227.
The complement is A335236.
Numbers whose binary indices are pairwise coprime are A326675.
Coprime partitions are counted by A327516.
All of the following pertain to compositions in standard order:
- Length is A000120.
- The parts are row k of A066099.
- Sum is A070939.
- Product is A124758.
- Reverse is A228351
- GCD is A326674.
- Heinz number is A333219.
- LCM is A333226.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,100],Length[stc[#]]==1||CoprimeQ@@stc[#]&]

A060142 Ordered set S defined by these rules: 0 is in S and if x is in S then 2x+1 and 4x are in S.

Original entry on oeis.org

0, 1, 3, 4, 7, 9, 12, 15, 16, 19, 25, 28, 31, 33, 36, 39, 48, 51, 57, 60, 63, 64, 67, 73, 76, 79, 97, 100, 103, 112, 115, 121, 124, 127, 129, 132, 135, 144, 147, 153, 156, 159, 192, 195, 201, 204, 207, 225, 228, 231, 240, 243, 249, 252, 255, 256, 259, 265, 268, 271
Offset: 0

Views

Author

Clark Kimberling, Mar 05 2001

Keywords

Comments

After expelling 0 and 1, the numbers 4x occupy same positions in S that 1 occupies in the infinite Fibonacci word (A003849).
a(A026351(n)) = A219608(n); a(A004957(n)) = 4 * a(n). - Reinhard Zumkeller, Nov 26 2012
Apart from the initial term, this lists the indices of the 1's in A086747. - N. J. A. Sloane, Dec 05 2019
From Gus Wiseman, Jun 10 2020: (Start)
Numbers k such that the k-th composition in standard order has all odd parts, or numbers k such that A124758(k) is odd. 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. For example, the sequence of all compositions into odd parts begins:
0: () 57: (1,1,3,1) 135: (5,1,1,1)
1: (1) 60: (1,1,1,3) 144: (3,5)
3: (1,1) 63: (1,1,1,1,1,1) 147: (3,3,1,1)
4: (3) 64: (7) 153: (3,1,3,1)
7: (1,1,1) 67: (5,1,1) 156: (3,1,1,3)
9: (3,1) 73: (3,3,1) 159: (3,1,1,1,1,1)
12: (1,3) 76: (3,1,3) 192: (1,7)
15: (1,1,1,1) 79: (3,1,1,1,1) 195: (1,5,1,1)
16: (5) 97: (1,5,1) 201: (1,3,3,1)
19: (3,1,1) 100: (1,3,3) 204: (1,3,1,3)
25: (1,3,1) 103: (1,3,1,1,1) 207: (1,3,1,1,1,1)
28: (1,1,3) 112: (1,1,5) 225: (1,1,5,1)
31: (1,1,1,1,1) 115: (1,1,3,1,1) 228: (1,1,3,3)
33: (5,1) 121: (1,1,1,3,1) 231: (1,1,3,1,1,1)
36: (3,3) 124: (1,1,1,1,3) 240: (1,1,1,5)
39: (3,1,1,1) 127: (1,1,1,1,1,1,1) 243: (1,1,1,3,1,1)
48: (1,5) 129: (7,1) 249: (1,1,1,1,3,1)
51: (1,3,1,1) 132: (5,3) 252: (1,1,1,1,1,3)
(End)
Numbers whose binary representation has the property that every run of consecutive 0's has even length. - Harry Richman, Jan 31 2024

Examples

			From _Harry Richman_, Jan 31 2024: (Start)
In the following, dots are used for zeros in the binary representation:
   n  binary(a(n))  a(n)
   0:    .......     0
   1:    ......1     1
   2:    .....11     3
   3:    ....1..     4
   4:    ....111     7
   5:    ...1..1     9
   6:    ...11..    12
   7:    ...1111    15
   8:    ..1....    16
   9:    ..1..11    19
  10:    ..11..1    25
  11:    ..111..    28
  12:    ..11111    31
  13:    .1....1    33
  14:    .1..1..    36
  15:    .1..111    39
  16:    .11....    48
  17:    .11..11    51
  18:    .111..1    57
  19:    .1111..    60
  20:    .111111    63
  21:    1......    64
  22:    1....11    67
(End)
		

Crossrefs

Cf. A003714 (no consecutive 1's in binary expansion).
Odd partitions are counted by A000009.
Numbers with an odd number of 1's in binary expansion are A000069.
Numbers whose binary expansion has odd length are A053738.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Compositions without odd parts are A062880.
- Sum is A070939.
- Product is A124758.
- Strict compositions are A233564.
- Heinz number is A333219.
- Number of distinct parts is A334028.

Programs

  • Haskell
    import Data.Set (singleton, deleteFindMin, insert)
    a060142 n = a060142_list !! n
    a060142_list = 0 : f (singleton 1) where
       f s = x : f (insert (4 * x) $ insert (2 * x + 1) s') where
           (x, s') = deleteFindMin s
    -- Reinhard Zumkeller, Nov 26 2012
    
  • Mathematica
    Take[Nest[Union[Flatten[# /. {{i_Integer -> i}, {i_Integer -> 2 i + 1}, {i_Integer -> 4 i}}]] &, {1}, 5], 32]  (* Or *)
    Select[Range[124], FreeQ[Length /@ Select[Split[IntegerDigits[#, 2]], First[#] == 0 &], ?OddQ] &] (* _Birkas Gyorgy, May 29 2012 *)
  • PARI
    is(n)=if(n<3, n<2, if(n%2,is(n\2),n%4==0 && is(n/4))) \\ Charles R Greathouse IV, Oct 21 2013

Extensions

Corrected by T. D. Noe, Nov 01 2006
Definition simplified by Charles R Greathouse IV, Oct 21 2013

A284005 a(0) = 1, and for n > 1, a(n) = (1 + A000120(n))*a(floor(n/2)); also a(n) = A000005(A283477(n)).

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480, 600, 720, 64, 96, 144, 192, 216, 288, 384, 480, 324, 432, 576, 720, 768, 960, 1200, 1440, 486, 648, 864, 1080, 1152, 1440, 1800, 2160, 1536, 1920, 2400, 2880, 3000
Offset: 0

Views

Author

Antti Karttunen, Mar 18 2017

Keywords

Crossrefs

Similar recurrences: A124758, A243499, A329369, A341392.

Programs

  • Mathematica
    Table[DivisorSigma[0, #] &@ Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]], {n, 0, 71}] (* Michael De Vlieger, Mar 18 2017 *)
  • PARI
    A284005(n) = numdiv(A283477(n)); \\ edited by Michel Marcus, May 01 2019, M. F. Hasler, Nov 10 2019
    
  • PARI
    a(n) = my(k=if(n,logint(n,2)),s=1); prod(i=0,k, s+=bittest(n,k-i)); \\ Kevin Ryde, Jan 20 2021
  • Scheme
    (define (A284005 n) (A000005 (A283477 n)))
    

Formula

a(n) = A000005(A283477(n)).
Conjecture: a(n) = 2*a(f(n)) + Sum_{k=0..floor(log_2(n))-1} a(f(n) + 2^k*(1 - T(n,k))) for n > 1 with a(0) = 1, a(1) = 2, f(n) = A053645(n), T(n,k) = floor(n/2^k) mod 2. - Mikhail Kurkov, Nov 10 2019
From Mikhail Kurkov, Aug 23 2021: (Start)
a(2n+1) = a(n) + a(2n) for n >= 0.
a(2n) = a(n) + a(2n - 2^A007814(n)) for n > 0 with a(0) = 1. (End)
Conjecture: a(n) = Sum_{k=0..n} (binomial(n, k) mod 2)*A329369(k). In other words, this sequence is modulo 2 binomial transform of A329369. - Mikhail Kurkov, Mar 10 2023
Conjecture: a(2^m*(2n+1)) = Sum_{k=0..m+1} binomial(m+1, k)*a(2^k*n) for m >= 0, n >= 0 with a(0) = 1. - Mikhail Kurkov, Apr 24 2023

Extensions

Made Mikhail Kurkov's Nov 10 2019 formula the new primary name of this sequence - Antti Karttunen, Dec 30 2020

A111528 Square table, read by antidiagonals, where the g.f. for row n+1 is generated by: x*R_{n+1}(x) = (1+n*x - 1/R_n(x))/(n+1) with R_0(x) = Sum_{n>=0} n!*x^n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 3, 6, 1, 1, 4, 13, 24, 1, 1, 5, 22, 71, 120, 1, 1, 6, 33, 148, 461, 720, 1, 1, 7, 46, 261, 1156, 3447, 5040, 1, 1, 8, 61, 416, 2361, 10192, 29093, 40320, 1, 1, 9, 78, 619, 4256, 23805, 99688, 273343, 362880, 1, 1, 10, 97, 876, 7045, 48096, 263313
Offset: 0

Views

Author

Paul D. Hanna, Aug 06 2005

Keywords

Examples

			Table begins:
  1, 1,  2,   6,   24,   120,    720,    5040,     40320, ...
  1, 1,  3,  13,   71,   461,   3447,   29093,    273343, ...
  1, 1,  4,  22,  148,  1156,  10192,   99688,   1069168, ...
  1, 1,  5,  33,  261,  2361,  23805,  263313,   3161781, ...
  1, 1,  6,  46,  416,  4256,  48096,  591536,   7840576, ...
  1, 1,  7,  61,  619,  7045,  87955, 1187845,  17192275, ...
  1, 1,  8,  78,  876, 10956, 149472, 2195208,  34398288, ...
  1, 1,  9,  97, 1193, 16241, 240057, 3804353,  64092553, ...
  1, 1, 10, 118, 1576, 23176, 368560, 6262768, 112784896, ...
Rows are generated by logarithms of factorial series:
log(1 + x + 2*x^2 + 6*x^3 + 24*x^4 + ... n!*x^n + ...) = x + (3/2)*x^2 + (13/3)*x^3 + (71/4)*x^4 + (461/5)*x^5 + ...
(1/2)*log(1 + 2*x + 6*x^2 + ... + ((n+1)!/1!)*x^n + ...) = x + (4/2)*x^2 + (22/3)*x^3 + (148/4)*x^4 + (1156/5)*x^5 + ...
(1/3)*log(1 + 3*x + 12*x^2 + 60*x^3 + ... + ((n+2)!/2!)*x^n + ...) = x + (5/2)*x^2 + (33/3)*x^3 + (261/4)*x^4 + (2361/5)*x^5 +...
G.f. of row n may be expressed by the continued fraction:
R_n(x) = 1/(1+n*x - (n+1)*x/(1+(n+1)*x - (n+2)*x/(1+(n+2)*x -...
or recursively by: R_n(x) = 1/(1+n*x - (n+1)*x*R_{n+1}(x)).
		

Crossrefs

Cf: A003319 (row 1), A111529 (row 2), A111530 (row 3), A111531 (row 4), A111532 (row 5), A111533 (row 6), A111534 (diagonal).
Similar recurrences: A124758, A243499, A284005, A329369, A341392.

Programs

  • Maple
    T := (n, k) -> coeff(series(hypergeom([n+1, 1], [], x)/hypergeom([n, 1], [], x), x, 21), x, k):
    #display as a sequence
    seq(seq(T(n-k, k), k = 0..n), n = 0..10);
    # display as a square array
    seq(print(seq(T(n, k), k = 0..10)), n = 0..10); # Peter Bala, Jul 16 2022
  • Mathematica
    T[n_, k_] := T[n, k] = Which[n < 0 || k < 0, 0, k == 0 || k == 1, 1, n == 0, k!, True, (T[n - 1, k + 1] - T[n - 1, k])/n - Sum[T[n, j]*T[n - 1, k - j], {j, 1, k - 1}]]; Table[T[n - k, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 18 2018 *)
  • PARI
    {T(n,k)=if(n<0||k<0,0,if(k==0||k==1,1,if(n==0,k!, (T(n-1,k+1)-T(n-1,k))/n-sum(j=1,k-1,T(n,j)*T(n-1,k-j)))))}
    for(n=0,10,for(k=0,10,print1(T(n,k),", ")); print(""))
    
  • PARI
    {T(n,k)=if(n<0||k<0,0,if(k==0,1,if(n==0,k!, k/n*polcoeff(log(sum(m=0,k,(n-1+m)!/(n-1)!*x^m)),k))))}
    for(n=0,10,for(k=0,10,print1(T(n,k),", ")); print(""))

Formula

T(n, 0) = 1, T(0, k) = k!, otherwise for n>=1 and k>=1:
T(n, k) = (T(n-1, k+1) - T(n-1, k))/n - Sum_{j=1..k-1} T(n, j)*T(n-1, k-j).
T(n, k) = (k/n)*[x^k] log(Sum_{m=0..k} (n-1+m)!/(n-1)!*x^m).
T(n, k) = Sum_{j = 0..k} A089949(k, j)*n^(k-j). - Philippe Deléham, Aug 08 2005
R_n(x) = -((n-1)!/n)/Sum_{i>=1} (i+n-2)!*x^i, n > 0. - Vladeta Jovovic, May 06 2006
G.f. of row R may be expressed by the continued fraction: W(0), where W(k) = 1 - x*(k+1)/( x*(k+1) - 1/(1 - x*(k+1+R)/( x*(k+1+R) - 1/W(k+1) ))). - Sergei N. Gladkovskii, Aug 26 2013
Conjecture: T(n, k) = b(2^(k-1) - 1, n) for k > 0 with T(n, 0) = 1 where b(n, m) = b(floor(n/2), m) + b(floor((2n - 2^A007814(n))/2), m) + m*b(A025480(n-1), m) for n > 0 with b(0, m) = 1. - Mikhail Kurkov, Dec 16 2021
From Peter Bala, Jul 11 2022: (Start)
O.g.f. for row n, n >= 1: R(n,x) = ( Sum_{k >= 0} (n+k)!/n!*x^k )/( Sum_{k >= 0} (n-1+k)!/(n-1)!*x^k ).
R(n,x)/(1 - n*x*R(n,x)) = Sum_{k >= 0} (n+k)!/n!*x^k.
For n >= 0, R(n,x) satisfies the Riccati equation x^2*d/dx(R(n,x)) + n*x*R(n,x)^2 - (1 + (n-1)*x)*R(n,x) + 1 = 0 with R(n,0) = 1.
Apply Stokes 1982 to find that for n >= 0, R(n,x) = 1/(1 - x/(1 - (n+1)*x/(1 - 2*x/(1 - (n+2)*x/(1 - 3*x/(1 - (n+3)*x/(1 - 4*x/(1 - (n+4)*x/(1 - ...))))))))), a continued fraction of Stieltjes type. (End)
Showing 1-10 of 31 results. Next