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 26 results. Next

A380697 Frobenius number of the set S = {e_i+2; 1 <= i <= m}, where the e_i are the exponents in the binary expansion n = Sum_{i=1..m} 2^e_i, or 0 if GCD(S) = A326674(2*n) > 1.

Original entry on oeis.org

0, 0, 1, 0, 0, 5, 1, 0, 3, 7, 1, 11, 3, 2, 1, 0, 0, 0, 1, 0, 0, 5, 1, 19, 3, 7, 1, 7, 3, 2, 1, 0, 5, 11, 1, 17, 5, 5, 1, 23, 3, 4, 1, 6, 3, 2, 1, 29, 5, 11, 1, 9, 5, 5, 1, 9, 3, 4, 1, 3, 3, 2, 1, 0, 0, 13, 1, 0, 0, 5, 1, 27, 3, 7, 1, 11, 3, 2, 1, 0, 0, 13, 1
Offset: 1

Views

Author

Pontus von Brömssen, Jan 30 2025

Keywords

Comments

The sequence gives the Frobenius numbers of all sets of integers greater than 1, encoded by the binary expansion of n.

Examples

			For n = 262288 = 2^4+2^7+2^18, a(n) is the Frobenius number of {6, 9, 20}, i.e., the last term of A065003, so a(262288) = 43.
		

Crossrefs

Formula

a(n) = 1 if and only if n == 3 (mod 4) (i.e., if and only if n is in A004767).
a(n) = 2 if and only if n == 14 (mod 16).
a(2^e+2^f) = (e+1)*(f+1)-1 for nonnegative integers e and f such that e+2 and f+2 are coprime.

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

			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)
		

References

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

Crossrefs

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.

Programs

  • C
    #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
    }
    
  • Haskell
    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
  • Maple
    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
  • Mathematica
    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 *)

Formula

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

Extensions

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

Views

Author

Keywords

Comments

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

Examples

			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
		

Crossrefs

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.
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.

