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

A325781 Heinz numbers of complete integer partitions.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 30, 32, 36, 40, 42, 48, 54, 56, 60, 64, 72, 80, 84, 90, 96, 100, 108, 112, 120, 126, 128, 132, 140, 144, 150, 160, 162, 168, 176, 180, 192, 198, 200, 210, 216, 220, 224, 234, 240, 252, 256, 260, 264, 270, 280, 288, 294, 300
Offset: 1

Views

Author

Gus Wiseman, May 21 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The sum of prime indices of n is A056239(n). A number is in this sequence iff its divisors have sums of prime indices covering an initial interval of nonnegative integers. For example, the divisors of 60 are {1, 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, 60}, with respective sums of prime indices {0, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 7}, so 60 is in the sequence.

Examples

			The sequence of terms together with their prime indices begins:
     1: {}
     2: {1}
     4: {1,1}
     6: {1,2}
     8: {1,1,1}
    12: {1,1,2}
    16: {1,1,1,1}
    18: {1,2,2}
    20: {1,1,3}
    24: {1,1,1,2}
    30: {1,2,3}
    32: {1,1,1,1,1}
    36: {1,1,2,2}
    40: {1,1,1,3}
    42: {1,2,4}
    48: {1,1,1,1,2}
    54: {1,2,2,2}
    56: {1,1,1,4}
    60: {1,1,2,3}
    64: {1,1,1,1,1,1}

Crossrefs

Cf. A002033, A056239, A103295, A112798, A126796, A188431, A299702, A325684, A325770, A325780, A325782.

Programs

  • Mathematica
    normQ[m_]:=Or[m=={},Union[m]==Range[Max[m]]];
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Select[Range[1000],normQ[hwt/@Rest[Divisors[#]]]&]

A188431 The number of n-full sets, F(n).

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 2, 2, 1, 2, 1, 2, 3, 4, 5, 7, 7, 8, 9, 11, 10, 13, 14, 17, 20, 25, 28, 34, 40, 46, 54, 62, 69, 80, 90, 102, 115, 131, 144, 167, 186, 213, 239, 273, 304, 349, 388, 441, 495, 563, 625, 710, 790, 890, 990, 1114, 1232, 1387, 1530, 1713, 1894, 2119, 2330, 2605, 2866, 3192, 3512, 3910, 4289, 4774, 5237, 5809, 6377, 7068, 7739
Offset: 0

Views

Author

Madjid Mirzavaziri, Mar 31 2011

Keywords

Comments

Let A be a set of positive integers. We say that A is n-full if (sum A)=[n] for a positive integer n, where (sum A) is the set of all positive integers which are a sum of distinct elements of A and [n]={1,2,...,n}. Then F(n) denotes the number of n-full sets.
Also the number of distinct and complete partitions of n, by definition, which are counted by A000009 and A126796. - George Beck, Nov 06 2017
An integer partition of n is complete (see also A325781) if every number from 0 to n is the sum of some submultiset of the parts. The Heinz numbers of these partitions are given by A325986. - Gus Wiseman, May 31 2019

Examples

			a(26) = 10, because there are 10 26-full sets: {1,2,4,5,6,8}, {1,2,3,5,7,8}, {1,2,3,5,6,9}, {1,2,3,4,7,9}, {1,2,3,4,6,10}, {1,2,3,4,5,11}, {1,2,4,8,11}, {1,2,4,7,12}, {1,2,4,6,13}, {1,2,3,7,13}.
G.f.: 1 = 1/(1+x) + 1*x/((1+x)*(1+x^2)) + 0*x^2/((1+x)*(1+x^2)*(1+x^3)) + 1*x^3/((1+x)*(1+x^2)*(1+x^3)*(1+x^4)) +...+ a(n)*x^n / Product_{k=1..n+1} (1+x^k) +...

Links

Crossrefs

Cf. A188429, A188430, A126796.
Cf. A002033, A103295, A108917, A276024, A325763, A325765, A325781, A325782, A325788, A325986, A325790, A325791.

Programs

  • Haskell
    import Data.MemoCombinators (memo2, integral, Memo)
a188431 n = a188431_list !! (n-1)
a188431_list = map
   (\x -> sum [fMemo x i | i <- [a188429 x .. a188430 x]]) [1..] where
   fMemo = memo2 integral integral f
   f _ 1 = 1
   f m i = sum [fMemo (m - i) j |
                j <- [a188429 (m - i) .. min (a188430 (m - i)) (i - 1)]]
-- Reinhard Zumkeller, Aug 06 2015
  • Maple
    sums:= proc(s) local i, m;
          m:= max(s[]);
         `if`(m<1, {}, {m, seq([i, i+m][], i=sums(s minus {m}))})
       end:
a:= proc(n) local b;
      b:= proc(i,s) local si;
            if i=1 then `if`(sums(s)={$1..n}, 1, 0)
          else si:= s union {i};
               b(i-1, s)+ `if`(max(sums(si)[])>n, 0, b(i-1, si))
            fi
          end; b(n, {1})
    end:
seq(a(n), n=1..40);  # Alois P. Heinz, Apr 03 2011
# second Maple program:
b:= proc(n, i) option remember; `if`(i*(i+1)/2n or i>n-i+1, 0, b(n-i, i-1))))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..80);  # Alois P. Heinz, May 20 2017
  • Mathematica
    Sums[s_] := Sums[s] = With[{m = Max[s]}, If[m < 1, {}, Union @ Flatten @ Join[{m}, Table[{i, i + m}, {i, Sums[s ~Complement~ {m}]}]]]];
a[n_] := Module[{b}, b[i_, s_] := b[i, s] = Module[{si}, If[i == 1, If[Sums[s] == Range[n], 1, 0], si = s ~Union~ {i}; b[i-1, s] + If[Max[ Sums[si]] > n, 0, b[i - 1, si]]]]; b[n, {1}]];
Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 80}] (* Jean-François Alcover, Apr 12 2017, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Union[Total/@Union[Subsets[#]]]==Range[0,n]&]],{n,30}] (* Gus Wiseman, May 31 2019 *)
  • PARI
    /* As coefficients in g.f. */
{a(n)=local(A=[1]); for(i=1, n+1, A=concat(A,0); A[#A]=polcoeff(1 - sum(m=1,#A,A[m]*x^m/prod(k=1, m, 1+x^k +x*O(x^#A) )), #A) ); A[n+1]}
for(n=0, 50, print1(a(n),", ")) /* Paul D. Hanna, Mar 06 2012 */

Formula

F(n) = Sum_(i=L(n) .. U(n), F(n,i)), where F(n,i) = Sum_(j=L(n-i) .. min(U(n-i),i-1), F(n-i,j) ) and L(n), U(n) are defined in A188429 and A188430, respectively.
G.f.: 1 = Sum_{n>=0} a(n)*x^n / Product_{k=1..n+1} (1+x^k), with a(0)=1. - Paul D. Hanna, Mar 08 2012
a(n) ~ c * exp(Pi*sqrt(n/3)) / n^(3/4), where c = 0.03316508... - Vaclav Kotesovec, Oct 21 2019

Extensions

More terms from Alois P. Heinz, Apr 03 2011
a(0)=1 prepended by Alois P. Heinz, May 20 2017

A325780 Heinz numbers of perfect integer partitions.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 18, 20, 32, 42, 54, 56, 64, 100, 128, 162, 176, 234, 256, 260, 294, 392, 416, 486, 500, 512, 798, 1024, 1026, 1064, 1088, 1458, 1936, 2048, 2058, 2300, 2432, 2500, 2744, 3042, 3380, 4096, 4374, 4698, 5104, 5408, 5888, 8192, 8658, 9620, 10878
Offset: 1

Views

Author

Gus Wiseman, May 21 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The sum of prime indices of n is A056239(n). A number is in this sequence iff all of its divisors have distinct sums of prime indices, and these sums cover an initial interval of nonnegative integers. For example, the divisors of 260 are {1, 2, 4, 5, 10, 13, 20, 26, 52, 65, 130, 260}, with respective sums of prime indices {0, 1, 2, 3, 4, 6, 5, 7, 8, 9, 10, 11}, so 260 is in the sequence.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}
      2: {1}
      4: {1,1}
      6: {1,2}
      8: {1,1,1}
     16: {1,1,1,1}
     18: {1,2,2}
     20: {1,1,3}
     32: {1,1,1,1,1}
     42: {1,2,4}
     54: {1,2,2,2}
     56: {1,1,1,4}
     64: {1,1,1,1,1,1}
    100: {1,1,3,3}
    128: {1,1,1,1,1,1,1}
    162: {1,2,2,2,2}
    176: {1,1,1,1,5}
    234: {1,2,2,6}
    256: {1,1,1,1,1,1,1,1}
    260: {1,1,3,6}

Crossrefs

Equals the sorted concatenation of the triangle A258119.
A subsequence of A299702 and A325781.
Cf. A002033, A056239, A103300, A108917, A112798, A126796, A188431, A325763, A325768, A325782.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Select[Range[1000],Sort[hwt/@Rest[Divisors[#]]]==Range[DivisorSigma[0,#]-1]&]

Formula

Intersection of A299702 (knapsack partitions) and A325781 (complete partitions).

A371452 Number of connected components of the prime indices of the binary indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 01 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The prime indices of binary indices of 281492156579880 are {{1,1},{1,2},{3,4},{4,4}}, with 2 connected components {{1,1},{1,2}} and {{3,4},{4,4}}, so a(281492156579880) = 2.

Crossrefs

Positions of first appearances are A080355, opposite A325782.
For prime indices of prime indices we have A305079, ones A305078.
For binary indices of binary indices we have A326753, ones A326749.
Positions of ones are A371291.
For binary indices of prime indices we have A371451, ones A325118.
A001187 counts connected graphs.
A007718 counts non-isomorphic connected multiset partitions.
A048143 counts connected antichains of sets.
A048793 lists binary indices, reverse A272020, length A000120, sum A029931.
A070939 gives length of binary expansion.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A326964 counts connected set-systems, covering A323818.
Cf. A000720, A019565, A087086, A096111, A325097, A326782, A368109, A371292, A371294, A371445, A371447.

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Subsets[Range[Length[s]],{2}], Length[Intersection@@s[[#]]]>0&]},If[c=={},s, csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[csm[prix/@bix[n]]],{n,100}]

A353743 Least number with run-sum trajectory of length k; a(0) = 1.

Original entry on oeis.org

1, 2, 4, 12, 84, 1596, 84588, 11081028, 3446199708, 2477817590052, 4011586678294188, 14726534696017964148, 120183249654202605411828, 2146833388573021140471483564, 83453854313999050793547980583372, 7011542477899258250521520684673165324
Offset: 0

Views

Author

Gus Wiseman, Jun 11 2022

Keywords

Comments

Every sequence can be uniquely split into a sequence of non-overlapping runs. For example, the runs of (2,2,1,1,1,3,2,2) are ((2,2),(1,1,1),(3),(2,2)), with sums (4,3,3,4). The run-sum trajectory is obtained by repeatedly taking the run-sum transformation (A353832, A353847) until a squarefree number is reached. For example, the trajectory 12 -> 9 -> 7 corresponds to the partitions (2,1,1) -> (2,2) -> (4).

Examples

			The terms together with their prime indices begin:
      1: {}
      2: {1}
      4: {1,1}
     12: {1,1,2}
     84: {1,1,2,4}
   1596: {1,1,2,4,8}
  84588: {1,1,2,4,8,16}

Crossrefs

The ordered version is A072639, for run-lengths A333629.
The version for run-lengths is A325278, firsts in A182850 or A323014.
The run-sum trajectory is the iteration of A353832.
The first length-k row of A353840 has index a(k).
Other sequences pertaining to this trajectory are A353841-A353846.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A300273 ranks collapsible partitions, counted by A275870.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353866 ranks rucksack partitions, counted by A353864.
Cf. A002033, A005117, A006939, A071625, A076954, A126796, A181819, A182857, A188431, A299702, A325780, A325781, A353834.

Programs

  • Mathematica
    Join[{1,2},Table[2*Product[Prime[2^k],{k,0,n}],{n,0,6}]]

Formula

a(n > 1) = 2 * Product_{k=0..n-2} prime(2^k).
a(n > 0) = 2 * A325782(n).

A325789 Number of perfect necklace compositions of n.

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, May 22 2019

Keywords

Comments

A necklace composition of n is a finite sequence of positive integers summing to n that is lexicographically minimal among all of its cyclic rotations. A circular subsequence is a sequence of consecutive terms where the last and first parts are also considered consecutive. A necklace composition of n is perfect if every positive integer from 1 to n is the sum of exactly one distinct circular subsequence.

Examples

			The a(1) = 1 , a(2) = 1, a(3) = 2, a(7) = 3, a(13) = 5, and a(31) = 11 perfect necklace compositions (A = 10, B = 11, C = 12, D = 13, E = 14):
  1  11  12   124      1264           12546D
         111  142      1327           1274C5
              1111111  1462           13278A
                       1723           13625E
                       1111111111111  15C472
                                      17324E
                                      1A8723
                                      1D6452
                                      1E4237
                                      1E5263
                                      1111111111111111111111111111111

Crossrefs

Cf. A002033, A008965, A103300, A108917, A126796, A325676, A325679, A325685, A325780, A325782, A325786, A325787, A325789.

Programs

  • Mathematica
    neckQ[q_]:=Array[OrderedQ[{q,RotateRight[q,#]}]&,Length[q]-1,1,And];
subalt[q_]:=Union[ReplaceList[q,{_,s__,_}:>{s}],DeleteCases[ReplaceList[q,{t___,,u___}:>{u,t}],{}]];
Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],neckQ[#]&&Sort[Total/@subalt[#]]==Range[n]&]],{n,10}]

Formula

For n > 1, a(n) = A325787(n) + 1.

A325986 Heinz numbers of complete strict integer partitions.

Original entry on oeis.org

1, 2, 6, 30, 42, 210, 330, 390, 462, 510, 546, 714, 798, 2310, 2730, 3570, 3990, 4290, 4830, 5610, 6006, 6090, 6270, 6510, 6630, 7410, 7590, 7854, 8778, 8970, 9282, 9570, 9690, 10230, 10374, 10626, 11310, 11730, 12090, 12210, 12558, 13398, 13566, 14322, 14430
Offset: 1

Views

Author

Gus Wiseman, May 30 2019

Keywords

Comments

Strict partitions are counted by A000009, while complete partitions are counted by A126796.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
An integer partition of n is complete (A126796, A325781) if every number from 0 to n is the sum of some submultiset of the parts.
The enumeration of these partitions by sum is given by A188431.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}
      2: {1}
      6: {1,2}
     30: {1,2,3}
     42: {1,2,4}
    210: {1,2,3,4}
    330: {1,2,3,5}
    390: {1,2,3,6}
    462: {1,2,4,5}
    510: {1,2,3,7}
    546: {1,2,4,6}
    714: {1,2,4,7}
    798: {1,2,4,8}
   2310: {1,2,3,4,5}
   2730: {1,2,3,4,6}
   3570: {1,2,3,4,7}
   3990: {1,2,3,4,8}
   4290: {1,2,3,5,6}
   4830: {1,2,3,4,9}
   5610: {1,2,3,5,7}

Crossrefs

Cf. A002033, A056239, A103295, A112798, A126796, A188431, A299702, A304793.
Cf. A325780, A325781, A325782, A325788, A325790, A325988.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p] k]];
Select[Range[1000],SquareFreeQ[#]&&Union[hwt/@Divisors[#]]==Range[0,hwt[#]]&]

Formula

Intersection of A005117 (strict partitions) and A325781 (complete partitions).

A326018 Heinz numbers of knapsack partitions such that no addition of one part up to the maximum is knapsack.

Original entry on oeis.org

1925, 12155, 20995, 23375, 37145
Offset: 1

Views

Author

Gus Wiseman, Jun 03 2019

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
An integer partition is knapsack if every submultiset has a different sum.
The enumeration of these partitions by sum is given by A326016.

Examples

			The sequence of terms together with their prime indices begins:
   1925: {3,3,4,5}
  12155: {3,5,6,7}
  20995: {3,6,7,8}
  23375: {3,3,3,5,7}
  37145: {3,7,8,9}

Crossrefs

Cf. A002033, A108917, A275972, A299702, A299729, A304793.
Cf. A325780, A325782, A325857, A325862, A325878, A325880, A326015, A326016.

Programs

  • Mathematica
    ksQ[y_]:=UnsameQ@@Total/@Union[Subsets[y]];
Select[Range[2,200],With[{phm=If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]},ksQ[phm]&&Select[Table[Sort[Append[phm,i]],{i,Max@@phm}],ksQ]=={}]&]

A371451 Number of connected components of the binary indices of the prime indices of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 1, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 01 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The binary indices of prime indices of 805 are {{1,2},{3},{1,4}}, with 2 connected components {{1,2},{1,4}} and {{3}}, so a(805) = 2.

Links

Crossrefs

For prime indices of prime indices we have A305079, ones A305078.
Positions of ones are A325118.
Positions of first appearances are A325782.
For prime indices of binary indices we have A371452, ones A371291.
For binary indices of binary indices we have A326753, ones A326749.
A001187 counts connected graphs.
A007718 counts non-isomorphic connected multiset partitions.
A048143 counts connected antichains of sets.
A048793 lists binary indices, reverse A272020, length A000120, sum A029931.
A070939 gives length of binary expansion.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A326964 counts connected set-systems, covering A323818.
Cf. A000720, A019565, A087086, A096111, A325097, A326782, A368109, A371292, A371294, A371445, A371447.

Programs

  • Mathematica
    csm[s_]:=With[{c=Select[Subsets[Range[Length[s]],{2}], Length[Intersection@@s[[#]]]>0&]},If[c=={},s, csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[csm[bix/@prix[n]]],{n,100}]
  • PARI
    zero_first_elem_and_bitmask_connected_elems(ys) = { my(cs = List([ys[1]]), i=1); ys[1] = 0; while(i<=#cs, for(j=2, #ys, if(ys[j]&&(0!=bitand(cs[i], ys[j])), listput(cs, ys[j]); ys[j] = 0)); i++); (ys); };
A371451(n) = if(1==n, 0, my(cs = apply(p -> primepi(p), factor(n)[, 1]~), s=0); while(#cs, cs = select(c -> c, zero_first_elem_and_bitmask_connected_elems(cs)); s++); (s)); \\ Antti Karttunen, Jan 29 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 29 2025

A325787 Number of perfect strict necklace compositions of n.

Original entry on oeis.org

1, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, May 22 2019

Keywords

Comments

A strict necklace composition of n is a finite sequence of distinct positive integers summing to n that is lexicographically minimal among all of its cyclic rotations. In other words, it is a strict composition of n starting with its least part. A circular subsequence is a sequence of consecutive terms where the last and first parts are also considered consecutive. A necklace composition of n is perfect if every positive integer from 1 to n is the sum of exactly one distinct circular subsequence. For example, the composition (1,2,6,4) is perfect because it has the following circular subsequences and sums:
1: (1)
2: (2)
3: (1,2)
4: (4)
5: (4,1)
6: (6)
7: (4,1,2)
8: (2,6)
9: (1,2,6)
10: (6,4)
11: (6,4,1)
12: (2,6,4)
13: (1,2,6,4)
a(n) > 0 iff n = A002061(k) = A004136(k) for some k. - Bert Dobbelaere, Nov 11 2020

Examples

			The a(1) = 1 through a(31) = 10 perfect strict necklace compositions (empty columns not shown):
  (1)  (1,2)  (1,2,4)  (1,2,6,4)  (1,3,10,2,5)  (1,10,8,7,2,3)
              (1,4,2)  (1,3,2,7)  (1,5,2,10,3)  (1,13,6,4,5,2)
                       (1,4,6,2)                (1,14,4,2,3,7)
                       (1,7,2,3)                (1,14,5,2,6,3)
                                                (1,2,5,4,6,13)
                                                (1,2,7,4,12,5)
                                                (1,3,2,7,8,10)
                                                (1,3,6,2,5,14)
                                                (1,5,12,4,7,2)
                                                (1,7,3,2,4,14)
From _Bert Dobbelaere_, Nov 11 2020: (Start)
Compositions matching nonzero terms from a(57) to a(273), up to symmetry.
a(57) = 12:
  (1,2,10,19,4,7,9,5)
  (1,3,5,11,2,12,17,6)
  (1,3,8,2,16,7,15,5)
  (1,4,2,10,18,3,11,8)
  (1,4,22,7,3,6,2,12)
  (1,6,12,4,21,3,2,8)
a(73) = 8:
  (1,2,4,8,16,5,18,9,10)
  (1,4,7,6,3,28,2,8,14)
  (1,6,4,24,13,3,2,12,8)
  (1,11,8,6,4,3,2,22,16)
a(91) = 12:
  (1,2,6,18,22,7,5,16,4,10)
  (1,3,9,11,6,8,2,5,28,18)
  (1,4,2,20,8,9,23,10,3,11)
  (1,4,3,10,2,9,14,16,6,26)
  (1,5,4,13,3,8,7,12,2,36)
  (1,6,9,11,29,4,8,2,3,18)
a(133) = 36:
  (1,2,9,8,14,4,43,7,6,10,5,24)
  (1,2,12,31,25,4,9,10,7,11,16,5)
  (1,2,14,4,37,7,8,27,5,6,13,9)
  (1,2,14,12,32,19,6,5,4,18,13,7)
  (1,3,8,9,5,19,23,16,13,2,28,6)
  (1,3,12,34,21,2,8,9,5,6,7,25)
  (1,3,23,24,6,22,10,11,18,2,5,8)
  (1,4,7,3,16,2,6,17,20,9,13,35)
  (1,4,16,3,15,10,12,14,17,33,2,6)
  (1,4,19,20,27,3,6,25,7,8,2,11)
  (1,4,20,3,40,10,9,2,15,16,6,7)
  (1,5,12,21,29,11,3,16,4,22,2,7)
  (1,7,13,12,3,11,5,18,4,2,48,9)
  (1,8,10,5,7,21,4,2,11,3,26,35)
  (1,14,3,2,4,7,21,8,25,10,12,26)
  (1,14,10,20,7,6,3,2,17,4,8,41)
  (1,15,5,3,25,2,7,4,6,12,14,39)
  (1,22,14,20,5,13,8,3,4,2,10,31)
a(183) = 40:
  (1,2,13,7,5,14,34,6,4,33,18,17,21,8)
  (1,2,21,17,11,5,9,4,26,6,47,15,12,7)
  (1,2,28,14,5,6,9,12,48,18,4,13,16,7)
  (1,3,5,6,25,32,23,10,18,2,17,7,22,12)
  (1,3,12,7,20,14,44,6,5,24,2,28,8,9)
  (1,3,24,6,12,14,11,55,7,2,8,5,16,19)
  (1,4,6,31,3,13,2,7,14,12,17,46,8,19)
  (1,4,8,52,3,25,18,2,9,24,6,10,7,14)
  (1,4,20,2,12,3,6,7,33,11,8,10,35,31)
  (1,5,2,24,15,29,14,21,13,4,33,3,9,10)
  (1,5,23,27,42,3,4,11,2,19,12,10,16,8)
  (1,6,8,22,4,5,33,21,3,20,32,16,2,10)
  (1,8,3,10,23,5,56,4,2,14,15,17,7,18)
  (1,8,21,45,6,7,11,17,3,2,10,4,23,25)
  (1,9,5,40,3,4,21,35,16,18,2,6,11,12)
  (1,9,14,26,4,2,11,5,3,12,27,34,7,28)
  (1,9,21,25,3,4,8,5,6,16,2,36,14,33)
  (1,10,22,34,27,12,3,4,2,14,24,5,8,17)
  (1,10,48,9,19,4,8,6,7,17,3,2,34,15)
  (1,12,48,6,2,38,3,22,7,10,11,5,4,14)
a(273) = 12:
  (1,2,4,8,16,32,27,26,11,9,45,13,10,29,5,17,18)
  (1,3,12,10,31,7,27,2,6,5,19,20,62,14,9,28,17)
  (1,7,3,15,33,5,24,68,2,14,6,17,4,9,19,12,34)
  (1,7,12,44,25,41,9,17,4,6,22,33,13,2,3,11,23)
  (1,7,31,2,11,3,9,36,17,4,22,6,18,72,5,10,19)
  (1,21,11,50,39,13,6,4,14,16,25,26,3,2,7,8,27)
(End)

Links

Crossrefs

Cf. A002033, A002061, A004136, A008965, A032153, A103300, A126796, A325680, A325682, A325780, A325782, A325786, A325788, A325789.

Programs

  • Mathematica
    neckQ[q_]:=Array[OrderedQ[{q,RotateRight[q,#]}]&,Length[q]-1,1,And];
subalt[q_]:=Union[ReplaceList[q,{_,s__,_}:>{s}],DeleteCases[ReplaceList[q,{t___,,u___}:>{u,t}],{}]];
Table[Length[Select[Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&],neckQ[#]&&Sort[Total/@subalt[#]]==Range[n]&]],{n,30}]

Extensions

More terms from Bert Dobbelaere, Nov 11 2020
Showing 1-10 of 13 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.