Programs

  • Haskell
    a029931 = sum . zipWith (*) [1..] . a030308_row
    -- Reinhard Zumkeller, Feb 28 2014
    
  • Maple
    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
  • Mathematica
    a[n_] := (b = IntegerDigits[n, 2]).Reverse @ Range[Length @ b]; Array[a,78,0] (* Jean-François Alcover, Apr 28 2011, after B. Cloitre *)
  • PARI
    for(n=0,100,l=length(binary(n)); print1(sum(i=1,l, component(binary(n),i)*(l-i+1)),","))
    
  • PARI
    a(n) = my(b=binary(n)); b*-[-#b..-1]~; \\ Ruud H.G. van Tol, Oct 17 2023
    
  • Python
    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

Formula

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

Extensions

More terms from Erich Friedman

A228351 Triangle read by rows in which row n lists the compositions (ordered partitions) of n (see Comments lines for definition).

Original entry on oeis.org

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

Views

Author

Omar E. Pol, Aug 30 2013

Keywords

Comments

The representation of the compositions (for fixed n) is as lists of parts, the order between individual compositions (for the same n) is (list-)reversed co-lexicographic. - Joerg Arndt, Sep 02 2013
Dropping the "(list-)reversed" in the comment above gives A228525.
The equivalent sequence for partitions is A026792.
This sequence lists (without repetitions) all finite compositions, in such a way that, if [P_1, ..., P_r] denotes the composition occupying the n-th position in the list, then (((2*n/2^(P_1)-1)/2^(P_2)-1)/...)/2^(P_r)-1 = 0. - Lorenzo Sauras Altuzarra, Jan 22 2020
The k-th composition in the list is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, and taking first differences. Reversing again gives A066099, which is described as the standard ordering. Both sequences define a bijective correspondence between nonnegative integers and integer compositions. - Gus Wiseman, Apr 01 2020
It follows from the previous comment that A000120(k) is the length of the k-th composition that is listed by this sequence (recall that A000120(k) is the number of 1's in the binary expansion of k). - Lorenzo Sauras Altuzarra, Sep 29 2020

Examples

			Illustration of initial terms:
-----------------------------------
n  j     Diagram     Composition j
-----------------------------------
.         _
1  1     |_|         1;
.         _ _
2  1     |_  |       2,
2  2     |_|_|       1, 1;
.         _ _ _
3  1     |_    |     3,
3  2     |_|_  |     1, 2,
3  3     |_  | |     2, 1,
3  4     |_|_|_|     1, 1, 1;
.         _ _ _ _
4  1     |_      |   4,
4  2     |_|_    |   1, 3,
4  3     |_  |   |   2, 2,
4  4     |_|_|_  |   1, 1, 2,
4  5     |_    | |   3, 1,
4  6     |_|_  | |   1, 2, 1,
4  7     |_  | | |   2, 1, 1,
4  8     |_|_|_|_|   1, 1, 1, 1;
.
Triangle begins:
[1];
[2],[1,1];
[3],[1,2],[2,1],[1,1,1];
[4],[1,3],[2,2],[1,1,2],[3,1],[1,2,1],[2,1,1],[1,1,1,1];
[5],[1,4],[2,3],[1,1,3],[3,2],[1,2,2],[2,1,2],[1,1,1,2],[4,1],[1,3,1],[2,2,1],[1,1,2,1],[3,1,1],[1,2,1,1],[2,1,1,1],[1,1,1,1,1];
...
For example [1,2] occupies the 5th position in the corresponding list of compositions and indeed (2*5/2^1-1)/2^2-1 = 0. - _Lorenzo Sauras Altuzarra_, Jan 22 2020
12 --binary expansion--> [1,1,0,0] --reverse--> [0,0,1,1] --positions of 1's--> [3,4] --prepend 0--> [0,3,4] --first differences--> [3,1]. - _Lorenzo Sauras Altuzarra_, Sep 29 2020
		

Crossrefs

Row n has length A001792(n-1). Row sums give A001787, n >= 1.
Cf. A000120 (binary weight), A001511, A006519, A011782, A026792, A065120.
A related ranking of finite sets is A048793/A272020.
All of the following consider the k-th row to be the k-th composition, ignoring the coarser grouping by sum.
- Indices of weakly increasing rows are A114994.
- Indices of weakly decreasing rows are A225620.
- Indices of strictly decreasing rows are A333255.
- Indices of strictly increasing rows are A333256.
- Indices of reversed interval rows A164894.
- Indices of interval rows are A246534.
- Indices of strict rows are A233564.
- Indices of constant rows are A272919.
- Indices of anti-run rows are A333489.
- Row k has A124767(k) runs and A333381(k) anti-runs.
- Row k has GCD A326674(k) and LCM A333226(k).
- Row k has Heinz number A333219(k).
Equals A163510+1, termwise.
Cf. A124734 (increasing length, then lexicographic).
Cf. A296774 (increasing length, then reverse lexicographic).
Cf. A337243 (increasing length, then colexicographic).
Cf. A337259 (increasing length, then reverse colexicographic).
Cf. A296773 (decreasing length, then lexicographic).
Cf. A296772 (decreasing length, then reverse lexicographic).
Cf. A337260 (decreasing length, then colexicographic).
Cf. A108244 (decreasing length, then reverse colexicographic).
Cf. A228369 (lexicographic).
Cf. A066099 (reverse lexicographic).
Cf. A228525 (colexicographic).

Programs

  • Haskell
    a228351 n = a228351_list !! (n - 1)
    a228351_list = concatMap a228351_row [1..]
    a228351_row 0 = []
    a228351_row n = a001511 n : a228351_row (n `div` 2^(a001511 n))
    -- Peter Kagey, Jun 27 2016
    
  • Maple
    # Program computing the sequence:
    A228351 := proc(n) local c, k, L, N: L, N := [], [seq(2*r, r = 1 .. n)]: for k in N do c := 0: while k != 0 do if gcd(k, 2) = 2 then k := k/2: c := c+1: else L := [op(L), op(c)]: k := k-1: c := 0: fi: od: od: L[n]: end: # Lorenzo Sauras Altuzarra, Jan 22 2020
    # Program computing the list of compositions:
    List := proc(n) local c, k, L, M, N: L, M, N := [], [], [seq(2*r, r = 1 .. 2^n-1)]: for k in N do c := 0: while k != 0 do if gcd(k, 2) = 2 then k := k/2: c := c+1: else L := [op(L), c]: k := k-1: c := 0: fi: od: M := [op(M), L]: L := []: od: M: end: # Lorenzo Sauras Altuzarra, Jan 22 2020
  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Table[Differences[Prepend[bpe[n],0]],{n,0,30}] (* Gus Wiseman, Apr 01 2020 *)
  • Python
    from itertools import count, islice
    def A228351_gen(): # generator of terms
        for n in count(1):
            k = n
            while k:
                yield (s:=(~k&k-1).bit_length()+1)
                k >>= s
    A228351_list = list(islice(A228351_gen(),30)) # Chai Wah Wu, Jul 17 2023

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

Views

Author

Eric W. Weisstein, Aug 19 2017

Keywords

Comments

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

Crossrefs

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[#]&]

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

Views

Author

Gus Wiseman, Jul 17 2019

Keywords

Examples

			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.
		

Crossrefs

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.

Programs

  • Maple
    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
  • Mathematica
    Select[Range[100],CoprimeQ@@Join@@Position[Reverse[IntegerDigits[#,2]],1]&]
  • PARI
    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

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[#]]&]

A124758 Product of the parts of the compositions in standard order.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 2, 1, 4, 3, 4, 2, 3, 2, 2, 1, 5, 4, 6, 3, 6, 4, 4, 2, 4, 3, 4, 2, 3, 2, 2, 1, 6, 5, 8, 4, 9, 6, 6, 3, 8, 6, 8, 4, 6, 4, 4, 2, 5, 4, 6, 3, 6, 4, 4, 2, 4, 3, 4, 2, 3, 2, 2, 1, 7, 6, 10, 5, 12, 8, 8, 4, 12, 9, 12, 6, 9, 6, 6, 3, 10, 8, 12, 6, 12, 8, 8, 4, 8, 6, 8, 4, 6, 4, 4, 2, 6, 5, 8, 4, 9, 6
Offset: 0

Views

Author

Keywords

Comments

The standard order of compositions is given by A066099.
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. - Gus Wiseman, Apr 03 2020

Examples

			Composition number 11 is 2,1,1; 2*1*1 = 2, so a(11) = 2.
The table starts:
  1
  1
  2 1
  3 2 2 1
  4 3 4 2 3 2 2 1
  5 4 6 3 6 4 4 2 4 3 4 2 3 2 2 1
The 146-th composition in standard order is (3,3,2), with product 18, so a(146) = 18. - _Gus Wiseman_, Apr 03 2020
		

Crossrefs

Cf. A066099, A118851, A011782 (row lengths), A001906 (row sums).
The lengths of standard compositions are given by A000120.
The version for prime indices is A003963.
The version for binary indices is A096111.
Taking the sum instead of product gives A070939.
The sum of binary indices is A029931.
The sum of prime indices is A056239.
Taking GCD instead of product gives A326674.
Positions of first appearances are A331579.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Table[Times@@stc[n],{n,0,100}] (* Gus Wiseman, Apr 03 2020 *)

Formula

For a composition b(1),...,b(k), a(n) = Product_{i=1}^k b(i).
a(A164894(n)) = a(A246534(n)) = n!. - Gus Wiseman, Apr 03 2020
a(A233249(n)) = a(A333220(n)) = A003963(n). - Gus Wiseman, Apr 03 2020
From Mikhail Kurkov, Jul 11 2021: (Start)
Some conjectures:
a(2n+1) = a(n) for n >= 0.
a(2n) = (1 + 1/A001511(n))*a(n) = 2*a(n) + a(n - 2^f(n)) - a(2n - 2^f(n)) for n > 0 with a(0)=1 where f(n) = A007814(n).
From the 1st formula for a(2n) we get a(4n+2) = 2*a(n), a(4n) = 2*a(2n) - a(n).
Sum_{k=0..2^n - 1} a(k) = A001519(n+1) for n >= 0.
a((4^n - 1)/3) = A011782(n) for n >= 0.
a(2^m*(2^n - 1)) = m + 1 for n > 0, m >= 0. (End)

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[#]&]
Showing 1-10 of 26 results. Next